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TORONTO MAIN STP MISA PILOT SITE STUDY 
COMPONENT REPORT 
WATER QUALITY 


EXECUTIVE SUMMARY 


As a component of the Ministry’s "Municipal-Industrial Strategy for Abatement" (MISA) 
program, the Ministry has evaluated several water quality assessment techniques at six "Pilot 
Sites" in Ontario. The site-specific water quality-based effluent limits consider the 
environmental sensitivity of receiving water bodies, and evaluate the potential need for 
effluent limits which are more stringent than those being developed by the technology track 


investigations. 


The purpose of this report is to present and discuss the results of water quality surveys 
conducted in 1986 and 1987 in support of the Toronto Main STP MISA pilot site study. 
Several model predictions of effluent reductions required to achieve Provincial Water 
Quality Objectives (PWQO) for parameters of concern are compared. It should be noted 


that this report is a technical report that is not intended to set final policy. 


Findings: 


Dilution factors for various chemicals calculated for surface samples in the vicinity of the 
Main STP outfall are given in Table 1. The average initial dilution at the outfall station 
(1419) is approximately 1:3 for conservative parameters and 1:4 to 1:7 for parameters 
subject to decay. Within 1 km to the south of the outfall, the dilution is approximately 1:20 
to 1:50. Exceedences of PWQO in the vicinity of the Main STP outfall are summarized in 
Table 2. 


In the immediate vicinity of the outfall, only the ammonia and residual chlorine 
concentrations approach levels which are acutely toxic to the biota. For ammonia, the 
threshold for acute toxicity depends on both temperature and pH, but is approximately 10 


mg/L total ammonia, for the temperature and pH observed during the study (Emerson, et 


al., 1975). The average outfall concentration is 7.0 mg/L, and the PWQO is exceeded 75% 
of the time. A slight pH depression at the outfall slightly favours ionization of ammonia, 


thereby suppressing toxicity to some extent. 


For chlorine, the acute toxic threshold for freshwater fish is approximately 50 to 200 ug/L 
(U.S. EPA, 1984). The average outfall concentration over the period of STP chlorination 
(June through September) is 330 ug/L, and the Provincial Water Quality Objective 
(PWQO) (2 ug/L) is exceeded 70% of the time. Based on observed water quality data, 


acute toxic effects would not be expected beyond a 500 m radius of the outfall. 


Total phosphorus concentrations, average 0.2 mg/L near the outfall, almost always 
exceeding the provincial water quality guideline (0.02 mg/L), and are still above the 
background level of 0.015 mg/L at a distance of 3 km south of the outfall. These levels are 


of concern primarily because of potential stimulation of nuisance algal growth. 


Iron, copper and cadmium exceed their provincial receiving water objectives at the outfall 
77, 38 and 31% of the time, respectively. The average concentration of iron (0.456 mg/L) 
exceeds the provincial objective of 0.3 mg/L. Other heavy metals do not exceed objectives 
on average. Provincial water quality objectives are designed to provide protection to 
aquatic life from chronic toxicity. Acute toxic levels are not approached by any heavy 


metal. 


Fecal coliform bacteria, suspended solids and conductivity are of interest primarily as 
surrogates for human pathogens or toxic contaminants that may be either adsorbed to 
particulates or dissolved. All these parameters approach background levels within 1 km of 


the outfall, but are two to three times higher than background at the outfall. 


The organochlorine pesticides lindane and heptachlor have exceeded provincial water 
quality objectives at the outfall. Only lindane is frequently in exceedence (17% of the time 
in 1986, 46% in 1987). However, even at the maximum concentrations observed, 
Structure-activity models suggest that chronic toxic effects from these chemicals on aquatic 
biota are almost certainly absent. Average lindane concentrations at the outfall (10 ng/L 


in 1987) are equal to the provincial water quality objective, and are an order of magnitude 


below the chronic toxicity threshold. 


Ultra-trace methods (16-L effluent and 200-L lake water samples) were used to measure 
the concentrations of PCBs, organochlorine pesticides and chlorobenzenes in the STP 
effluent and receiving water. Within the lake, few contaminants are present at 
concentrations above 1 ng/L. Such contaminants include the dichlorobenzenes, a-BHC, 
7-BHC (lindane) and, in October 1987, total PCBs in 23 of 30 samples collected. This 
finding is most significant for total PCBs for which the PWQO is 1 ng/L. Several other 
chlorinated benzenes were present at concentrations between 0.1 and 1 ng/L. Other 


contaminants averaged below 0.1 ng/L. 


In the STP effluent, the various chlorinated benzenes and PCBs were always found using 
ultra-trace methods. Mean loadings for these compounds were as high as 630 g/day 
(1,4-dichlorobenzene); the mean loadings for total PCBs and lindane were 17 and 21 g/day, 
respectively. Average effluent concentrations of these three compounds were 710, 20 and 
24 ng/L, respectively. However, these loadings are minor compared with those from 
external sources, in particular the atmosphere and Niagara River. In addition to loading 
reductions from Toronto sources, PCB loadings from other sources discharging to Lake 
Ontario must be reduced considerably in order to achieve compliance in the Toronto 
Waterfront. 


The size and direction of the effluent plume may easily be established by tracking using 
parameters which are easily measured on board the survey vessel, such as conductivity, 
ammonia and chlorine residual. The magnitude of the effluent plume is strongly influenced 
by the lake currents. Weak currents (<5 cm/s) result in rapid dilution within several 
hundred metres of the outfall. On the other hand, strong currents (> 15 cm/s) produce 
poor dilution, and the effects of the plume can be observed for as much as 6 km from the 
outfall, even in the vicinity of the R.C. Harris WIP intake. However, as the plume is at 
the surface, increased concentrations are unlikely to be found at the depth of the WIP 


intake. 


Historical data suggest that conductivity, turbidity and total ammonia have been generally 
decreasing on the Toronto Waterfront in the recent past; however, the Don River mouth 


and Main STP outfall have shown opposite trends since 1983. Total phosphorus has been 
generally increasing on the Toronto Waterfront, especially in the vicinity of the Don River 
mouth and Main STP. Thus, these two sources are further implicated by long-term 


temporal trends as major influences on local water quality. 


Total ammonia in the Main STP effluent shows a consistent pattern of increase over the 
course of the day from about 10:00 a.m. to 5:00 p.m. Phosphorus is more tightly controlled, 
though extremes tend to occur late in the day. Flow may either increase or decrease over 


the course of the day. Both ammonia and flow exhibit large day-to-day variations. 


Long-term prediction of daily average effluent outfall concentrations may be possible with 
some success based on daily average concentrations in the Main STP effluent. However, 
short-term prediction based on instantaneous effluent measurements requires a 
hydrodynamic model which explicitly accounts for receiving water conditions. Existing data 


are inadequate for development of these predictive relationships. 


Derivation of Effluent Limits 


For parameters of concern, the allowable load is the load for which the defined objective 
(in most cases the PWQO) will be satisfied for a specific fraction of the time at the edge 
of any defined mixing zone. For the purpose of this exercise, two scenarios were modelled: 
a mixing zone of 500 m radius about the outfall, and no mixing zone (i.e. zero radius - 
achieve PWQO at the edge of near-field diffuser mixing). Also for the purposes of this 
exercise, the water quality-based effluent limits were derived on the basis of 95% 
compliance with objectives - i.e. no parameter will be allowed to be out of compliance more 
than 5% of the time, given the statistical variability of contaminant loadings. Parameters 


of concern include those which frequently exceed objectives near the outfall. 


Initial dilution factors above the outfall (i.e. station 1419 - edge of near-field mixing) have 
been defined from average water quality data (Table 1). Using only average to worst-case 
current conditions, the mean far-field dilution rate was found from plume-tracking 
experiments to be 9%/100 m for conductivity and ammonia, and 19%/100m for chlorine 
residual. The higher rate for chlorine residual reflects its rapid chemical and biochemical 


decay. These rates were combined with loading data to estimate the extent of reduction 
required in order to meet the stated objectives or guidelines at zero or 500 m radius from 
the outfall (Tables 3, 4). 


In addition, two models were used to predict allowable loadings at the present and 
proposed outfall sites. The dilution model combines the effects of diffuser geometry, water 
depth and mixing depth to estimate initial dilution, plus receiving water hydrodynamics to 
estimate far-field dilution. The fate model uses modelled lake hydrodynamics for the 
Toronto Waterfront nearshore area plus processes such as volatilization, sorption, bio- 
uptake and sedimentation to estimate in-lake concentrations at various distances from the 
outfall, from assumed loading values. However, it suffers from the limitation that 
concentrations are calculated on the basis of 500 m grid squares; thus the "zero" mixing 
zone result is actually an average of up to 250 m from the outfall. Similarly for the 500 m 
mixing zone, it represents an average of 250-750 m. In addition, modelled hydrodynamics 
used in the fate model near the edge of the modelled area (i.e. proposed outfall) are known 
to be less accurate. For these reasons, only the dilution model results are presented in this 


summary with the plume track results. 


The dilution model was run for average conditions at the existing outfall location, plus a 
proposed outfall with a 1000 m diffuser located about 3.7 km from shore in a water depth 


of at least 15 m. Results of this model are also summarized in Tables 3 and 4. 


For a zero mixing zone at the existing location (Table 3), required loading reductions range 
from 50-70% (iron, lindane) to 96-98% for total phosphorus and greater than 99% for 
chlorine residual. At 500 m mixing zone radius (Table 4), the required reductions are only 
slightly lower for most parameters due to slow far-field dilution which occurs under most 
current conditions. However, if the proposed outfall is adopted, objectives should be met 
for all parameters except total phosphorus, chlorine residual and (for zero mixing zone) 
ammonia. This is due to the improved diffuser design as well as greater initial dilution 


resulting from the improved water depth. 


It must be kept in mind that these loading reductions are estimates which are subject to 


inherent variability due to short-term changes in loadings or receiving water physical 


conditions which are hard to define. However, they represent a range from average to 
worst-case conditions in order to provide realistic estimates of required dilutions for 
conditions which prevail at both the existing and proposed locations for the largest part of 


the time. 


Conclusions and Recommendations: 


The extent of loading reductions required to meet Provincial Water Quality Objectives at 
various locations, has been estimated using the plume tracking results and dilution model 
explained above. It should be noted that these are only estimates of required reductions; 
the final model for loading reduction calculations has not been selected, nor has the 


allowable distance from the outfall been finalized. 
The main conclusions and recommendations of this study are as follows: 


i: It is important that site-specific water quality-based effluent limits be developed for 
individual dischargers. This is evident due to the difference in effluent reductions 


required for the present and proposed outfalls. 


2 Dechlorination of the final effluent will be necessary in order to meet Provincial 
Water Quality Objectives (PWQO) in the receiving water. 


° At the existing outfall location, in order to meet PWQO at the end of the 
pipe 95% of the time, residual chlorine levels in the final effluent must be 
reduced by 99.6% to 5 ug/L. In order to meet PWQO at the edge of a 500m 
mixing zone 95% of the time, final effluent chlorine levels must be reduced 
by 99% to 16 ug/L. 


e At the proposed offshore outfall location, in order to meet PWOQO at the end 
of the pipe 95% of the time, chlorine levels in the effluent must be reduced 
by 97% to 47 ug/L. In order to meet PWQO at the edge of a 500m mixing 
zone 95% of the time, chlorine levels in the effluent must be reduced by 94% 
to 86 ug/L. 


Ammonia removal in the final effluent will be necessary to meet PWQO in the 


receiving water under summer conditions. | 


e At the existing outfall, in order to meet PWQO at the end of the pipe 95% 
of the time, ammonia levels in the final effluent must be reduced by 88-93% 
to 2.3-4.3 mg/L. Similarly, for the edge of a 500m mixing zone, the levels 
must be reduced by 79-87% to 4.3-6.9 mg/L. 


e At the proposed offshore location, in order to meet PWQO at the end of the 
pipe 95% of the time, ammonia levels in the final effluent must be reduced 
by 28% to 23.5 mg/L. However, no reduction in effluent ammonia levels is 


required to meet PWQO 95% of the time at the edge of a 500m mixing zone. 


Enhanced phosphorus removal (tertiary treatment) of the final effluent will be 
necessary to meet the Provincial Water Quality Guidelines (PWQG) for prevention 
of nuisance algal growth. 


@ At the existing outfall location, in order to meet the PWQG at the end of the 
pipe 95% of the time, final effluent total phosphorus levels must be reduced 
by 96-98% to 24-46 pg/L. Similarly for the edge of a 500m mixing zone, the 

_ level must be reduced by 95% to 68 ug/L. Based on dilution model 
predictions, distances of 17 and 35 km are required for dilution to PWQG of 
effluent concentrations of 0.5 and 1.0 mg/L, respectively. 


e At the proposed offshore outfall location, in order to meet PWQG 95% of 
the time at the end of the pipe, final effluent total phosphorus levels must be 
reduced by 80% to 245 yg/L. Similarly for a 500m mixing zone, the levels 
must be reduced by 65% to 440 ug/L. Based on dilution model predictions, 
distances of 0.7 and 3.6 km are required for dilution to PWQG of effluent 
concentrations of 0.5 and 1.0 mg/L, respectively. 


> At the existing outfall location, loading reductions of heavy metals and organic 
contaminants will be necessary in order to meet PWQO in the receiving water. 
However, improved dilution characteristics at the proposed offshore outfall location 


will allow PWQO for these parameters to be achieved after initial dilution. 


e In order to meet PWQO for copper at the end of the pipe 95% of the time 
at the existing outfall location, the final effluent level must be reduced by 65- 
75% to 10-15 yg/L. Similarly, for the edge of a 500m mixing zone, the 
required reduction is 43-58% to 18-24 pg/L. Required reductions for other 
metals and organics such as lindane are lower and should be achieved if the 


above reductions are met. 


e In the case of PCBs, loading reductions of 93 to 95% are estimated for zero 
and 500 m mixing zones, based on average observed background levels. 
However, the highest background PCB levels are similar to the PWQO of 1 
ng/L; hence required loading reductions may be as high as 100% for either 


outfall location. 


Recommendations for Further Study 


A number of recommendations for design of future pilot studies were suggested. In 
general, the study design should be derived from the study objectives. When there are 
multiple study objectives, each sampling activity should be related to one of these objectives 
so that sampling design can be optimized. 


While short-term (i.e. daily cycle) studies may be valuable for modelling purposes, study 
objectives should address the long-term question of objective compliance. For these 
purposes, daily 24-hour composite sampling, supported by more frequent (hourly to 3-h) 
grab sampling, should be conducted for one or two weeks per survey. Based on the results 
of the short-term survey, effort should be made to limit temporal sampling bias, both 
between days and between regions. For these purposes, a specific sampling design should 
be adopted, including a series of fixed stations and depths. 


In addition, further refinement of mixing zone models for adequate description of near- 
field dilution (initial mixing) is required. Such refined models could be used to run 
conditions of existing treatment at the new outfall location to determine the size of mixing 


zone required for various parameters of concern, such as total phosphorus and ammonia. 


When large amounts of "below detection" data are obtained, careful consideration to 
handling such data is required. Statistical methods such as the maximum likelihood method 
provide the best estimate of overall conditions; however, the best solution to this problem 


is improved detection limits in the laboratory for critical substances such as PCBs. 


TABLE 1: DILUTION FACTORS’ FOR DIFFERENT CHEMICALS AT SELECTED 
INDIVIDUAL STATIONS IN 1987 


Station Number’ 


Chemical 
Parameter? 1419 PHT) 2909 2878 2910 2209 

Field Ammonia: 

0.5 m depth sic ul SW) 32% 14. 45. 80. 

6.0 m depth 80. 300. 300. 45. 300. = 
Lab Ammonia 2.9 40. 110. 8.7 19. 110. 
Residual Chlorine’ 3.0 216 28. - - 400. 
Phosphorus 5.0 PNSIE 40. PAT) = 40. 80. 
Conductivity 370 18. Zan 14. 40. 100. 
DOC 333 20. 27. 16. 50. 130. 
Fecal Coliforms 6.6 470. 700. 27. 50. 700. 
Copper 5.4 36. 100. Zon 29. 150. 
Zinc 7.1 55. 150. 75. 70. 400. 
Iron 4.1 20. 28. 16. 40. 30. 
Manganese 55 74 aes 28. 12. 24. 75. 

* Dilution Factor = (Source-Background) / (Concentration-Background) 


Chemical measured at 0.5 m water depth in receiving water, unless 
otherwise indicated. 


Source concentration of free chlorine assumed = 1000 ug/L; other 
concentrations averaged over period of Main STP chlorination 
(June-September). 


Station locations are shown on Figure 1.2. 


TABLE 2: EXCEEDENCE OF CONVENTIONAL AND METAL WATER QUALITY 
OBJECTIVES AT STATION 1419 OVER THE MAIN 8TP OUTFALL 


(see Appendix 2 for complete statistical summary) 


1986 1987 
Parameter N % N % PWQO’ 
0.5 m Depth 
Residual Chlorine*** 7 100 13 70 2.0 pg/L 
Field Ammonia 6 100 12 72 0.02* 
Lab Ammonia 6 100 - - 0.02* 
Phosphorus 10 100 14 100 0.02 
Iron 10 40 13 77 0.30 
Copper 10 30 13 38 0.005 
Cadmium 10 20 13 31 0.2 po/L 
Lead = me 13 0 0.025** 
Nickel 10 0 13 0 0.025 
Mercury 2 11 13 8 0.2 pg/L 
6 m Depth 
Residual Chlorine*** 6 . 100 5 0 2.0 pg/L 
Field Ammonia 6 67 8 0 0.02* 
Lab Ammonia 6 50 - - 0.02* 
Phosphorus 10 80 4 100 0.02 
Iron 9 22 4 25 0.30 
Copper 9 11 4 25 0.005 
Cadmium 9 11 4 25 0.2 ug/L 
Lead = = 4 0 0.025** 
Nickel 9 0 4 25 0.025 
Mercury 8 0 4 0 0.2 pg/L 


Note: N = Number of samples taken. 


PWQO in mg/L unless otherwise indicated; phosphorus value is a 
guideline only. 

* Ammonia objective is for unionized fraction. 

** Lead objective is alkalinity-dependent; alkalinity = 80 mg/L 
was assumed. 

**k*k Chlorination period (to end of September) only. 
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1.0 INTRODUCTION 


Under the Municipal and Industrial Strategy for Abatement (MISA) being developed by 
Ontario, all dischargers in the province will be required to update their pollution control 
technology at least to the level of BATEA (Best Available Technology Economically 
Achievable). If receiving water quality conditions so dictate, the dischargers may be 
required to achieve even higher levels of pollution control. Regulations to govern the 
allowable discharge from various industrial and municipal sectors are being developed by 
MOE based on BATEA,; intensive studies are being conducted at six "pilot sites" situated 
in various parts of the province to examine effluent requirements based on receiving water 
conditions. The Toronto Main Sewage Treatment Plant (STP) was selected as a pilot site, 
as it is one of the largest municipal sewage treatment plants discharging to the Great Lakes 
and receives influents from a large variety of sources (i.e., industrial, domestic, storm 


overflow). 


Given that protection of the environment is the ultimate goal of the MISA initiative, the 
MOE additionally undertook to examine site-specific, water-quality-based effluent limits at 
the six pilot sites, to protect environmentally sensitive areas, to develop assessment 
protocols and to eventually assess those effluent requirements relative to those prescribed 
by BATEA. Water quality based effluent limits can be derived through the use of 
sophisticated modelling techniques which establish a predictive link between effluent quality 
and specific effects on the aquatic environment. Additionally, the results of the Main STP 
study are expected to contribute significantly to the development of the Toronto Waterfront 
Remedial Action Plan (RAP). The RAP is being prepared for the International Joint 
Commission (IJC) by a combined federal-provincial work team, in cooperation with the 
municipalities and a formal public advisory committee, and will consider all possible 


measures of ecosystem quality restoration for the Metro Toronto Waterfront. 


This report presents the results of the Main STP Pilot Study for 1986 and 1987, in relation 


to historical data (back to 1976), and in relation to concurrent data for the entire Toronto 


Waterfront (Etobicoke Creek to Rouge River). Sampling and chemical analysis were 
performed by the Ontario Ministry of the Environment (MOE). A portion of the data 
analysis and evaluation as described here was performed by Beak Consultants Ltd. under 
contract to MOE. 


1.1 Study Objectives 


In order to develop standard procedures for the derivation of receiving-water-based effluent 
limits, the MOE has evaluated numerous traditional and novel water quality assessment 
techniques. These techniques have been applied on a site-specific basis at the Toronto 
Main STP pilot site, as well as at the other Pilot Sites throughout Ontario, to determine 


their suitability in establishing water quality-based effluent limits. 


Two distinct goals were expressed at the outset of these Pilot Site investigations, which were 
initiated in 1986. These goals are: 


e To develop site-specific effluent requirements for hazardous and conventional 


pollutants from specific dischargers, based on receiving-water impact assessment. 


® To develop and document, for application to other comparable discharger and 
receiving-water situations in the province, an acceptable set of receiving-water 


assessment techniques and procedures. 


As a basis for determining the impacts of toxic contaminants on the aquatic environment, 
the Toronto Main STP study team has relied on Provincial Water Quality Objectives 
(PWQO’s) and Guidelines (PWQGs), and advisories as well as water- and biologically-based 


criteria, developed by other government agencies to determine acceptable conditions. 


Specifically to the Main STP pilot site, the objectives have also included determination of 
the relationship between in-lake contaminant concentrations and causative factors such as 


local effluent inputs, wet weather inputs, advection, dispersion, external sources and 
sedimentation; and the estimation of the impact of input reductions on the zones of 
non-compliance with PWQOs and PWQGs in the vicinity of the Main STP. 


1.2 Site Description and Inputs 


The Toronto Main STP is located on the east side of the Metro Toronto Waterfront, which 
stretches from Etobicoke Creek at the western limit of Metro Toronto to the Rouge River 
at the eastern limit of Metro Toronto (Figure 1.1). With a design flow of 8.18 x 10°m’/d 
(180 MGD), the Main STP is the largest single municipal discharge, not only to all of 
Metro Toronto, but to the entire Lake Ontario shoreline. 


The 45 km Metro Toronto Waterfront shoreline provides many attractions to the 2.1 million 
population of Metro Toronto as well as a large tourist population. Public attractions 
include Exhibition Park, Ontario Place, Centre Island and Harbourfront, as well as many 
swimming beaches and protected mooring for over 7,000 small craft (RAP, 1988). The 
watershed that drains to the lake includes the drainage basins of Etobicoke Creek, Mimico 
Creek, the Humber and Don Rivers, Highland Creek and the Rouge River. These are 
separated into western and eastern sectors by the Port of Toronto which encompasses the 
Don River, Ontario Place, the inner and outer harbours, the Toronto Islands and the 
Eastern Headland. Information on the geology, soils and land use of these watersheds is 
contained elsewhere (RAP, 1988); in brief, the lower reaches of all these watersheds are 
highly urbanized, while some agricultural area is still found in the upper (northernmost) 
areas. 


Much of the western and eastern portions of the Lake Ontario shoreline in Metropolitan 
Toronto is relatively inaccessible to the public due to extensive private ownership, as well 
as rapid drops to the water’s edge (especially at the Scarborough Bluffs). Two lakefill 
projects in the west (Colonel Sam Smith Park and Humber Bay Park), and one in the east 
(Bluffers Park) have provided increased public access to these parts of the waterfront. East 
of the Humber River are found recreational areas such as the Western Beaches as well as 


sailing and rowing clubs and marinas. 
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The Central Waterfront is dominated by Toronto Inner and Outer Harbours, which are 
enclosed by natural and man-made features. While the Toronto Islands, which enclose the 
Inner Harbour, are largely natural, the Outer Harbour is delineated by the man-made East 
Headland, which has been constructed since 1956 by using largely clean fill from 
construction projects. The East Headland forms the western boundary of the Main STP 
area (Figure 1.2) and deflects westerly lake currents southward in the vicinity of the Main 
STP outfall. Near the Main STP, Ashbridges Bay is a popular area for pleasure boating, 
and the Eastern Beaches, a popular area for body-contact recreation, stretches eastward to 
the edge of the Scarborough Bluffs. 


1.2.1 Description of the Toronto Main STP 


With a capacity of 8.18 x 10°m*/day, the Toronto Main STP is presently the largest 
complete sewage treatment plant in Canada, and has an estimated value of approximately 
$400,000,000. 


The initial primary treatment plant was constructed in 1910. Since 1943, the plant has been 
upgraded in stages, and its ultimate capacity could reach 1.36 x 10°m’/day before the end 
of the century. (Gore and Storrie, 1986b). 


The wastewater treatment processes at the Toronto Main STP are described briefly here: 


e Grit removal and screening - Grit and non-organic material are removed from the 
wastewater in 16 grit channels. Additional suspended debris are then removed by 


14 automatic self-cleaning bar screens. Debris is hauled to a sanitary landfill site. 


e Primary settling - Solids in the raw wastewater are settled in 12 primary 
sedimentation tanks with a total volume of approximately 1.13 x 10°m’*. The resulting 
sludge is pumped to sludge blending tanks when it is mixed with a pre-selected ratio 
of digested waste activated sludge, or to anaerobic digestion tanks. 
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e Secondary (activated sludge) treatment - The primary effluent is mixed with return 
activated sludge and aerated in nine aeration tanks (total volume 1.6 x 10°m’), where 
air from centrifugal blowers is diffused in to provide an aerobic environment for 
biological growth. The sewage is then settled in nine sedimentation tanks (total 
volume 1.05 x 10°m’). The settled sludge is partially returned to the aeration tanks 
as a biological seed, while the remainder is pumped back to the head of the plant 
or to dissolved air flotation tanks for thickening. 


e Chlorination - The secondary effluent is disinfected with a chlorine solution, and is 
then discharged to the lake. Contact time required for chlorination to take effect 
is provided in the long outfall pipe. Chlorination is utilized only between April 1 
and September 30 of each year. 


® Sludge treatment - Waste activated sludge is first thickened from about 1% solids 
to 5% solids in ten dissolved air flotation units. It is then digested for several days 
in anaerobic digestion tanks at a controlled temperature of about 35°C. The sludge 
is then conditioned at a temperature of 215°C and a pressure of 2500 kPa in order 
to break down the sludge cells and release bound water. This thermally conditioned 
sludge (approximately 12% solids) is then dewatered on rotary drum vacuum filters 
or stainless steel coil media vacuum filters on belt filter processes. The dewatered — 
Sludge cake is then incinerated at approximately 800°C in multiple-hearth 


incinerators. The resultant ash is removed by a private firm to reclaim the minerals. 


Waste gases from the activated sludge treatment or sludge incineration are scrubbed with 
water including ozone or sodium hypochlorite and exhausted to the atmosphere through a 
185 m high stack. 


12.2 Other Inputs 


The following description of inputs to Lake Ontario in Metro Toronto is taken from the 
draft Stage 1 Toronto Waterfront Remedial Action Plan (RAP, 1988) with minor 


modifications. 


Metro Toronto has an extensive sewer system consisting of combined sewer and storm 
sewer networks. Storm sewers serve to facilitate efficient overland drainage. Urban storm 
runoff is conveyed by storm sewers directly to the Metro Waterfront and its tributaries. 
The volume and rate of flow varies with the duration, intensity and areal extent of storms, 
and the time interval between each successive event. The locations of major storm and 


combined sewer outfalls along the central waterfront are shown in Figure 1.3. 


Older areas of development within Metropolitan Toronto were initially serviced with 
combined sewers which were designed to carry both sanitary sewage and stormwater runoff. 
During dry weather, sanitary sewage is conveyed to area sewage treatment plants (STPs) 
for treatment. During wet weather, stormwater runoff enters the combined sewers and 
mixes with the sanitary sewage and is termed combined sewage. During small storms, which 
do not exceed the capacity of the system, the combined sewage is to be conveyed to the 
STPs for treatment. Combined sewer overflows (CSO’s) occur when the capacity of the 
system is exceeded. Excess combined sewage is then discharged directly to the receiving 
waters of the Humber and Don rivers and Lake Ontario. Although the City of Toronto, 
the City of Scarborough and the Borough of East York have undertaken major storm sewer 
separation programs, a large portion of the stormwater in older areas still enters combined 
sewers. The Humber and Don Rivers (mouths at stations 2043 and 1375 respectively) are 
the two most significant point sources of pollution other than the STPs described below. 


Storm and combined sewers may also carry flows during dry weather as a result of minor 
tributaries, groundwater infiltration, illegal sanitary connections or malfunctioning CSO 
regulators. 
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Domestic sewage from within the study area is conveyed by sanitary sewers to four STPs: 


e Humber STP (outfall station 2072) 

e Toronto Main STP (outfall station 1419) 

e Highland Creek STP (outfall station 2902) 

e North Toronto STP 
The Humber, Main and Highland Creek STPs discharge treated effluent directly to Lake 
Ontario. The North Toronto STP discharges to the Don River. These four STPs are 
operated by the Regional Municipality of Metropolitan Toronto. The Lakeview STP (63 
MGD) lies outside the study area, west of Etobicoke Creek in the city of Mississauga. 


The four STPs serving the Toronto area are conventional activated sludge plants with 
continuous phosphorus removal. Additional information on these plants is provided in 
Table 1.1. 


There are no direct discharges of industrial process waters to the waterfront (only cooling 
water discharges are allowed), but there are periodic spills of oils and other material into 
the Spadina and Simcoe Street Slips, probably from the waterfront railway lands and 
adjacent urban areas. 
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2.0 STUDY OUTLINE AND METHODOLOGY 


The water quality component of the Main STP MISA pilot site study consisted of collection 
of samples from the Main STP effluent as well as the receiving water in Lake Ontario. 
These were done during a series of intensive sampling periods, each of one to three weeks 
duration in 1986 and 1987. (Table 2.1). | 


2.1 1986 Study 


The 1986 studies were considered as preliminary to the more detailed work done in 1987. 
As such, the sampling effort was concentrated in the immediate vicinity of the Main STP 
outfall (Station 1419, plus stations 1431, 1432, 1433 and 1434, which were located on a 
circle of 250 m radius centered on the outfall; Figure 1.2). Two "control" stations were 
located outside the Main STP study area: 1536 (Island WTP intake) and 2029 (R.C. Harris 
WTP intake) (figure 1.1). All stations near the STP outfall were sampled in duplicate two 
to three times daily; on most study days only one "control" station was sampled. In addition, 
an hourly series of six to eight duplicate final effluent grab samples was collected from the 


manhole prior to discharge to the lake, on the same days as lake water collections. 


At each visit to a lake station, a profile for temperature, conductivity and turbidity (as 
estimated by fluorescence) was first obtained. The duplicate water samples were collected 
at the surface (generally 0.5 or 1.0 m water depth) and at the bottom (6 m) or depth of 
turbidity maximum. The depth of turbidity maximum was thought to indicate possibly 
higher contaminant concentration, as bottom turbidity plumes have frequently been found 


in previous work (Boyd, 1988). 


At the in-plant final effluent sampling location (station 03-0022), station 1419 and 1431 to 
1434, a complete set of parameters was collected: general chemistry, bacteriology, heavy 
metals and toxic organics (PCBs, organochlorine pesticides, chlorobenzenes, chlorophenols, 
organo- phosphorus insecticides, phenoxy acid herbicides, and volatile organics). At stations 
1431 to 1434, only a limited number of general chemistry and bacteriology parameters 


were sampled. 


TABLE 2.1; SAMPLING PERIODS, TORONTO MAIN STP MISA SURVEY 


In-Plant 


(a) 1986 June 6 
August 12,13,14 
September 8,9,10,11 
November 4,6,7 


(b) 1987 June 2-5, 15-19 
July 28-31 
August 4, 11-12 


October 5-8, 13-16 


Note: W = Western portion of waterfront 
C = Central portion of waterfront 
E = Eastern portion of waterfront 


Lake Ontario 


June 6 

August 12,13,14 
September 8,9,10,11 
November 4,6,7 
June 4,7,14 (W) 
June 3,6,11 (C) 
June 2,5,10 (E) 
August 1,7,11 (W) 


July 30, August 2,6 (C) 
July 29, August 4,8 (E) 


October 7,11,13 (W-C) 
October 6,10,12 (C-E) 


2.2 1987 Study 


In 1987, the survey area encompassed the entire Metro Toronto waterfront area (Figures 
1.1 and 1.2), in support of both the Main STP MISA study and the Toronto Waterfront 
RAP. The purposes of this survey were: (a) to provide further baseline data on water 
quality across the entire Metro Toronto Waterfront area, (b) to update water quality 
zonation in the Toronto Waterfront; (c) to provide an estimate of zones in non-compliance | 
with PWQO, in order to estimate the extent of loading reductions required, and (d) to 
provide input for an organic contaminant modelling study being undertaken by Dr. E. 
Halfon of NWRI (see section 4.3). The survey also included in-plant sampling (daily 
duplicate grab samples) at each of the four Metro Toronto area STPs discharging to Lake 
Ontario (Section 1.2), as well as tributary mouth sampling at all Metro area tributaries 
(section 1.2). Loading information from these components was required for the 


contaminant model, as well as for the Toronto Waterfront RAP. 


Three survey intervals of three weeks’ duration were performed; these occurred in June, 
July-August and October, 1987 (Table 2.1). During each survey, each station was visited 
on three separate dates; however, two to three days were required to visit all stations once 
due to the large area and large number of stations to be sampled. This also resulted in 
profiling activity being limited to a few selected stations. At each station, duplicate samples 
were collected for general water chemistry, bacteriology and heavy metals. Based on 1986 
data (section 3.2.2), only stations closest to dischargers were sampled for organic scans, and 
these scans were limited to PCBs, organochlorine pesticides and chlorobenzenes. However, 
a larger range of stations were sampled for organics using ultra-trace methodology (see 
below). The in-plant effluent samples were also analyzed for PCBs, organochlorines, 
chlorobenzenes, chlorophenols, triazine herbicides and volatile organics on a series of dates 
during each of the three survey periods. 

Analysis for organic contaminants at ultra-trace levels was also conducted on lake and 
effluent samples during the three 1987 survey periods. This was required because nearly 


all the parameters sought in 1986 were rarely if ever, found in the lake water and even only 


occasionally in the STP effluent (Sections 3.1.1 and 3.1.2). Quantitative data were also 


required for model calibration and validation of the fate and transport model (section 4.3). 


At about 20 lake stations in each survey period, 200-L of lake water was pumped through 
a centrifuge into an APLE (Aqueous Phase Liquid-liquid Extractor; Oliver and Nichol, 
1986), in which the organic contaminants were extracted into 8 L of dichloromethane and 
the dichloromethane returned to the laboratory for analysis. On eight dates during the 1987 
survey periods, 16 L final effluent samples from the four STPs were also collected (Eli-Eco 
Laboratories, 1988). Daily composites were obtained by collecting 8 L in the morning and 
8 L in the afternoon. These samples were sequentially extracted with three aliquots of 
dichloromethane, which were then combined for analysis. Analysis of extracts from either 
the sewage or lake samples involved clean-up and concentration steps followed by analysis 
using a capillary gas chromatograph with an electron capture detector. Quantification was 
achieved by correcting the measured results for the recovery of brominated organic 
compounds with which the samples had been spiked immediately after collection. Further 
analytical details are given elsewhere (Eli-Eco Laboratories, 1988; Oliver and Nicol, 1986). 
Field work (on lake) for this portion of the study was done by Lakes Research Branch, 
NWRYI, collection of STP samples and all analytical work was done by B.G. Oliver (Eli-Eco 
Laboratories) under contract. Similar 200 L samples of tributary water were also collected 
in August and October 1987 and analyzed under contract by Zenon Laboratories, 
Burlington, Ontario. 


3.0 RESULTS AND DISCUSSION 


The results of the 1986-87 Toronto MISA water quality sampling program are discussed 
along with data obtained in related studies and historical data in this section where 
appropriate. The discussion is divided into subsections dealing with the Main STP itself, 
Lake Ontario in the Main STP vicinity, and Lake Ontario across the entire Toronto 
Waterfront area. The latter section places the relative impact of the Main STP discharge 
into context with the entire waterfront, and is also important in terms of the Toronto 
Waterfront RAP. 


3.1 Main STP Effluent 


This section is primarily concerned with loadings of conventional and toxic pollutants from 
the Main STP to Lake Ontario as determined during the Great Lakes Section surveys in 
1986 and 1987. Additional heavy metal loadings were also obtained from data collected by 
the MISA Municipal Section (MOE, 1989). 


3.1.1 Monthly and Annual Loadings to Lake Ontario 


Concentrations of all water quality parameters (conventionals, bacteria, metals and organics) 
measured during the 1986 and 1987 surveys were multiplied by the effluent flow data 
corresponding to the time of collection of each effluent sample (nearest 15 minutes in 1986, 
nearest hour in 1987). These produced time series of 7 to 8 hourly loadings per day in 1986 
and duplicate daily results in 1987. Within-day variations in loadings as obtained from 1986 
data are discussed in section 3.1.3; in this section, monthly and annual loadings are 
discussed. These are compared with loadings from the other Toronto-area STP’s 
discharging to Lake Ontario during 1987 in section 3.3.3. 


Annual statistics for loadings of conventionals, bacteria, heavy metals and organic 
compounds (conventional and ultra-trace dectection limits) are given in Tables 3.1 for 1986 
and 3.2'for 1987. Data for organic compounds at conventional detection limits are limited 


to those found in 50% or more of the effluent samples. 


3-1 


For calculation of means and standard deviations in these samples, the maximum likelihood 
estimator (MLE) method (El-Sharaawi, 1989) was used to estimate the mean loadings and 
their confidence limits. Other methods of handling such samples as discussed in Appendix 
4, 


Organic compounds which were infrequently or never detected in the Main STP effluent 
at conventional detection limits are summarized in Table 3.3. In general, seasonal and 
short-term variability in effluent loadings (as a result of variations in both flow and 
concentration) were found to be high; consequently no statistically significant (P>0.05) 
differences were found between years in loadings results. 


Loadings for conventionals and heavy metals have been summarized on a monthly basis for 
the survey months in 1986 and 1987 in Tables 3.4 - 3.10 and for selected parameters in 
Figure 3.1. Monthly bacterial loadings are summarized in Table 3.1b and 3.2b. Nitrate 
+ nitrite-N was found to be highest during July to September of both years, suggesting 
some nitrification of effluent during the warmest months; however, the extent of nitrification 
was very small at any time compared to Lakeview and Highland Creek STPs (section 3.3.3). 
Total phosphorus, BOD, COD and suspended solids show some tendency towards higher 
loadings early in the season; however, variability is great. Heavy metal loadings were 
slightly higher in 1986 than in 1987, but variability in 1986 was greater. Bacteria loadings 
during fall post-chlorination periods (November 1986 and October 1987) are one to three 
orders of magnitude higher than for the other survey months, as expected for chlorine 
disinfection; the difference was greatest for fecal coliforms and E.coli. Unexpectedly, high 
bacterial loadings were found on some dates in September 1986 (Table 3.1b); the reason 


for this is unknown. 


TABLE 3.1(a): LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
MAIN STP, ANNUAL MEANS, JUNE - NOVEMBER 1986 


STD. 

PARAMETER MEAN DEV MIN MAX N 
Flow (10°m?/day) 1260 210 790 2080 ts 
(a) Conventionals (10° kg/day) 

NH3-N 26.0 7.4 7.0 42.6 148 
TKN 30.8 8.8 10.0 63.8 148 
(NO2+NO3) -N 1 Y 1.70 0.06 7.3 148 
Total P eg 220 O23 125 148 
Total Dissolved Solids* 776. 135. 488. 1221. 90 
Suspended Solids 30.4 53.0 Bal 338 140 
BOD 15.9 Ly 0.4 81 126 
COD 51.9 2145 20.3 T352. 119 
DOC 14.3 4.3 0.12 Jap EEN 142 
Phenols (kg/day) 5.8 5 a7 ae 9 34. 129 
Cyanide (kg/day) 14.2 17. 4l. 108. 126 
(b) Heavy Metals (kg/day) 
Cadmium ‘Teg! 31.8 0. 11 146 
Chromium BWA. 49, 1. 282 146 
Copper 4l. ey ae 1G 290 146 
Iron 1300. 2000. SH PS 12400. 146 
Mercury (g/day) 107. 136. 11. 700 146 
Manganese 68. 18. vhs 135 146 
Nickel $30 13). ae 68. 146 
Lead oe Fah <3. 118. 135 
Zinc 88. 116. 1. 614. 146 


* Estimated as conductivity (25°C) x 0.65 


TABLE 3.1(b): BACTERIAL LOADINGS (organisms/day) 


MAIN STP, 1986 


(a) Whole year 
Aug. -Nov. 
Fecal coliforms 
ES coti 
Fecal streptococci 
P. aeruginosa 


(b) August 
Fecal coliforms 
E. colj 
Fecal streptococci 


Pp, aeruginosa 


(c) September 
Fecal Coliforms 
E. coli 


Fecal streptococci 
Pp, aeruginosa 


(d) November* 
Fecal coliforms 
E, coli 
Fecal streptococci 


Pp. aeruginosa 
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Note * - effluent not chlorinated 


TABLE 3.1(c): LOADINGS OF ORGANIC COMPOUNDS (CONVENTIONAL DETECTION 
LEVELS), MAIN STP, ANNUAL MEANS, 
JUNE - NOVEMBER 1986 (g/day) 


MLE ESTIMATES OBSERVED 
CONF. INT. 
PARAMETER MEAN MIN MAX MEDIAN MAX N % 
Dichloromethane 14700. 7400. 29000. 8400. 58000. 28 79 
Chloroform 4700. 2800. 7900. 2600. 29000. 28 71 
y-BHC (Lindane) oo. 19. 66. 222 88. 34 79 
1,2,4-Trichlorobenzene 2268 70. 740. IDG 290. 33 73 
1,2,3-Trichlorobenzene 105. 43. 250. SLs 330. 34 74 
1,2,3,4-Tetrachlorobenzene 58. SHEN 100. 34. 91. 34 85 
Pentachlorobenzene 20. 10. 44. 11. 60. 34 70 
2,4,6-Trichlorophenol 165. 115. 240. 140. 700. 36 TZ 
2,3,4,5-Tetrachlorophenol 102. 4 150. BIS h 1290. 36 50 
Pentachlorophenol 490. 330. 730. 360. 1900. 26 88 
Diazinon 230. 180. 290. 230. 370. 17 100 


Note: % = percent of effluent samples in which the indicated constituent 
was detected (including samples reported as "trace" - <T) 


TABLE 3.2(a): LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
MAIN STP, ANNUAL MEANS, JUNE - OCTOBER 1987 


STD. 
PARAMETER MEAN DEV MIN MAX N 
Flow (10°m?/day) 960 250 560 1730 48 
(a) Conventionals (10° kg/day) 
NH3-N 20.0 6.5 8.9 34.6 47 
TKN 24.1 6.9 1233 37 aL 45 
(NO2+NO3) -N 0.86 2.19 0.03 110. 48 
Total P 0.79 1.69 0.10 10.89 47 
Total Dissolved Solids* eis: 142. 355. 961. 47 
Total Filtered Residue Avion 122 2152 734. 46 
Suspended Solids path Ak 75.6 SIRS 398. 47 
coD 87.2 81.0 33.9 442. 41 
DOC 13.4 3.4 6.1 Pda ak 47 
(b) Heavy metals (kg/day) 
Chromium 223 We, 6. 74. 46 
Copper L6- L52 0. =). 46 
Iron 980. 670. 270. 3700. 46 
Mercury (g/day) 39. Zon 10. Lie 46 
Manganese 70. 22s 36. 124. 46 
Nickel 34. TES he 7. 65. 46 
Lead 9. 10. 0. 34. 46 
Zine 34. oe 0. 94. 46 


* Estimated as conductivity (25°C) x 0.65 


TABLE 3.2(b): BACTERIAL LOADINGS (organisms/day) , 
MAIN STP, 1987 


GEOM. SD SD 
PARAMETER N MEAN LOWER UPPER MIN MAX 
(a) Whole year 
Fecal coliforms 46 2.8x1015 8 .7x1012 9.5x1014 2.0x101? 1.7x101§ 
E, coli 46 2.0x10}% 5eoxl0** JP 7K102 1.7x10}2 1.5x10? 
Fecal streptococci 46 1.1x1075 7.6x1013 1.8x10** 4.5x10?° 9.8x101* 
Heterotrophs 10 6.3x1014 Qmix Os 1.8x103> 1.1x10** 1.9x1015 
P, aeruginosa 46 5.0x1012 4..3x10?° 5.8x10!2 2.4x101° 6 .0x103 
(b) June 
Fecal coliforms 16 1.1x101 1.4x10}4 9.7x1015 1.9x10}2 2.2x1015 
E, coli 16 9.5x10}4 ee xLo* 7.6x10}3 1.6x10}4 1.6x10}> 
Fecal streptococci 16 1.9x1035 = 55..8x 1032 6.5x1019 2° 7210-4) 6 3 ox 107 
Heterotrophs 10 6.3x104 2.1x1014 1.8x101° 1.1x1014 1.9x10}5 
P, aeruginosa 16 4.7x1011 5.3x10?° 4..1x1012 4..1x10?° 5.7x10}5 
(c) July-August 
Fecal Coliforms 14 1.0x102 3.9x101! 2.6x1012 2.0x1012 4.9x10}2 
E, coli 14 5.2x1012 2.1x101 1.3x10?4 1.7x1012 2.7x1012 
Fecal streptococci 14 4e6xL0+*) eligx107* 1.7x1012 4.5x107° 2.7x10}2 
P, aeruginosa 2) aR 5.3x101° 2.7x10?° 1.1x1012 2.4x10?° 1.7x101! 
(d) October* 
Fecal coliforms 16 1.3x1015 1.7x10+4 1.1x1016 6.4x1018 1.7x101§ 
E, coli 16 elx10- 1.2x10%4 9.8x10?° 4.6x10}% 1.5x10?6 
Fecal streptococci 16 121 Ol 1. 9x 1.02° 7.9x1014 7.8x1012-9..8x101* 
P, aeruginosa 16 3.8x10}2 5.4x101! 2.7x1035 4.1x1011 6.0x10?9 


Note * - effluent not chlorinated 


TABLE 3.2(c): LOADINGS OF ORGANIC COMPOUNDS (CONVENTIONAL DETECTION 
LEVELS), MAIN STP, ANNUAL MEANS, JUNE - OCTOBER 1987 (g/day) 


MLE ESTIMATES OBSERVED 
CONF. INT. 

PARAMETER MEAN MIN MAX MEDIAN MAX N % 
a-BHC 4.1 E35 Lin 1.6 26. 20 60 
7y-BHC (Lindane) 4l. 26. 65. Sho 78. 20 95 
1,2,3-trichlorobenzene 37a 2.6 530 ie 81. 20 55 
1,2,4-trichlorobenzene 28. 1.0 800 Loe y fe ye 20 55 
1,2,3,4-tetrachlorobenzene 315 ay, 126 Le 63. 19 74 
pentachlorobenzene Sag S22 11 3.9 18. 20 75 
hexachlorobenzene Sao 1.8 16. LieZ 51s 19 63 
2,4,6-trichlorophenol Hates 382 132% 42. 200. 18 83 
pentachlorophenol 160. 80. 320. 62. 820. 18 95 
chloroform 2700. 2100. 3600. 2400. 9300. 44 86 
dichloromethane 5300. 4100. 6800. 4500. 13000. 44 89 
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Note: % = percent of effluent samples in which the indicated consituent 
was detected. 


TABLE 3.2(d): LOADINGS OF ULTRA-TRACE ORGANIC COMPOUNDS 
MAIN STP, ANNUAL MEANS, 1987 (g/day) 


STD. DETECT. % OF 
PARAMETER MEAN DEV. MIN MAX N LIMIT SAMPLES 


(ng/L) DETECTED 


1,3-dichlorobenzene 78. 38. 25. 155. 10 1 100 
1,4-dichlorobenzene 630. 260. 100. 990. 10 1 100 
1,2-dichlorobenzene 280. 130. 63. 480. 10 1 100 
1,3,5-trichlorobenzene 2.5 2.1 0.4 7.4 10 0.1 100 
1,2,4-trichlorobenzene 116. 118. 15 350. 10 0.1 100 
1,2,3-trichlorobenzene 54. 79. 6.1 210. 10 0.1 100 
1,2,3,5-tetrachlorobenzene 9.2 15.8 0.1 42. 10 0.05 100 
1,2,4,5-tetrachlorobenzene 19. 15. 4.3 47 10 0.05 100 
1,2,3,4-tetrachlorobenzene 68. 101. 12. 271. 10 0.05 100 
pentachl orobenzene . 6.9 5.1 0. 17 10 0.02 80 
hexach lLorobenzene 4.1 6.1 0. 17. 10 0.02 80 
hexachloroethane 0.16 0.33 0. 1.1 10 0.02 50 
2,4,5-trichlorotoluene 0.80 ene 0. 0.70 10 0.1 50 
2,3,6-trichlorotoluene 0.25 0.42 0. 1.2 10 0.1 40 
pentachlorotoluene 0.67 0.84 0. eat 10 0.02 70 
a-BHC 1.3 0.8 0. 2.5 10 0.02 90 
4-BHC (lindane) 21. 17. 0. 50 10 0.02 70 
heptachlor + heptachlor epoxide 0.11 0.17 0. 0.49 10 0.02 - 40 
endosulfan - I 0.42 0.43 0. 1.3 10 0.02 60 
endosulfan - II 0.10 0.22 0. 0.63 10 0.02 20 
4-chlordane 0.34 0.18 0. 0.67 10 0.02 90 
a-chlordane 0.43 0.42 0. 1.3 10 0.02 70 
nonachlor 0.27 0.24 0. 0.79 10 0.02 70 
pp-DDE 0.62 0.42 0. 1.5 10 0.02 90 
pp-DDD 0.37 0.64 0. 2.1 10 0.02 60 
op-DDT 0.13 0.32 0. 1.0 10 0.02 20 
Total PCB 16.8 10.3 4.7 40. 10 0.05 100 


TABLE 3.3: ORGANIC PARAMETERS SELDOM OR NEVER DETECTED 
AT CONVENTIONAL DETECTION LEVELS, MAIN STP, 1986 


DETECTION FREQ. OF MAXIMUM 
PARAMETER LIMIT DETECTION CONCENT. 
(ng/L) (%) (ng/L) N 
(a) PCBs, organochlorine pesticides 
Total PCB 20 0 - 29 
Hexachlorobenzene 1 43 6 29 
Heptachlor 1 0 - 29 
Aldrin 1 0 - 29 
pp-DDE 1 20 4 29 
op-DDT 5 0 - 29 
pp-DDD 5 0 - 29 
pp-DDT 5 0 : 29 
Mirex 5 0 - 29 
a- BHC 1 13 2 29 
&-BHC 1 7 5 29 
a-chlordane 2 3 22 29 
y-chlordane 2 7 20 29 
Oxychlordane iz 0 - 7h) 
DMDT -methoxychlor 5 0 - 29 
Heptachlor epoxide 1 0 - 29 
Endosulfan I Z 0 - 29 
Endosulfan II 4 13 14 29 
Endosulfan Sulfate 4 0 - 29 
Dieldrin 2 3 6 29 
Endrin 4 7 8 29 
(b) Chlorobenzenes 
Hexachloroethane 1 27 17 30 
1,3,5-Trichlorobenzene 5 17 iM 30 
2,4,5-Trichlorotoluene 5 0 ~ 30 
2,3,6-Trichlorotoluene 5 0 - 30 
Hexachlorobutadiene 1 0 - 30 
1,2,3,5-Tetrachlorobenzene 1 23 30 30 
1,2,4,5-Tetrachlorobenzene 1 23 40 30 
2,6,a-Trichlorotoluene 5 0 - 30 
(c) Chlorophenols 
2,4,5-Trichlorophenol 50 5 50 30 
2,3,4-Trichlorophenol 100 3 120 30 
2,3,5,6-Tetrachlorophenol 50 Zi 1600 30 


TABLE 3.3: (cont’d) 


DETECTION 
LIMIT 


(ng/L) 


FREQ. OF 
DETECTION 
(%) 


(d) 


(e) 


(f) 


Phenoxyacid Herbicides 

Dicamba 

2,4-D-propionic acid 
2,4-Dichlorophenoxyacetic acid 
Silvex 
2,4,5-Trichlorophenoxyacetic acid 
2,4-Dichlorophenoxybutyric acid 
Picloram 


Organophosphorus Insecticides 
Dichlorovos 
Dursban 

Ethion 

Guthion 
Malathion 
Mevinphos 
Methylparathion 
Methyltrithion 
Parathion 
Phorate 

Reldan 

Ronnel 


Carbamate Herbicides 
Carbofuran 

CiPC 

Diallate 

Eptam 

IPC 

Propoxur 

Sevia 

Sutan 
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TABLE 3.3: (cont’d) 


DETECTION FREQ. OF MAXIMUM 
PARAMETER LIMIT DETECTION CONCENT. 
(ng/L) (%) (ng/L) N 
(g) Volatile Organics 

1,1-Dichloroethylene 1 0 - 24 
TRS -1,2-Dichloroethylene 1 0 - 24 
1,1-Dichloroethane 1 0 - 24 
1,1,1-Trichloroethane 1 24 Fas y 24 
1,2-Dichloroethane 1 0 - 24 
Carbon Tetrachloride 1 0 - 24 
Benzene 1 5 Leer 24 
Trichloroethylene D 0 - 24 
Dichlorobromomethane 1 0 - 24 
Toluene 1 10 17 24 
1,1,2-Trichloroethane 1 0 - 24 
Chlorodibromomethane 1 0 - 24 
Chlorobenzene 1 0 - 24 
Ethylbenzene 1 0 - 24 
m- and p- Xylenes 1 14 75 24 
0-Xylene 1 14 17 24 
Bromoform 1 0 - 24 
1,1,2,2,-Tetrachloroethane 1 14 Sul 24 
1,4-Dichlorobenzene 1 14 2a 24 
1,3-Dichlorobenzene 1 sis} Sat 24 
1,2-Dichlorobenzene 1 5 1T 24 


TABLE 3.4: LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
MAIN STP, JUNE 1986 
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STD. 
PARAMETER MEAN DEV MIN MAX N 
Flow (103m?/day) 1190. 50. 1130. 1250. 7 
(a) Conventionals (10° kg/day) 
NH3-N 28.8 Ze 25.1 32.6 10 
TKN 33.6 3.2 28.1 40.2 10 
(NO2+NO3) -N 0.06 0.01 0.05 0.07 10 
Total P 0.92 0.08 0.75 1.01 10 
Total Dissolved Solids* A922. 36. 744. 854. 10 
Suspended Solids 27.4 3.4 23.4 35.1 10 
DOC 21.4 a hiy 16.2 25.2 10 
(b) Heavy Metals (kg/day) 
Cadmium 1.3 0.4 1.0 2.0 10 
Chromium 38. th 35. 46. 10 
Copper 52/, 9. 35x 66. 10 
Iron 1900. 300. 1600. 2700. 10 
Mercury (g/day) 117. 44, 90. 234. 10 
Manganese 61. 4. 57. 70. 10 
Nickel 39. 6. 34. 50. 10 
Zinc 64. 6. 56. 78. 10 


* Estimated as conductivity (25°C) x 0.65 


TABLE 3.5: LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
: MAIN STP, AUGUST 1986 


STD. 
PARAMETER MEAN DEV MIN MAX N 
Flow (10°m?/day) 1200. 180. 790. 1600. 20 
(a) Conventionals (10° kg/day) 
NH3-N 24.1 4.0 16.9 30.3 40 
TKN 2759 4.5 18.3 34.6 40 
(NO2+NO;) -N 2.06 1.87 0.49 (in Tf 40 
Total P 2.36 257 0.45 12.5 40 
Total Dissolved Solids* 683. 101. 488. 883. 24 
Suspended Solids 54.2 53.6 10.8 195. 40 
BOD 26.6 2a 9.0 49.3 32 
COD 73.6 30.3 47.5 192s: 21 
DOC L2ie Si] 0.1 PPE SE 34 
Phenols (kg/day) 70 4.1 2.8 Lae 40 
Cyanide (kg/day) 7a 26. i 108. - 30 
(b) Heavy Metals (kg/day) 
Cadmium 2.0 7s Pd 0. iD 40 
Chromium pon 46. Le 194. 40 
Copper 80. 54. ie 228. 40 
Iron 1700. 1600. OMA 7200. 40 
Mercury (g/day) 191. 153. 63. 63.7 38 
Manganese 67. 20. L 102. 40 
Nickel 4l. Ze 2: 63. 40 
Lead 20. 19 oir 74, 40 
Zinc 118. 109 ke, 423. 40 


* Estimated as conductivity (25°C) x 0.65 


TABLE 3.6: 


LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
MAIN STP, SEPTEMBER 1986 


Flow (10°m?/day) 


(a) Conventionals (10° kg/day) 


(b) 


NH3-N 
TKN 


(NO2+NO3) -N 


Total P 


Total Dissolved Solids* 
Suspended Solids 


BOD 
COD 
DOCG 


Phenols (kg/day) 
Cyanide (kg/day) 


Heavy Metals (kg/day) 


Cadmium 


Chromium 


Copper 
Iron 


Mercury (g/day) 
Manganese 


Nicke 


Om 
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* Estimated as conductivity 


(25°C) x Of65 


TABLE 3.7: LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
MAIN STP, NOVEMBER 1986 


STD. 
PARAMETER MEAN DEV MIN MAX N 

Flow (10%m?/day) 1250. 2102 890. 1800. 20 

(a) Conventionals (10° kg/day) 
NH3-N S00 B.S 227A 39.6 40 
TKN 37.0 4.1 25ne 47.2 40 
(NO2+NO3) -N 0.67 0.87 0.18 5.5 40 
Total P 0.32 0.05 0.23 0.43 40 
Total Dissolved Solids* 865. Pole 579% 1220. 40 
Suspended Solids 6.4 iso Bho 9.3 40 
BOD Gn pL! 13 9.9 40 
COD 46.7 8.5 28.6 64.9 40 
DOC 16.0 2.2 sb aak Zone 40 
Phenols (kg/day) one 0.6 19 4.6 37 
Cyanide (kg/day) Lon 13%.3 4. 65. 40 

(b) Heavy Metals (kg/day) 
Cadmium 0.6 0.1 0.0 1.0 39 
Chromium 18. LOZ iy 60. 39 
Copper ok, Ze ih 18. 39 
Iron 600. 100. 500. 900. 39 
Mercury (g/day) 523 16. 233 86. 40 
Manganese sO: 14. 46. 108. By] 
Nickel 322 8. 18. 47. 39 
Lead <4 0.6 < 3 35 39 
Zinc 83 Toss 20 469 39 


* Estimated as conductivity (25°C) x 0.65 


TABLE 3.8: LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
MAIN STP, JUNE 1987 


STD. 
PARAMETER MEAN DEV MIN MAX N 
Flow (10%m?/day) 910. 170. 560. 1120. 18 
(a) Conventionals (10° kg/day) 
NH3-N 18.5 3.9 10.2 234.5 17 
TKN 234) °6.3 12.3 39.0 18 
(NO2+NO3) -N 0.15 0.19 0.03 0.79 18 
Total P 1.19 2.63 0.17 10.89 18 
Total Dissolved Solids* D9: 109. 356. 750. 17 
Total Filtered Residue 507. 119. 300. 704. 17 
Suspended Solids 50.2 119. 5a2 398. 18 
COD Li2, L2Le 373 442. 17 
DOC 14.0 ihe 8.7 19.7 D7 
(b) Heavy Metals (kg/day) 

Chromium 16. 26 9 22e 16 
Copper : 16. oh 0 ite 16 
Iron 1050. 280. 600 1640. 16 
Mercury (g/day) aly Le 21 45. 16 
Manganese 89. 18. 62 124. 16 
Nickel 28. 9. 17 45. 16 
Lead ib De 30 21% 16 
Zinc 16 9. 0 34. 16 
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* Estimated as conductivity (25°C) x 0.65 


TABLE 3.9: LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
MAIN STP, JULY-AUGUST 1987 


STD. 
PARAMETER MEAN DEV MIN MAX N 
Flow (10%m?/day) 1090. 330. 610. 1730; 14 
(a) Conventionals (10° kg/day) 
NH-N 1992 7139 8.9 34.6 14 
TKN 2243 6.7 1235 Sie 13 
(NO2+NO3) -N 2.49 3.63 0.67 11.1 14 
Total P 0.22 0.12 0.10 0.42 14 
Total Dissolved Solids* 615. 180. 350% 961. 14 
Total Filtered Residue 486. 138. 2 oe 734. 14 
Suspended Solids 6.6 S36 oheve) 14.2 14 
COD 68. 222 34. 98. 9 
DOC lla ei 6.1 18.8 14 
(c) Heavy Metals (kg/day) 

Chromium LSe oF 6. 353 14 
Copper 28. 132 10. 52k 14 
Iron 640. 370. 270. 1300. 14 
Mercury (g/day) 24. 132 LOF 32% 14 
Manganese 68. Loe 36. 104. 14 
Nickel 312 132 18. 65. 14 
Lead 14 14. 0. 34 14 
Zinc 19 Toe 0 52 14 


* Estimated as conductivity (25°C) x 0.65 


TABLE 3.10: LOADINGS OF CONVENTIONALS AND HEAVY METALS, 
MAIN STP, OCTOBER 1987 
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STD. 
PARAMETER MEAN DEV MIN MAX N 
Flow (103m?/day) 900. 210. 620. 1210. 16 
(a) Conventionals (10° kg/day) 
NH3-N 2254 7.0 12.8 Pe) 16 
TKN 26.8 7.6 15.7 38.0 13 
(NO2+NO3) -N 0.21 0.14 0.03 0.51 16 
Total P 0.85 0.65 0.14 283 2 15 
Total Dissolved Solids* 530. 132% 374. 748. 16 
Total Filtered Residue 428. 102. 297. 594. 15 
Suspended Solids 18.5 24 ak) 70. 15 
coD 70.6 23.0 38.0 120. 15 
DOC 14.2 2 hy 10.6 Pah 72 16 
(b) Heavy Metals (kg/day) 
Chromium sya 18. ily es 74. 16 
Copper 18. 16. Me 55. 16 
Iron 1200. 1000. 500. 3700. 16 
Mercury (g/day) 61. 24. 35. LT. 16 
Manganese Do 10. 41. i: 16 
Nickel 42. 13. 26. ee 16 
Lead 4 Way 6. 7 25 16 
Zinc 66. 51. 14 194 16 


* Estimated as conductivity (25°C) x 0.65 
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Additional loadings data for heavy metals were available from a 37-STP pilot monitoring 
study (MOE, 1989). The results of that study are compared to the present results in Table 
3.11. For most metals, the present study showed higher loadings; however, nickel loadings 
were similar and chromium loadings were lower. The latter study also included removal 
efficiencies; it is seen that these are 90% or higher except for chromium (69%) and nickel 
(63%). Metal loadings from the Main STP were the highest of all STPs studied, likely due 
to the large throughput of sewage. However, for most metals, Main STP shows a higher 
than average removal efficiency compared to the other eight STPs for which detailed 
information was available; the lower removal efficiencies for chromium and nickel were 
about average. Mass balance studies in the same report showed that most of the metal 
removed was present in the sludge; they also concluded that variability in removal 
efficiencies between metals was likely due to differences in metal chemistry and 


biochemistry. 


3.1.2 Probability Statistics 


As is evident from the data presented in section 3.1.1, loadings of chemicals (conventionals 
or trace contaminants) from the Main STP are highly variable in nature. This section 
presents additional data on loading variability. This is important for two purposes: (1) 
Investigators working on toxicity studies, either in-situ or using laboratory serial dilutions 
of final effluent, need to know under what conditions their experiments were conducted - 
ie. were they using an average concentration or perhaps one whose concentration was 
higher than that observed 95% of the time; and (2) load allocation procedures to be 
developed based upon effluent limits derived from various studies must relate to this 
loading variability. For the purposes of load allocation, the allowable load is considered 
to be the load for which the defined objective will be satisfied for a specified fraction of the 
time at the edge of the mixing zone, given the statistical variability of effluent quality. The 
water quality-based effluent limits were derived on the basis of 95% compliance with 
objectives - i.e. effluent quality must not be out of compliance more than 5% of the time 


on a monthly average basis. From the various models used in their studies, they produced 


TABLE 3.11: COMPARATIVE HEAVY METAL LOADINGS (kg/day), 
TORONTO MAIN STP 


THIS STUDY MOE, 1989 
METAL 1986 1987 LOADING % REMOVAL 
Cr 37 am 58 68.8 
Cu 41 26 8 92982 
Hg (g/day) 107 64 31 O22 
Ni 33 34 27 63.2 
Pb - 12 8 89.5 
Zn 88 56 23 933 
e 


cumulative probability curves for loadings of various parameters, from which the fraction 


of non-compliance for any given load reduction scenario could be obtained. 


For conventional parameters, ammonia-N, total P and suspended solids were chosen, and 
probability plots for these parameters are presented for 1986 in Figure 3.2 and for 1987 in 
Figure 3.3. Total P and suspended solids are plotted as logarithmic transformations; as the 
distribution of ammonia is approximately normal, no transformation was required. The 
y-axis of each figure represents the percent of observation above the specified point on the 
line - ie. if one desires the 95th percentile ammonia loading (the loading exceeded by 5% 
of observations), this value is 3.8 x 10* kg/d in 1986 (Figure 3.2) and 3.0 x 10° kg/d in 1987 
(Figure 3.3). It should be noted that the minima and maxima of all parameters reflect the 
values given in Tables 3.1(a) and 3.2(a); in addition, the median loading (not given in the 
tables) is equal to the 50th percentile value, 2.6 x 10° and 2.0 x 10* kg/d for ammonia in 
1986 and 1987, respectively. 


Similar plots for several important heavy metals in 1987 are reproduced in Figures 3.4. All 
data were logarithmically transformed; usage of these plots is similar to that for 
conventionals. It should be noted that probability plots for two toxicologically important 
metals - cadmium and lead - were not prepared as most of the data were below detection 
for these metals. As for conventionals, a large range of loadings, due to the range in both 
concentration and flow, is evident from the figures. However, use of these figures will allow 
a clear perspective of the relative conditions occurring on any particular day of interest 


(e.g. date of toxicity experiment in the receiving water.) 


3.1.3 Daily Variation in Concentration and Loading 


The ammonia concentration in Main STP effluent tends to increase with time of day 
between 10:00 a.m. and 5:00 p.m. (Figure 3.5). However, there is a great deal of scatter 
in this relationship, due partly to variation in the timing of the daily cycle from day-to- 
day, and due also to amplitude variation between days. Amplitude variation between days 
can be large relative to within-day variation. For example, ammonia varied from 4 to 9 
mg/L over the course of the day on 11 September 1986, and from 25 to 30 mg/L on 07 


3 - 23 


cumulative percent 


‘cusulative percent 


cumulative percent 


Norma) Probability Plot 
Ammonie-N loading, 1986 Main STP 


i1 ioe’ gpGisay Vet 
| Ammonie-N loading, ko/day 


| Norma) Probability Plot 
' Log (total P loading), 1986 Main STP 


2.6 12.9 3.2 35 3.8 
' log (total P loading) kg/day 


Norma] Probability Plot 


' Log (susp.solids loading), 1966 Main STP 


3.8 4.2 4.6 5 5.4 
Log (suspended solids loading) kg/day 


(1986 | 


Probability distributions of loadings, 


FIGURE 3.2 


Toronto Main STP, 


‘cumulative percent 


cuaulative percent : 


Normal Probability Piot 
Log (Cu loading), 1987 Moin STP 


| 09 13 U7 | 2.1 2.5 
Leg (Cu loading), kg /day 


| Normal Probability Plot 
| Log (Hg loading), 1987 Main STP 


cugsulative percent 


1.4 1g 2.4 29 EWS RICE oo hie 
Log (Hg leading), kg/day 


Normal Probebility Plot 
Log (Zn loading), 1987 Main STP 


Log (Zn loading), kg/day 


Probability distribution of metal loadings 


1987 


Toronto Main STP, 


FIGURE 3.4 


cumulative percent 


cumulative percent 


cumulative percent 


99.9 


99.9 


95; 


Normal] Probability Plot 
Ammonia-N loading, 1987 Main STP 


Ammonia-N loading, kg/day 


Normal Probability Plot 
Log (Total P loading), 1987 Main STP 


2.2 2.5 28 3.1 
log (total P loading) kg /day 


Normal Probability Plot . 
Log (susp.solids loading), .1987 Main STP 


3.8 4.1 4.4 | 4.7 
log (susp. solids loading) kg/day 


1987 . 


Probability distributions of loadings, 


Toronto Main STP, 


FIGURE 3.3 


JOqWaAON £0 °6 Jaquisidas oT °9 isn3Bny hI °¢ 
JOqWaACN 90° Joquiaidas 690 «°S isn8Bny.¢] ‘Z 
Jaquiaides TT °Z Joqweizdas go . °h isn38ny ZI ‘J 

(ajdwes Jo ayep azedtpur spoqwiAs }o]d) 


9861 dLS NIVW NI DNITdWYVS SO AWILL HLIM DNIGVOT GNY NOLLVUYLNADNOO VINOWWY TVLOL LNaN 1444 >$°¢€ AYNSIA 


Say outs * say auty 
02 rel gl ra ol Sk S ae 0 0z Cer | 5°21 vl yd g es i 
Geo aces) Siti sez 1T eee Sie g2°9 GL'E z'T o2°St 62°91 eee] Se Stuy G38 62'9 or 52°) 

Sa Naa aa aaa Salat ated Dalal Related Deleted Dette Dette Det Dee Senet Sees eee eee ees Oe ee et Peo oot ot os os et ee ee Pee See Pee 
+ +G°E H t t 
' (i H H fb ' 
t ' + t +N 
H E H H “ l f 
J ‘ ‘ : 
+ js + ; cigs 
H H jt fh 
: : : +0058 
con; ‘ : 
+ #501 { 
' ; : e 
' ! ¢ z § 
: } i. ae 
' i 6 H ! 
¢ +o] ws : Z p 
3 i : z $ $ §S / 
9 $ € ! u + y tas 3S Res +ngt7 & 
i 13 H b) H F 5 a ~~ 
€ ! 1 ‘ C Sms <9 € N 
: € $ ey A : ces 9 t : 
i 0.7% e $$ é ' 
r oY ' J + y $ 1 +00L7 4 o 
i e 1 $ g : 3 ; $ $ 7 I 
7 1 $3 $ Es Ss : u : 1 6 g p 
+ f po § $. 8 8 +1Z > ' $ $ 6S 2 8 e 
yore, (Ce, 8 ' B) ; I $ 8 9 
z cts : u + ) $ +pnste 7 
ek: ! f) : 6 : 
$ $ 9 ' 7 : 8 
‘ E 3 6 #6°9Z a $2 
: 9 6 ‘ : 6 i 
; 5 8 & ' + $ +E 
6 : i 3 
: : 9 : 
; 8 +97 : 9 

g 6 : : 
6 54 ' a +HOSNF 
6 6 6 : i 9 : 
6 : i 9 9 
' ‘SIE ; 
aa Naat elt Del Related tale tote ee Dee Peete Pere Pee Pe! Pees Pood eee 2) Py ees een eee ee Ses ee eee ee ee eee See Dee Bee Pee Pe eee 


WANT) 315 UTey Z7 uOTjeI.~ $USN1S92 dIS URW 27 wotsezS 


see ‘e 


November 1986. The pattern of variation in TKN closely resembles that for ammonia. 
The difference in ammonia and organic nitrogen concentrations in the effluent are 
seasonally related to temperature because the conversion of ammonia to nitrate nitrogen 
is bacteriological. 


Phosphorus shows a different pattern of variation (Figure 3.6). While the extreme 
concentrations, representing upset conditions, do tend to be more extreme later in the day, 
most of the samples cluster at the low end of the concentration scale, even late in the day. 
Sometimes a pattern of mid-afternoon upset followed by recovery is indicated (for example, 
on 13 August 1986). This reflects the fact that phosphorus concentrations are one of the 
effluent concentrations specified in the permit. In addition, the within-day variations shown 
here indicate that grab sampling of the effluent is of only limited usefulness in determining 


loadings, as the result obtained may be dependent on the time of sampling. 


Flow in the Main STP effluent shows no general pattern of variation with time of day 
(Figure 3.7). Sometimes flow increases over the course of the day (for example, on 10 
September 1986, and 13 August 1986). There is also considerable variation between days 
(e.g. 800 to 1,200 x 10°m*/d on 13 August and 1,400 to 2,100 x 10°m°/d on 11 September 
1986). Stormwater is included in the Main STP effluent; consequently, flow is greater after 
a rainfall. There were large rainfall events on 10 and 11 September. 


Instantaneous loadings were computed over the course of the day using flow data from the 
Metro Works Department (readings every 15 minutes). These flows (m*/d) were multiplied 
by 1986 concentration data from the Main STP study, matched to the nearest 15 minutes. 


The ammonia loadings (Figure 3.5) are less consistent in pattern than the ammonia 
concentrations. Generally, ammonia loadings tend to increase over the course of the day, 
following the concentration pattern (for example, on 11 September 1986). However, on 
other days, the loading pattern fluctuates without reflecting the concentration increase over 
the course of the day (for example, on 13 August 1986). 


JSQUISAON /0 °6 Jaquisidas oT °9 ysn3ny HI °¢ 
JOQqW3AON 90 °8 Jaquiaidss 690 «°C isn3ny ¢] °z 
Jaqwaidss TT °Z Jaqwaidas 390° ysn3ny Z] oT 

(a}dwes Jo a}3ep ayedIpul sjoquiAs jo]d) 


9861 dLS NIVW NI ONITdWYS dO AWIL HLIM SNIGVOT GNV NOILVULNADNOO SNUYOHdSOHd TVLOL LNAN1SS4 29°€ FJYaNdIS 


SJy owt) : 544 SW] 

ve S°dt cI aA | OT c‘l S c‘z 0 02 a a SI ar | Ot cz 5 2 

GL'3] Gz°9T cL°ey fae 8 c.°8 oz°9 cL’ 62°} SL’8l 2°91 cL’et rane o2°8 s°9 g2°E 37°] 
aa aa Salah Salted Dalal etal telat tee Dee et Set Soot Cnet Cees a ee oY thewnn torn mnbonmn pomnn panne panna panna panna pon fan pan nn pn nn pom n pnb - = bg 
, +0 + he 
H sy $25 $5 § = S. § ! i let te eP te te th ‘ 
H Pe SE Va 9 GSe>- $26 ! i yoy 3 8 I 
H J Canc mC REE, : ‘ I Beek z 2 
H t Zz f H z cae 
$ I ‘ +62" + zt rossi 
H fh : ' et H 
- z l : H ' 
H rs ! H Eee , 
H : H is § it 
, J +6°2 ¢ @ 2 rise 
' td ! ' Sey, 
$ : ce 
i é z F ! H 
H E ! 6 ' 
3 g 2 a + a 
: u i : 
i ' fe) t t i 
' ' { H 
d 46 3 z $522 
' 9 ' e i = 
H OS ! 4 ' S : . 
H ' 5 H t 
i H u ' et ‘ S 
¢ £ +798 + 9 {ime Z 
H ' bt) i 9 3 
' 3 e : u ' * 
‘ f ° i i 
H ' R) i 
: +6°L il bens 
H ! ” 9 
J +61 "8 ¢ si877) 
: ! ' t 
+ +0] + . 
$b--~-4----4----4----5----4----4---- Pieces Pei einsise' iene Sasi ox Seeaatumsatn-aabLaasetoase ys th al a le th oa ESC I Oo Soe Sy Ree ee eS, ee ee 


+ + + 
WANT JF g1S UteN IZ UOTzERS WEN JIT 31S -12W 27 uCtyes: 


35-129 


Station 22 Main STP Effiuent 
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The loading pattern for phosphorus (Figure 3.6) more closely resembles the concentration 
pattern. There is a tendency for upset loadings to occur later in the day, although loadings 
at most times are near the low end of the scale, even late in the day. There is also 
evidence of recovery following upset (for example, on 13 August 1986). These patterns 
probably reflect effort to control loadings on the basis of phosphorus. 


Daily variation in concentrations of ammonia and phosphorus in the lake near the outfall 
(Station 1419) do not closely follow either each other, or the corresponding concentrations 
in the STP. Figure 3.8 illustrates daily variation in both parameters at the outfall in 1986. 
Figure 3.9 compares the STP and station 1419 concentration patterns for ammonia over a - 
three-day period. This lack of correlation suggests a strong influence of receiving water 
conditions on outfall concentration. 


3.1.4 Interparameter Correlations 


Correlations between chemical parameters in the effluent may be used to identify potential 
surrogate parameters for trace organics or certain metals (e.g., cadmium) that are seldom 
detectable in the receiving water. Dispersion of the surrogate can be used to infer the 
dispersion of the trace contaminant. However, the dispersion properties of the surrogate 
in receiving water should be shown to resemble those of the trace contaminant. This can 
be tested if the trace contaminant is sometimes detected in receiving water over a 
meaningful concentration range. Receiving water concentrations above the detection limit 
can then be checked for correlation with the surrogate. 


(a) Parameters Measured at Conventional Detection Levels 


Interparameter correlations within the STP effluent are summarized in Table 3.12. Only 
measurements above the detection limit contribute to these correlation coefficients. The 
large positive coefficients for suspended solids with all metals suggest that suspended solids 
may be a useful surrogate for metals. The largest coefficients are for the metals of greatest 
concern with respect to potential toxicity (Cd, Cr, Cu, Hg, Pb, Zn). Total phosphorus is 
almost as good, and in fact, superior as a surrogate for nickel. DOC is the best predictor 
for nickel and Pb, though Pb sample sizes were small. 
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The only organic parameter measured above detection limits with sufficient frequency to 
compute a meaningful correlation to conventional parameters was lindane (7-BHC). No 
surrogate for lindane was evident from the correlation coefficients. Asignificant correlation 
(r = 0.54) of «-BHC with conductivity and total solids was suggested, although the sample 
size (n = 12) was small. 


The Metro Works Department measures a few conventional chemical parameters in daily 
composite samples of STP effluent. The 1986-1987 data, over periods coincident with the 
Main STP Study sampling, were used to examine correlations among conventional 
parameters. Table 3.13 shows these interparameter correlation coefficients. The 
correlations between ammonia and TKN (r = 0.88) and total phosphorus and suspended 
solids (r = 0.85) are significant (p<0.05). The first coefficient reflects the fact that TKN 
constitutes organic N+ ammonia and, in sewage, is approximately 90% ammonia. The 
second coefficient explains the similarity of surrogate relationships for phosphorus and 
suspended solids in Main STP effluent (Table 3.12). 


b) Ultra-trace Organic Parameters 


The improved detection limits obtained using the 16-L STP samples (Section 2.0; Eli-Eco 
Laboratories, 1988) permitted an improved search for surrogates of toxic organics which 
might be especially useful in modelling of the toxics. A subset file containing data for only 
those dates in 1987 for which ultra-trace samples were collected was prepared, and the 
following were selected as potential surrogate parameters: ammonia-N, total Kjeldahl N, 
nitrate + nitrite-N, total P, pH, conductivity, chemical oxygen demand, dissolved organic 
carbon, total filtered solids and suspended solids. Computer files containing paired 
observations of conventional and ultra-trace organic data were prepared, and normal 
distribution plots were prepared using STATGRAPHICS in order to determine the best 
transformation for approaching normal distribution. In most cases, logarithmic trans- 
formation was used. Correlations were then computed using the SPSSPC statistical 
package. 


Initial attempts at correlation revealed many significant negative correlations, and scatter- 
plots revealed that several outlier points were influencing the correlations. Examinations 
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TABLE 3.13: CORRELATIONS AMONG CONVENTIONAL CHEMICAL PARAMETERS 
IN MAIN STP EFFLUENT, 1987 (Metro Works Data) 


Total 
Chemical Total Suspended Total Fecal 
Parameter Ammonia TKN Solids Phosphorus pH Coliforms 
Total Residual -0.2894* -0.2954* 0.0650 0.0411 OF13 77 -0.5466* 
Chlorine (57) (57) (104) (105) (107) (46) 
Total Ammonia 0.8836* -0.1141 -0.1667 0.0388 0.2224 
(57) (56) (56) (57) (31) 
TKN 0.1678 0.1386 -0.1618 0.2853 
(56) (56) (57) (31) 
Total Suspended 0.8478* -0.1708* -0.0158 
Solids (103) (104) (45) 
Total Phosphorus -0.2168* 0.1294 
(104) (45) 
pH -0.2047 
(46) 


Fecal Coliforms 


(n) indicates number of samples with paired measurements above detection limits 
contributing to the correlation coefficient. 


* indicates statistical significance based on one-tailed test (p<0.05). 


TABLE 3.14: 


Correlations: 


ORGANIC PARAMETERS, MAIN STP, 1987 


DcB13 


DCB14 


DCB12 


CORRELATIONS BETWEEN CONVENTIONAL AND ULTRA-TRACE 


TECB1235 


TECB1245 


TECB1234 


a 


NNHTFR 


NNOTFR 


NNTKUR 


PPUT 


PH 


coo 


COND25 


DOC 


RSF 


RSP 


5741 
« 11) 
P=.065 


ei Sola 
« 11) 
P=.010 


-6072* 
Cth) 
P=.048 


-5180 
« 11) 
P=.103 


- 1114 
« 11) 
P=.744 


-5604 
¢ 9) 
P=.117 


.8777** 
Comet) 
P=.000 


-9305** 
( 11) 
P=.000 


-7169* 
Coetk) 
P=.013 


-4426 
« 11) 
P=.173 


.2510 
Creat) 
p=.457 


- .0730 
« 11) 
P=.831 


- .3038 
« 11) 
P=.364 


-1515 
« 11) 
P=.657 


-2040 
( 11) 
P=.547 


- 3794 
os) 
P=.314 


- .0089 
( 11) 
P=.979 


- 1458 
« 11) 
P=.669 


- 2400 
« 11) 
P=.477 


-3116 
( 11) 
P=.351 


-1717 
Gif) 
P=.614 


- .6074* 
¢ 11) 
P=.047 


-3232 
Koa bh) 
P=.332 


-9800** 
« 11) 
P=.000 


-.7322* 
Crt) 
P=.010 


-8267** 
( 9) 
P=.006 


-6362* 
( 11) 
P=.035 


-7104* 
« 11) 
P=.014 


-2582 
¢ 11) 
P=.443 


-8966** 
(qos fel 
P=.000 


-8666** 
« 11) 
P=.001 


2441 
GFF 11) 
P=.469 


-8395** 
¢ 11) 
P=.001 


-2305 


« 11) 
P=.495 


- 3011 
Care tt) 
P=.368 


-7189* 
¢ 9) 
P=.029 


3136 
( 11) 
P=.348 


-4728 
« 11) 
P=.142 


-0091 
@ ail) 
P=.979 


-2509 
« 11) 
P=.457 


-8369** 
« 11) 
P=.001 


- 1654 
« 11) 
P=.627 


-8037** 
(11) 
P=.003 


-2691 
« 11) 
P=.424 


-.2570 
« 11) 
P=.446 


-6763* 
( 9) 
P=.045 


-3658 
Coals) 
P=.269 


£5441 
Goll) 
P=.084 


- 1049 
(11) 
P=.759 


-3114 
¢ 1) 
P=.351 


-6837* 
« 11) 
P=.020 


0936 
« 11) 
P=. 784 


- 7418** 
« 11) 
P=.009 


-5881 
CRT) 
P=.057 


- 6013 
Gana WD 
P=.050 


.8606** 
( 9) 
P=.003 


-3827 
Gaal) 
P=.245 


9313 
(11) 
P=.093 


- .0240 
« 11) 
P=.944 


-5877 
( 11) 
P=.057 


eee 


Note: 


First number is correlation coefficient 
Bracketed number is no. of observations 
P = probability of random correlation. 


Organic parameters are as follows: 
DCB13 = 1,3-dichlorobenzene 

pcB14 = 1,4-dichlorobenzene 

DCB12 = 1,2-dichlorobenzene 
1,3,5-trichlorobenzene 
1,2,4-trichlorobenzene 
1,2,3-trichlorobenzene 


TCB135 
TCB124 
TCB123 


TECB1235 
TECB1245 
TECB1234 


* + p<Q. 
wk - p<Q. 


TCB135 TCB124 TCB123 
- .1932 -5048 -7342* 
© Ul) .Ge 11) Gee 17D 
P=.569 P=.113 P=.010 
- .2746 . 1980 - 1207 
« 11) ¢ 11) ¢ = 11) 
P=.414 P=.559 P=.724 
- .3018 -5523 -7111* 
Gan LD Gan 1) (Cong l) 
P=.367 P=.078 P=.014 
- 6008 -5283 -2945 
Gv. 11) COA.11) | (CURE) 
P=.051 P=.095 P=.379 
-9085** - .7429** -.1935 
« 11) Cytit) C7 iT) 
P=.000 P=.009 P=.569 
- 6314 . 7133* -6292 
¢ 9) ( 9 ¢ 9) 
P=.068 P=.031 P=.069 
.0923 - 1733 -4002 
« 11) GOS yh oC GO I 
P=.787 P=.610 P=.223 
- .0560 3121 5484 
eh aa hk ear Gem 
P=.870 P=.350 P=.081 
-5920 - .2992 -2171 
Ceti) Cali)» «11 
P=.055 P=.371 P=.521 
- 5098 -5066 - 4008 
« 11) € fate? G owht) 
P=.109 P=.112 P=.222 
= 1,2,3,5-tetrachlorobenzene 
= 1,2,4,5-tetrachlorobenzene 
= 1,2,3,4-tetrachlorobenzene 
05 
01 
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TABLE 3.14: (cont'd) 
acB HCB HCE a-BHC 7-BHC GCHLOR ACHLOR PCB 

- .7070* 5999 -.4540 -.1457 - .5662 -.0100 4275 3776 

« 11) ¢ 11) ¢ uw) 11) ( 11) « 11) « 11) « 11) 
P=.015 P=.051 P=. 161 P=.669 P=.069 P=.977 P=.190 P=.252 
-1750 -.5910 -.5071 - .6545* - .3619 -.1525 -.8957** .1226 

« 11) ¢ 11) ¢ 11) ¢ 11) ¢ 11) « 11) ¢ 11) ¢ 11) 
P=.607 P=.056 P=.111 P=.029 P=.274 P=.654 P=.000 P=.719 
-.8112** .7193**  -.4505 - .0256 - .4632 - 1589 5073 -4300 

¢ 11) CARAT a6 11) ¢ 11) ¢ 11) « 11) « 11) ¢ 11) 
P=.002 P=.013 P=.164 P=.940 P=.151 P=.641 P=.111 P=. 187 
-.7439** .8453** .0627 -0996 - .0408 -8190** -6211* -4203 

« 11) ¢ 11) ¢ 11) ¢ 11) ¢ 11) ¢ 11) « 11) Cott) 
P=.009 P=.001 P=.855 P=.771 P=.905 P=.002 P=.041 P=.198 
-1692** -.6095** .4446 - 1437 -0681 -.9105** -.2052 - 2695 
« 11) ¢ 11) ¢ 11) 4 11) Ge ati) « 11) ¢ 11) ¢ 11) 
P=.006 P=.047 P=.171 P=.673 P=.842 P=.000 P=.545 P=.423 
- .9776** .9505** -.5114 -0556 - 2664 -6147 -5542 5541 

( 9) ¢ 9 ¢ 9) ¢ 9) ( 9) ¢ 9) ( 9) ¢ 9) 
P=.000 P=.000 P=.159 P=.887 P=.488 P=.078 P=.122 P=.122 
- 5677 -8489** .2206 -4838 -.0025 -2081 95 24** 2240 

« 11) ¢ 11) ¢ 11) ( 11) ( 11) Gopal) (Cae 11) (oe-1) 
P=.068 P=.001 P=.515 P=.132 P=.994 P=.539 P=.000 P=.508 
- .6664* -8932** .1729 - 1902 -.2919 2191 -9258** - 3435 
Catt) (Ge ay ke 11) ¢ 11) ¢ 11) ( 11) ¢ 11) « 11) 
P=.025 P=.000 P=.611 P=.575 P=.384 P=.517 P=.000 P=.301 
- .0353 - .4286 - .6084* - 5537 .0823 -.2560 -.7655** .0823 

« 11) ( 11) ¢ 11) ¢ 11) ¢ 11) ¢ 11) ¢ 11) Cait 
P=.918 P=. 188 P=.047 P=.077 P=.810 P=.447 P=.006. P=.810 
- .6475* .-7528** .1859 - .0359 -.1331 - 7328* 5387 -4863 

« 11) ( ty 9 £4 11) ¢ 11) ( 11) ¢ 11) ( 11) «= 11) 
P=.031 P=.007 P=.584 P=.917 P=.696 P=.010 P=.087 P=.129 
Note: First number is correlation coefficient 


Bracketed number is no. of observations 
P = probability of random correlation 


Organic parameters are as follows: 


QCB = PeuLecn ohovenzene 
HCB = hexachlorobenzene 
HCE = hexachloroethane 
ABHC = a-BHC 


GCHLOR = y-chlordane 
ACHLOR = a-chlordane 


PCB = Total PCBs 
GBHC = 7-BHC (lindane) 


of the raw data for conventional parameters in the Main STP effluent indicated excessively 
high (> 100 mg/L) suspended solids concentrations in the samples collected on June 2 and 
October 5. These samples also had abnormally high total phosphorus and COD 
concentrations. These two dates were excluded from the calculations; in addition the data 
for COD only on July 28 (trace to non-detect at all four STPs) were excluded and the 
correlations rerun. 


Interparameter correlations thus obtained are summarized in Table 3.14. The largest 
number of significant positive correlations (P<0.05) are found for ammonia, total Kjeldahl 
N and chemical oxygen demand. No more than three toxic parameters are correlated with 
any conventional parameters at the P<0.01 level. The low number of correlations is not 
surprising given the small number (N = 11) of paired observations; this number is even 
somewhat of an artifact as the single ultra-trace results were duplicated to pair with the 
replicate conventional results, and in fact only eight ultra-trace samples were actually used 
in the correlations. 


In contrast to metals, where suspended solids was found to be a surrogate for most 
parameters, no one parameter acts as a surrogate for all trace organics; in addition, no 
significant surrogate was found for several parameters (most notably PCBs). Besides the 
small number of samples available for correlation, the range of values for some parameters 
(e.g. conductivity) was limited. Short-term variability of effluent concentrations may be 
great, in particular the amount of suspended material, as was observed on June 2 and 
October 5. This may have a strong impact on organic concentrations, as the conventional 
and trace organic samples were not necessarily collected at exactly the same time and it is 
known that hydrophobic organic compounds may sorb strongly to particulate matter. 


In summary, the data collected from the Main STP were not adequate to establish 
correlations between conventional and ultra-trace compounds to a high level of confidence. 
It is recommended that a larger number of samples be obtained (including daily replicates 
for ultra-trace organic analyses); furthermore, conventional and ultra-trace samples should 
be collected together and extremes in suspended solids concentrations be avoided where 
possible as the heterogeneous nature of such samples renders it difficult to establish the 
representative nature of such samples. 


3.2 Lake Ontario, Main STP Ar 


Water quality data were collected at stations in the vicinity of the Main STP during 1986 
and 1987 (Figures 1.1 and 1.2) in support of the overall Main STP MISA water quality 
component described in this report, as well as the various toxicity tests described in other 
component reports. This included some mid-depth and bottom samples as well as surface 


samples; additional details of the sampling program are given in Section 2.0. 


The water quality summary for each station includes the mean, median, standard deviation, 
minimum and maximum concentrations of each measured chemical parameter. The 
maximum likelihood (MLE) method (El-Sharaawi, 1989; El-Sharaawi and Dolan, 1989) was 
used for parameters such as metals and organics, which are frequently below the analytical 
detection limit. Summaries for these parameters include the upper and lower confidence 
limits for the mean, as well as the frequency of samples below the detection limit, or above 
the provincial water quality objectives for protection of aquatic life (PWQOs). For residual 
chlorine, which decays rapidly in the lake, concentrations below the detection limit were 
considered to be zero. A comparison of several methods of handling "non-detect" results 
is given in Appendix 4. A separate summary is presented for surface and bottom samples 
in Appendix 1 (toxicity test periods) and Appendix 2 (annual statistics). 


In 1986, the water depth of "surface" samples varied between 0.1 and 1.0 m in an apparently 
random fashion. Profiling results (Section 3.2.5) indicates that significant depth variations 
do occur, with highest concentrations at the surface. However, it was felt that horizontal 
as well as short-term temporal variability was at least as great as depth variability in the top 
meter; consequently, surface data within this depth interval were pooled before calculation 


of statistical summaries. 
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3.2.1 Water Quality During Toxicity Tests 


Toxicity tests were performed in situ by the MOE at stations near the Main STP outfall. 
Rainbow trout were placed in fish cages to detect any acute mortality on 12-14 August, 
08-10 September and 04-07 November 1986; and on 15-19 June and 06-10 July 1987 (Flood, 
1990). In addition, clams were placed in cages to detect any chronic mortality or 
bioaccumulation on 16 June - 08 July and 13 October - 05 November 1987 (Aquatic 
Sciences, 1987). Water quality in the vicinity of Main STP during these periods is 
summarized briefly in this section. Detailed summary statistics for each station are included 
in Appendix 1. Results of these toxicity tests are reported in other Toronto MISA 
component reports (Flood, 1990; Hayton, 1990, Neville, 1990, etc.) 


During toxicity test periods, each sample was considered a separate unit of observation, 
even though some samples were field replicates. This has the effect of biasing mean 
concentrations toward dates and times of replication. However, these studies did not span 
long time periods. 


The chemical parameters which exceed PWQOs (MOE, 1984) near the Main STP outfall 
include ammonia (presented as a total concentration, but used as unionized ammonia in 
exceedence evaluation), phosphorus, chlorine, iron, copper, zinc and cadmium. Ammonia 
and phosphorus were usually in exceedence at the outfall (Station 1419) and often in 
exceedence near the outfall (Stations 1431 to 1434). The PWQO for both phosphorus and 
unionized ammonia is 0.02 mg/L. These nutrient concentrations were consistently greater 
near the surface than near the bottom, and tended to decline from August through 
November (Table 3.15). 


Total residual chlorine, as measured in the field, was consistently detected at a detection 
limit of 6 ug/L at the Main STP outfall, and was frequently detected within 500 m of the 
outfall. This parameter was always in exceedence of its PWQO (2 ug/L) when detected. 
However, a single measurement at a control station (20 ug/L at Station 2029 in September 


1986) also exceeded this objective. When residual chlorine was not detected, the concen- 
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TABLE 3.15: WATER QUALITY MEANS AT STATION 1419 DURING TOXICITY TEST 
PERIODS (see Appendix 1 for complete statistical summary) 


mee ee a a a a ew ww a a a we we ae a a a a a a = 


12-14 08-10 04-07 15-19 16 June- 

Parameter? Aug 86 Sept 86 Nov 86 June 87 08 July 87 
Surface (0-1 m depth) 

Sampling Time (hrs) 14.2 1249 13.6 11.0 Live? 
Temperature (°C) 14.6 11s ab Li6 12.9 13.0 
Conductivity (pmho/cm) - 547 502 - - 
pH 7.88 - 7.84 dun? Eat 
Residual Chlorine(pg/L) 551 1028 - 513 496 
Field Ammonia 7.30 Gras 5.17 5.46 4.91 
Lab Ammonia 6.79 6.17 4.55 - - 
Phosphorus 0.296 0.116 0.078 0.14 0.10 
Iron Ones 0.143 0.434 0.284 0.277 
Copper 0.007 0.004 0.004 0.005 0.005 
Zinc 0.017 0.012 0.008 0.006 0.006 
Turbidity (FTU) 5.94 2.45 9.5 3.431 3.23 
Chlorophyll a (ug/L) - Lee 29 - - 


Bottom or Depth of Turbidity Max. (4-6 m) 


Sampling Time (hrs) 14.2 12.9 13.6 10.8 11.0 
Temperature (°C) Lae? 9.2 8.8 8.5 8.9 
Conductivity (pmho/cm) - 367 322 - - 

pH 7.96 - 8.22 7.8 738 
Residual Chlorine(pg/L) 377 137 - 0 0 
Field Ammonia 4.40 1370 0.91 0.29 0.36 
Lab Ammonia 4.16 1.42 On - - 
Phosphorus 0.154 0.069 0.046 0.42 0.04 
Iron 0.208 0.080 0.581 0.182 0.191 
Copper 0.006 0.002 0.003 0.003 0.001 
Zinc 0.016 0.003 0.006 0.004 0.002 
Turbidity (FTU) 4.16 1.8 Laz 2.60 2.60 
Chlorophyll a (ug/L) . 35 3.0 - - 


1 Units are mg/L unless otherwise indicated. 
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tration was considered to be zero. It should also be noted that the total residual chlorine 
result, when measured in the presence of ammonia (as is the case here) includes mono- and 


di-chloramines as well as free chlorine. 


Heavy metals (iron, copper, zinc and cadmium) exceeded PWQOs occasionally at Station 
1419 during toxicity test periods. Cadmium was usually either non-detectable or in 
exceedence, since its PWQO (0.2 ng/L) was equal to the detection limit. The PWQOs for 
iron, copper and zinc are 0.3, 0.005 and 0.03 mg/L, respectively. Heavy metals were not 
present in receiving water at levels likely to have caused acute toxicity (MOE, 1990). 


The chemical parameters most likely to be responsible for toxic responses near the Main 
STP outfall were residual chlorine and/or ammonia. Acute toxicity to rainbow trout may 
be expected at approximately 50 to 200 ug/L free chlorine (U.S. EPA, 1984) or 10 mg/L 
total ammonia. Ammonia toxicity depends on pH and temperature which dictate the 
fraction of ammonia in unionized form (Emerson et al. 1975; MOE, 1984): 


1 
Fraction Unionized Ammonia = 1 + 19(PKa-pH) 


where: pKa = 0.09 + 2,730 
| 273 + Temperature (°C) 


Residual chlorine concentrations appear to be closer to acute toxic levels in general, and 
often considerably in excess of acute toxic thresholds for fish. 
3.2.2 Annual Summaries, Concentration Isopleths and 

Dilution Factors 


The annual water quality summaries for stations in the vicinity of the Main STP in 1986 and 
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1987 follow the same format as the summaries for toxicity test periods. They include the 
concentration mean, median, standard deviation, minimum and maximum, and frequen-cies 
below detection limits or above PWQOs (Appendix 2). However, in this case, N represents 
the number of sampling days and the daily mean is the unit of observation. This was 
considered preferable for representation of general spatial/temporal pattern since it prevents 
bias due to weighting of local means towards dates of greatest replication; however, standard 
deviations and exceedences may not reflect the full range of diurnal variation in 1986 when 
several samples per day were sometimes taken. Replication was not performed in 1987, and 


stations were seldom sampled more than once per day. 


Annual station means and standard deviations were plotted on maps of the Main STP area 
in order to illustrate dilution patterns in the vicinity of the Main STP outfall. Approximate 
concentration isopleths in 1987 were estimated from plots of concentration vs. distance from 
source on northeast and southwest transects. The form and extent of the concentration 
contour in the southeasterly direction is based entirely on speculation, since no sampling was 
performed in this direction. Rapid dilution in this direction might be expected since currents 


generally run toward the northeast or southwest (see Section 4.3). 


Concentration isopleths were also expressed as dilution isopleths (i.e. 1: x dilution) as follows: 


Xe= Source - Background 


Concentration - Background 


where: Concentration is estimated from the plot described above. 


* 


Effluent concentrations in the STP were used as source concentrations. Background was 
estimated from the control station means. The value of X increases with dilution. 
Concentrations can also be expressed as a relative concentration equal to 1/X which 
decreases with dilution. The X values at some individual stations, based on observed station 
means, are listed in Table 3.16. 
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TABLE 3.16: DILUTION FACTORS! FOR DIFFERENT CHEMICALS AT SELECTED 
INDIVIDUAL STATIONS IN 1987 


Ghemicalmee sae mL mele 27 + ee eee he ck 
Parameter? 1419 2877 2909 2878 2910 2209 


Field Ammonia: 


0.5 m depth Shad S26 52% 14. 45. 80. 

6.0 m depth 80. 300. 300. 45. 300. - 
Lab Ammonia Phe) 40. 110. 8.7 19. LL Oe 
Residual Chlorine? 3.0 2133 28. - - 400. 
Phosphorus 5.0 25 40 Die 40 80 
Conductivity 3.0 18 22 14. 40 100 
DOC Maen 20 27 16. 50 130 
Fecal Coliforms 6.6 470. 700. 27. 50. 700. 
Copper 5.4 36 100 208 29 150 
Zinc Tol she 150 75. 70 400 
Iron 4.1 20 28 16. 40 30 
Manganese 3.2 Li 28 12 24 75 


a a a a cm a ae ne a a a er ARTI RE Oe! EET: 


1 Dilution Factor = (Source-Background) / (Concentration- Background) 


2 Chemical measured at 0.5 m water depth in receiving water, unless 
otherwise indicated. 


3 Source concentration of free chlorine assumed = 1000 pg/L; other 


concentrations averaged over period of Main STP chlorination 
(June-September). 
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The spatial distributions of ammonia (field and lab measured) and total phosphorus are 
llustrated in Figure 3.10. Concentration isopleths of field-measured ammonia at 0.5 m 
depth suggest dilution to approximately 0.5 mg/L within about 1 km from the outfall. 
Using lab-measured ammonia as a source concentration, 1:40 dilution is found near the 
headland 1.5 km southwest of the outfall. At a 6 m water depth, dilution to 0.1 mg/L is 
suggested within 1 km of the outfall. Acute toxic levels depend on temperature and pH, 
but are likely in the 10 mg/L range (MOE, 1990). While Flood (1990) did not find acute 
lethality to rainbow trout during exposures at times when the effluent was nonchlorinated, 
Neville (1990) found sublethal effects (growth reduction) at station 1419 in such times. 


Lab-measured ammonia at 0.5 m depth is diluted to approximately 0.1 mg/L (1: 200 
dilution) within 1 km of the outfall. Contours extend slightly toward the northeast in 
contrast to the field ammonia pattern. Lab measurements suggest slightly more rapid 
dilution than field measurements, possibly due to high bias in field measurements, or low 


bias in lake measurements, at low concentrations. 


pH is slightly depressed near the Main STP outfall. The mean at Station 1419 is 7.30.3 
units. A range of 6.5 to 8.5 is considered to be within PWQO. Average pH in Lake 
Ontario near Toronto is around 8.5. pH is probably more important as a modifier of 
ammonia toxicity than for any direct toxic effects. As pH decreases, the fraction of 


unionized ammonia (and the ammonia toxicity) also decreases. 


Phosphorus concentration isopleths at 0.5 m depth extend slightly toward the southwest, as 
for most other chemical parameters. An unusually high mean at Station 2910 was ignored 
in isopleth calculation. The data suggest dilution to 0.03 mg/L (1:20 dilution) at a distance 
of 2 to 3 km from the outfall. 
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Figure 3.10 


Spatial distribution of ammonia and phosphorus in 1987 


1:21 1 lecpleth for 1 mg/L 
) Concentration = 


Lab Ammonia (mg/L) at 0.5m depth In 1987 


Figure 3.11 shows the spatial distributions of residual chlorine, fecal coliforms, suspended 
solids and conductivity at 0.5 m depth in 1987. Isopleths of chlorine residual were 
determined on the assumption of zero free chlorine in when none was detected; however, 
October and November data were excluded from the calculation since STP chlorination 
ceases in October. The data suggest chlorine as high as 330 ug/L near the STP outfall, 
declining to about 40 ng/L 1 km south of the outfall (station 2909). Acute toxic levels for 
rainbow trout are approximately 80 to 150 ug/L (MOE, 1990); Flood (1990) found acute 
lethality to rainbow trout when caged at locations within 500 m of the outfall. Acute 
lethality to fathead minnows was also found with station 1419 water (Neville, 1990). 


Fecal coliform isopleths at 0.5 m depth extend slightly to the northeast, as for lab-measured 
ammonia. Concentrations are expressed as logs; however, dilution factors are computed 
using geometric means in Figure 3.11. Dilution occurs to approximately ten organisms per 
100 mL within 1 to 2 km of the outfall (1:300 dilution). 


Suspended solids show very symmetric concentration isopleths, with dilution to about 2 
mg/L within 1 km of the outfall (a dilution factor of 1:25). This dilution rate is similar to 
that of other chemical parameters. 


Conductivity is diluted essentially to background levels within 1 km of the outfall. As for 
suspended solids, the dilution contours are quite symmetrical, at least in the vicinity of the 
STP outfall. The high natural background may obscure current and wind effects further 
away from the outfall. 


A few heavy metals are sufficiently detectable to examine dilution patterns. Average 
concentrations have been estimated using the maximum likelihood estimator (MLE; 
Appendix 2 and 4) (Figure 3.12). Copper concentration isopleths at 0.5 m depth suggest that 
copper is diluted to background levels (0.001 mg/L) within 1 km of the outfall. Zinc is 
diluted to approximately 0.002 mg/L within 1 km of the outfall (1:100 dilution), slightly 
above the 0.001 mg/L background. Iron is diluted to approximately 0.06 mg/L within 


3 7- 48 


Figure 3.11 


Spatial distribution of several conventionals in 1987 


1.14 \ 3 leopleth for 3 Lt 
oe 
2 1:14 dilution 


Suspended solids (mg/L) at 0.5m depth In 1987 Conductivity (umhos/cm) at 0.5m depth In 1987 


Figure 3.12 


Spatial distribution of several heavy metals in 1987 


Manganese (mg/L) at 0.5m depth In 1987 


1 


0.1 leopleth for 0. 
Concentration = 
1:25 dilution 


1:26 ) 
Ir 


0.5m depth In 1987 


Iron (mg/L) at 


Zinc (mg/L) at 0.5m depth in 1987 


Copper (mg/L) at 0.5m depth in 1987 


about 1.5 to 2 km of the outfall (1:50 dilution), slightly above the 0.03 mg/L background. 
Manganese is diluted to approximately 0.006 mg/L within 1 km of the outfall (1:20 
dilution). These metals do not approach acute toxic levels, even at the Main STP outfall. 


Concentration maps of these chemical parameters for 1986 do not permit isopleth 
construction since sampling in the Main STP area was confined to the immediate vicinity 
of the outfall (Stations 1419, 1431, 1432, 1433 and 1434). Maps illustrating those data, 
and computed X values, are included in Appendix 3. Generally, sampling was later in the 
day in 1986 (1300. to 1500 hrs.) Total ammonia concentrations at Station 1419 were similar 
to those in 1987 at 0.5 m (approximately 7 mg/L), but higher than in 1987 at 6 m (1.9 vs. 
0.28 mg/L measured in the field). Phosphorus concen-trations were essentially the same 
(0.19 vs. 0.18 mg/L at 0.5 m). Residual chlorine concentrations were nominally higher in 
1986 than 1987 (973 vs. 332 ug/L at 0.5 m); however, the difference was not statistically 
significant). While the main contributors to toxicity (chlorine and ammonia) were found at 


slightly higher average con-centrations in 1986, this may be due to sampling later in the day. 


3.2.3 Plume Tracking 


Surface movement of the Main STP outfall plume was tracked during the summers of 1984 
and 1986 by Gore and Storrie using multiple drogue releases (Gore and Storrie, 1986b). 
These drogue trackings show the movement of water particles near the water surface and 
permit the calculation of dispersion characteristics at the existing outfall site as well as 
potential locations for the proposed outfall. It was found that for the existing outfall 
location, the predominant surface movement was northeastward (towards the Eastern 
Beaches), although southward movement parallel to the shore of the Eastern Headland was 
frequently observed. At the proposed offshore locations, the surface movement towards the 
Eastern Beaches was much less and the offshore movement component was much larger 
than at the existing location. The offshore location was also characterized by better 
dispersion than the present location. Further details of the studies are available elsewhere 
(Gore and Storrie, 1986b). 
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In order to map the extent of the Main STP plume in greater detail than that obtained by 
sampling at fixed locations, a plume tracking survey was undertaken for four days in June 
1986. The plume- tracking vessel was equipped with a submersible pump attached to a 
towed "fish". This enabled a continuous record of temperature, conductivity and 
fluorescence (calibrated to represent turbidity) to be measured by on-board sensors and 
stored on an on-board computer. In addition, frequent water samples obtained from the 
submersible pump were analyzed on board ship for ammonia and chlorine residual. Vessel 
position was monitored by Loran, the readings of which were also stored on the on-board 


computer. 


Of the four survey days, June 18 and 19 have been chosen for discussion here. June 18 
represents stagnation conditions (average currents = 1.3 cm/s towards the north as 
measured by a drogue launched near the outfall and tracked during the survey). June 19 
represents jetting conditions (average currents = 12 to 21 cm/s to the northeast as 
measured by a current meter located about 4 km southeast of the outfall). Spatial contour 
maps of conductivity, ammonia and chlorine residual on the two survey days are presented 
in Figures 3.13 - 3.15. In these figures, the marker at the southern tip of the east arm of 
the headland is used as a reference location for distance measurements. Contour maps for 


June 16 and 17 are presented in Appendix S. 


Concentrations of the three parameters contoured here in the Main STP effluent were 
similar on June 18 and 19; hence average values are used here in the discussion 
(conductivity = 1015 ymho/cm; ammonia = 21.0 mg/L; chlorine residual measured in the 
effluent line before discharge to the lake = 1050 pg/L). Average background conductivity 
was found to be 313 uymho/cm; background ammonia and chlorine residual were both 
assumed to be zero. Using these average values, the dilution factors (X in section 3.2.2) 
for the contour values used in the figures are as follows: 
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Figure 3.13b 


Figure 3.13a 
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Figure 3.14a 
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Conductivity Ammonia Chlorine residual 


(umho/cm) (mg/L) (ug/L) 
Contour Contour Contour 
Value X Value X Value 4 
320 100. 0.5 “2s 10 100 
350 19. 0 was: 20 53 
380 10.5 Lad 14. 50 21 
410 TE! 220 10.5 100 10.5 
440 ca haa 2 Now 8.4 200 553 
470 NAS) 320 7.0 
500 She fs! 37 6.0 
4.0 5e3 
4.5 47 
520 aul 


Under stagnation on June 18, dilutions for conductivity and ammonia increase to 1:100 
within 1 km of the outfall close to shore and 0.5 km of the outfall away from shore. As a 
consequence of its rapid decay, dilutions for chlorine residual (Figure 3.15a) reach 1:100 
within about 0.5 km to the northeast and 0.3 km to the southeast of the outfall. A mild 
elevation of conductivity towards Coatsworth Cut (X = 5.5 to 7.3) suggests entrapment of 
a portion of the plume in this area, or possible input from Coatsworth Cut itself. 


Under jetting conditions towards the northeast (Figures 3.13b - 3.15b), a plume lying a few 
hundred metres offshore from the Eastern Beaches was observed. For conductivity and 
ammonia, X values of about 20 extend for 5 to 5.5 km from the outfall (and including the 
area of the R.C. Harris WTP intake); 1:10 dilution is observed up to 3 km from the outfall, 
although its zone is somewhat patchy in nature between 1.5 and 3 km. This may reflect 
either short-term variation in outfall discharge and lake currents, or limitations associated 
with the area surveys, or both. As expected, the chlorine residual plume is again shorter 
(1:20 dilution at 1 km and 1:100 at 3 to 4.5 km from the outfall; figure 3.15b). 


These results provide a greater detail of resolution over short-interval time scales (order 
of hours) than the averaged results presented in section 3.2.3. In addition, they show the 


variability in plume size resulting from extremes of current regimes. (See also results for 
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June 16 and 17, appendix 5). The averaged results for conductivity (figure 3.11d) and 
ammonia (figure 3.10a-c) show zones for 1:20 and 1:100 dilution similar to or slightly larger 
than those observed under stagnation (figures 3.13a - 3.15a); however, due to the 
combination of northeastward and southwestward current conditions, the averaging 
procedure combines lower (downstream) with higher (upstream) X values. In particular, 
the reduced dilution factors resulting from jetting (figures 3.13b - 3.15b) which occurs only 
a small proportion of the time are completely obscured by the averaging process. 
Additionally, the results in this section provide improved definition of the rate of dilution 
in the offshore (southeast) direction. Dilution in this direction is generally good, with X 
values of 1:20 to 1:100 being observed within 0.5 km of the outfall on most survey days. 


3.2.4 Exceedence Frequencies 


Exceedence frequencies are included in Appendix 2 for all stations sampled in 1986 and 
1987. Tables 3.17 and 3.18 present a summary of exceedence frequencies at Station 1419, 
the Main STP outfall, at 0.5 and 6 m water depths. All chemical parameters which 


exceeded their PWQOs in receiving water are included. 


1987 


The PWQO for unionized ammonia (0.02 mg/L) was exceeded in 75% of 0.5 m samples 
over the Main STP outfall (Station 1419) in 1987 (Figure 3.16). This percentage is based 
only on samples for which pH and temperature were also determined, since both these 
parameters are needed to calculate unionized ammonia from the total measured value. 
Percent exceedence declines with distance from the outfall and is zero at both control 
stations (1536 and 2029). Exceedence isopleths (Figure 3.16) suggest that less than 50% 
of samples exceeded the PWQO at any point greater than 1 km from the outfall. No 
ammonia exceedences were detected in samples from a 6 m depth at the outfall, although 


only three samples had concurrent ammonia, pH and temperature measurements. 
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TABLE 3.17: EXCEEDENCE OF CONVENTIONAL AND METAL WATER QUALITY 
OBJECTIVES AT STATION 1419 OVER THE MAIN STP OUTFALL 


(see Appendix 2 for complete statistical summary) 


1986 1987 
Parameter N % N % PwQo? 
0.5 m Depth 
Residual Chlorine*** /7 100 13 70 2.0 pg/L 
Field Ammonia 6 100 i ifs) 0.02* 
Lab Ammonia 6 100 - - 0.02* 
Phosphorus 10 100 14 100 0.02 
Iron 10 40 13 di, 0.30 
Copper 10 30 S| 38 0.005 
Cadmium 10 20 13 31 0.2 pg/L 
Lead - - 13 0 0.025** 
Nickel 10 0 13 0 0.025 
Mercury 9 is 13 8 0.2 pg/L 
6 m Depth 
Residual Chlorine*** 6 100 5 0 2.0 pg/L 
Field Ammonia 6 67 8 0 0.02* 
Lab Ammonia 6 50 - - 0.02* 
Phosphorus 10 80 4 100 0.02 
Iron 9 22 4 25 0.30 
Copper 9 ie 4 ZO 0.005 
Cadmium 9 11 4 Zo 0.2 pg/L 
Lead - - 4 0 0.025** 
Nickel 9 0 4 25 0.025 
Mercury 8 0 4 0 0.2 pg/L 


1 PWQO in mg/L unless otherwise indicated; phosphorus value is a 


guideline only. 

* Ammonia objective is for unionized fraction. 

** Lead objective is alkalinity-dependent; alkalinity = 80 mg/L 
was assumed. 

*x* Chlorination period (to end of September) only. 
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TABLE 3.18: EXCEEDENCE OF PROVINCIAL WATER QUALITY OBJECTIVES FOR 
ORGANOCHLORINE PESTICIDES AT STATIONS 1419, 1536 AND 2029 


1986 1987 PWQO 
Parameter N % N % (ug/L) 
STATION 1419 
Surface (0 to 1 m Water Depth) 
Lindane 10 30 13 46 0.01 
Heptachlor 10 O* 13 8 0.001 
a-Chlordane 10 O* 13 0 0.06 
Dieldrin 10 10 13 O* 0.001 


Bottom (6 m or Depth of Turbidity Maximum) 


Lindane 9 11 5 0 0.01 
Heptachlor 9 O* 5 O* 0.001 
a-Chlordane 9 O* 5 O* 0.06 
Dieldrin 9 11 5 O* 0.001 
STATION 1536 | 
Surface (0 to 1 m Water Depth) 
Lindane 6 0 0 - 0.01 
Heptachlor 6 O* 0 - 0.001 
a-Chlordane 6 O* 0 - 0.06 
Dieldrin 6 O* 0 - 0.001 
STATION 2029 
Surface (0 to 1 m Water Depth) 
Lindane 6 0 0 - 0.01 
Heptachlor 6 O* 0 - 0.001 
a-Chlordane 6 O* 0 - 0.06 
Dieldrin 6 O* 0 - 0.001 


* Indicates chemical was never detected. 
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The phosphorus guideline (PWQG) for prevention of nuisance of algal growth (0.02 mg/L) 
was exceeded in all 0.5 m and 6 m samples at Station 1419 in 1987 (Figure 3.16). The 
frequency of exceedence in 0.5 m samples declined with distance from the outfall to 20% 
and 11% at control Stations 1536 and 2029, respectively. Exceedence isopleths (Figure 
3.16) suggest that less than 50% of samples should exceed the PWQG at any point more 
than 1 km from the outfall. Few 6 m samples were analyzed for phosphorus at stations 
other than 1419. 


Residual chlorine exceeded the PWOQO of 2 ug/L in all station 1419 samples collected 
between June and September of both 1986 and 1987. Exceedence frequency declined with 
distance from the outfall, to about 5% of samples more than 1 km south of the outfall. 
STP chlorination ceases at the end of September. The PWQO is designed for chronic 
protection of aquatic life. Acute toxicity to rainbow trout was found to occur within 100 
m from the outfall in as short a period as 4 h (Flood, 1990), and in 96 h at 500 m from the 
outfall. However, toxicity would not be expected beyond a radius of approximately 500 m 
from the outfall. | 


Industrial Waste Discharge Guidelines (MOE, 1987) were exceeded in 1987 only by 
cadmium in the Main STP effluent. However, most heavy metals (except mercury and 
chromium) exceeded PWQOs in the receiving water adjacent to the Main STP outfall 
(Station 1419) in 1987. Frequency of PWQO exceedence in 0.5 m receiving water samples 
was highest at the outfall Station 1419 for iron (77%), copper (38%) and cadmium (31%). 
Nickel exceeded its PWQO once only in a 6.0 m sample. In general, PWQO exceedences 
were rare and showed little spatial pattern beyond a 500 m radius from the outfall. Very 
few heavy metal exceedences (one each for iron and cadmium) were found in the clean 
water zone of the waterfront (Cluster 4, Section 3.3.1). None of the metal exceedences 


were in the acutely toxic concentration range. 
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1986 


In 1986, exceedence frequencies showed the same general pattern as in 1987, although 
sampling in the Main STP area was confined to the immediate vicinity of the outfall 
(Station 1419). Residual chlorine objectives were consistently exceeded near the outfall; 
a single value at control Station 2029 (near the R.C. Harris WTP intake) was also in 
exceedence. The ammonia objective and phosphorus guideline for prevention of nuisance 
algal growth were usually exceeded near the outfall, with higher exceedence frequencies at 
surface samples than at 6 m water depths. Forty percent of the iron values were in 
exceedence at the surface, with less frequent exceedence at 6 m. Copper and cadmium 
values were occasionally in exceedence at the outfall. Main STP effluent occasionally 


exceeded the 1 mg/L objective for STP phosphorus. 


Organochlorine pesticides (lindane and heptachlor) occasionally exceeded PWQOs in 
receiving water at the Main STP outfall (Station 1419). Table 3.18 lists all the pesticides 
detected in receiving water, and their percent exceedence at Stations 1419, 1536 and 2029 
(the R.C. Harris control station). Exceedence frequencies were greater in 1987 than in 
1986, and greater in surface samples (0 to 1 m) than bottom samples (6 m or depth of 
turbidity maximum) at the outfall. Lindane exceedence in the STP effluent (Station 22) was 
86% in 1986 and 100% in 1987. Heptachlor exceedence at 0.5 m in the effluent was 29% 


in 1987. There were no heptachlor exceedences in 1986. 


3.2.5 Depth Profiles 


The results presented for the plume area in the previous two sections for water depths of 
0.5 m and 6 m suggest consistently better water quality near the bottom than at the surface. 
In order to determine whether there is any significant subsurface plume (maximum 
concentration below 0.5 m), the monthly mean conductivities were plotted against water 


depth to give a depth profile for each contaminant at each station. 
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The depth profile for conductivity at the Main STP outfall (Station 1419) shows a distinct 
surface maximum at all times of year (Figure 3.17), while conductivity at the control 
Stations 1536 and 2029 shows no consistent pattern. The stations immediately surrounding 
the STP outfall (1431, 1432, 1433, 1434) tend to show a surface maximum, at most times, 
and thus resemble the outfall (Figure 3.18). 


In 1987, the depth profiles for conductivity at the Main STP outfall (Station 1419) and the 
Humber Bay STP outfall (Station 2072) both displayed a distinct surface maximum (Figure 
3.19). Profiling at stations immediately surrounding the Main STP outfall was limited in 
1987, but the results obtained showed similar surface maximum to those in 1986 (Figure 
3.18). 


Other stations on the waterfront which were subjected to profile sampling failed to show 
this pattern, or any other consistent pattern. The STP effluent apparently rises quickly to 
the surface, since it is warmer than the receiving water at depth, but the plume loses any 
consistent depth dependence within a few hundred metres. 


3.2.6 Fugacity Analysis of Toxicity Potential 


As indicated in Sections 3.2.2 and 3.2.4, the major sources of toxicity in the vicinity of the 
Main STP outfall are free chlorine and ammonia, and toxic effects are likely to be confined 
to within a 500 m radius of the outfall. While there are exceedences of heavy metal and 
organochlorine pesticide water quality objectives at the outfall (Station 1419), concentrations 
of these contaminants in water do not approach toxic levels. 


Mackay et al. (1987) have devised an index (y) of toxicity potential suitable for application 
to non-polar organic contaminants: 
y= (K ®) C,, MV 10° 


where: K Aa octanol-water partition coefficient, 
C, = chemical concentration in water (g/m°), and 
MV = molar volume (cm’/g). 
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FIGURE 3.18 


Conductivity Depth Profiles at Stations 500 m 
From Main STP Outfall in 1986 
(plot symbols indicate month) 


Conductivity vs. Water Depth at Station 1431 in 1986 Conductivity vs. Water Depth at Station 1432 in 1986 
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Conductivity vs. Water Depth at Station 1433 in 1986 Conductivity vs. Water Depth at Station 1434 in 1986 
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Interpretation of this index is as follows: 


y greater than 0.01: acute effects virtually certain, 


0.01-0.001: acute effects possible, 

less than 0.001: acute effects probably absent, 

greater than 0.001: chronic effects certain, 

0.001-0.0001: chronic effects probable, 

0.0001-0.00001: chronic effects unlikely, but possible, 
and 


less than 0.00001: chronic effects almost certainly absent. 


We computed this index for lindane which was the only organic contaminant frequently in 
exceedence of the PWQO for receiving water at the Main STP outfall. Using the maximum 
concentration recorded at the outfall (Station 1419) in 1987: 


(2238.7)°* (0.000021 mg/L) (224) 10° 
0.000002 


< 
tt 


Thus, chronic toxic effects on biota would be almost certainly absent at the maximum 
concentration observed. Using the average concentration at the outfall (10.02 + 3.92 ng/L) 
in 1987: 

(2238.7)°* (0.000010 mg/L) (224) 10° 

0.000001 


< 
I 


This average concentration, equal to the PWQO, is an order of magnitude lower than the 


concentration at which chronic effects would be possible. 


Perusal of a limited quantity of APLE data also shows that lindane is the only organic 
contaminant present in the vicinity of the Main STP at levels approaching PWQOs for 
receiving water or other levels of concern to aquatic biota. Average concentrations at 
offshore stations were approximately 1 ng/L, or an order of magnitude lower than at the 
Main outfall. Maximum concentrations observed for Main and Humber STP effluents were 
56 and 66 ng/L, respectively. 


Spar 97, 


3.2.7 Interparameter Correlations With Main STP Effluent Data 


The obvious proximity of water quality impact zones (Section 3.3.1) to known point sources 
implies that the associated point source is responsible. However, quantitative prediction 
of lake concentra- tions of specific contaminants from the concentrations of those 
contaminants measured in the point source effluent is complicated by variation in source 
flow and loading, variation in lake current speed and direction, and a variable time lag 
between contaminant discharge and arrival at a given location in the lake. As a result, 
simple lake vs. source correlations or regression equations may have little predictive power. 


Lake vs. source correlation coefficients are tabulated in Table 3.19. These coefficients 
represent a simple linear model. The first row represents all lake stations, without a time 
lag, based on Main STP study data. Under these conditions, only conductivity was 
significantly correlated with the STP source (r = 0.59). There was a seasonal trend in 
conductivity on the waterfront (Section 3.3.4) which may be reflected, to some extent, in the 
Main STP. 


There are insufficient Main STP study data to examine correlations between the STP 
effluent and the outfall Station 1419 above (only two samples match to the nearest hour). 
However, the Metro Works Department daily composite samples were collected 
morefrequently and can be matched with daily average concentrations of contaminants at 
the outfall. In addition, daily loadings (average flow x composite concentration) can be 
matched with outfall concentrations. 


No source vs. outfall correlation was found for phosphorus, using daily average data, 
because there is much less day-to-day variation in phosphorus. 


Most of the phosphorus variation is related to fluctuating conditions of relatively short 
duration. Thus, source and outfall concentrations measured frequently and at matching 
times over the course of the day are needed to develop a predictive model. 
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Data suitable for development of a short-term predictive model are presently quite limited. 
To extend such a model beyond the outfall, appropriate time lags must be identified. A 


model based on available data is discussed in Section 4.2. 


3.2.8. hort-Term Correlation of Sour ding and Receivin 


Water Conductivi 


Correlating the instantaneous source loading and receiving water conductivity establishes 
whether receiving water conductivity is predictable over the short-term from instantaneous 
source loading. Given an approximately steady-state source, it is possible to correlate the 
source load to the receiving water conductivity after subtracting an appropriate time-lag 
from the receiving water time based on the average velocity from 1419 to 2891 and 1419 
to 2209 and the distance to a particular station. The conductivity data for receiving water 
in 1987 were thus shifted back in time and correlated to the source loadings at the same 


times (within a one-hour time window): 


Receiving Water Conductivity _ +; ,eMeasured-TimeLa g VS. 
Source Loady. oT 


The source load and receiving water conductivity were correlated for: 


o a data set extending south of the outfall (1419, 2877, 2909), and 
Oo a data set extending east of the outfall (1419, 2878, 2910). 


The expected result would be a positive correlation along each direction which decreases 
with increasing distance. The results obtained show a non-significant correlation over this 
entire spatial domain. We then reduced the domain to include only stations 1419, 2877 
and 2878. Again, there was no significant correlation. The correlations are negative for 


loading, indicating that as the loading increases the receiving water conductivity decreases. 


Je — w/o 


The correlations are larger in the eastern direction which is the direction of the most 
persistent current (i.e. largest mean vector component, Gore and Storie, 1986a). Finally, 
we examined the source vs. outfall correlation (no time lags involved), and obtained similar 


results. These results are summarized in Table 3.20. 


Since the source loadings are quite variable, a possible reason for the lack of source vs. 
receiving water correlation is the low frequency of sampling and the paucity of 
time-matched source and receiving water data. For example, in 1987 at the outfall itself, 
approximately 40 measurements were taken irregularly over a four-month period. Sampling 
rate should be directly proportional to temporal variation in source loading (i.e. 
concentrated at the most dynamic time of day), and sampling should include the entire daily 


range in order to optimize the study design for time series correlation. 


Since receiving water conditions (current, wind, weather) may vary considerably from 
day-to-day, the assumption of a constant time lag may also be inappropriate for the 
available data, which includes few samples on any one day. Intensive sampling over a short 
period of time (e.g., one day) would be more appropriate for development of short-term 


predictive models. 


Another possible reason for the lack of correlation at the outfall itself is that the outfall 
flow dominates the water volume in this vicinity. In such a system, receiving water 
concentrations should be correlated with source concentrations rather than loadings. 
Short-term prediction of receiving water concentrations from source concentrations is 
discussed in Section 3.2.9. 
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TABLE 3.20: SHORT-TERM CORRELATION OF MAIN STP EFFLUENT AND 
ADJACENT RECEIVING WATER CONDUCTIVITY, 1987 
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Receiving Water Time 

Data Set Lag (hrs.) Loading Conductivity 
Stn 1419 (outfall) 0 -0.93 (5)? -0.25 (7) 
Stn 1419-2877 (N-S) de -0.54 (8) 0.41 (10) 
Stn 1419-2878 (E-W) 1.8 -0.71 (7) 0.59 (9) 
Stn 1419-2909 (N-S) DG -0.34 (10) 0.36 (12) 
Stn 1419-2910 (E-W) ot) -0.51 (8) 0.47 (10) 


1 (n) indicates number of time-matched sample pairs within a one-hour 
time window. 
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3.2.9 Short-Term Correlation of Source and Receiving Water Conductivity 


Correlating the instantaneous source conductivity and receiving water conductivity 
establishes whether receiving water conductivity is predictable from source conductivity. As 
a stable chemical parameter, conductivity is more likely to be predictable than a chemical 
subject to degradation, volatilization or precipitation. The same methods described in 
Section 3.2.8 were used, except that source loading was replaced with source conductivity. 
Correlations were examined for: 


Receiving Water Conductivity; _ +: ,eMeasured-Tim eLag vs. Source 
Conductivityy: veT 


The results are included in Table 3.20. Again, no significant correlations were found either 


at the outfall, or over a broader spatial domain in the outfall vicinity. 


As outlined in Section 3.2.8, sample size were small, and allocation of sampling effort was 
sub-optimal for time series correlation. Moreover, the sampling effort captured only a small 
portion of the possible variation in conductivity. For example, neither the source or outfall 
water varied by more than 100 umhos/cm, in the few time-matched samples. However, the 
STP effluent varied by more than 300 umhos/cm over the long-term. 


SPS 


3.3 Lake Ontario, Toronto Waterfront Area 


In 1987, a fixed-grid survey encompassing the entire Metro Toronto area from Pickering to 
Lakeview was undertaken. The purpose of this survey is to provide further baseline data 
on water quality across the entire Metro Toronto Waterfront, including the impact on the 
lake of all pollution sources, not just the Main STP. Additionally, this survey provided an 
update on water quality zonation (previously obtained only for the central waterfront area 
- Poulton and Griffiths, 1986), as well as input data for the RAND currents and water 
quality model and TOXFATE fate and transport model of organic contaminants. In this 
section, the analysis of spatial zonation in water quality is discussed, including the average 


water quality of each zone and its variability. 


Data collected in 1986 are not considered here, as sampling was much more restricted to 
the immediate vicinity of the Main STP. 


3.3.1 Spatial Zonation 

The analysis of variance model first developed by El-Shaarawi and Kwiatkowski (1977) was 
used as a Starting point in the zonation analysis. This model breaks down the variance in 
the observation for any one chemical parameter (variable) into between- station and 
between-cruise portions. The following set of transformations (Box and Cox, 1964) are first 
applied to the data in order to satisfy the assumption of normal distribution required for 
the use of parametric methods such as analysis of variance: 


= A -l we 
Zi = (Yi; y/A (a = 0) 
Zi = 1h Yi (A = 0) 
where: Y 


ij = observed water quality at station j during cruise 


Z.. = transformed observation, and 


= transformation parameter. 


The method of maximum likelihood is used to estimate the value of 4 which produces the 
best approach to normal distribution of Zi for each parameter. 


The mean value of Zi is then resolved into additive components as follows: 


where: u = general mean, 
a, between-cruise (temporal) effect, and 
B, between-station (spatial) effect. 


For each parameter for which the station effects 8, are significant according to an F- test | 
with (j-1) and (ij-1) degrees of freedom, the 8, were used as water quality observations in 
multivariate analysis of water quality zones. Other parameters were excluded from the 


analysis. 


The water quality zonation analysis involved use of principal components analysis followed 
by cluster analysis, using a methodology similar to that first applied to Lake Ontario by 
Neilson and Stevens (1986) and Poulton and Griffiths (1986). 


Many of the chemical variables (station effects) are correlated with each other as they are 
often controlled by similar causative factors. Principal components (PC) analysis is used 
to transform such correlated data by finding linear additive functions of these variables 
which give a smaller number of new uncorrelated variables (PCs). Each of these PCs is a 
linear additive function of the original variables, with the first variable accounting for the 
largest proportion of variance in the original data, and the second variable accounting for 
the next largest amount of variance, and so forth. The PCs are derived so as to be 
independent of (ie. uncorrelated to) previously derived components; therefore, each 
component provides unique information about the original data set. They also tend to be 
more normally distributed than the original variables (station effects), although these station 
effects, as derived from El-Shaarawi and Kwiatkowski’s (1977) model, tend to be reasonably 


3°- 75 


normally distributed, provided an appropriate transformation index has been selected. 


Factor scores are calculated from the components for each station j and principal 


component k as follows: 


where: m = parameter index, 
Pp = total number of parameters, 
8. iss station effect in parameter m at station j (as calculated earlier), and 
cars = loading coefficient for parameter m on factor k. 


Cluster analysis was used to classify the stations into groups with similar factor scores from 


the principal components analysis. 


Alternatively, cluster analysis can be applied directly to the station effects; however, the 
result may be biased by the presence of highly correlated parameters as explained above. 
The method used was a hierarchical agglomerative method, i.e. it begins with each station 
occupying an individual cluster. At each step, the closest pair of clusters are combined into 
a single cluster; this proceeds until all the stations are in a small number of clusters which 


are readily interpretable. 


Squared Euclidean distance was used as a measure of dissimilarity between individual 
clusters. The distances between members of combined clusters may be evaluated as the 
minimum, average or maximum distances between individual members of the clusters being 
combined. These methods are referred to as single linkage, average linkage or complete 
linkage. 


The spatial distribution of station groups (i.e. water quality zonation) derived using 
complete linkage is shown in Figures 3.20 and 3.21. This clearly shows three groups of 
fig. 3.20 
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stations associated with known point sources. The Main STP outfall (Cluster 1) is quite 
distinct in terms of water quality. Other sources of moderate impact are represented by 
Cluster 2 (e.g. Humber Bay STP outfall). A zone of lesser impact (Cluster 3) is found 
beyond these moderate impact zones and in association with the Don River, Toronto 
Harbour and Rouge River. Cluster 4 represents relatively unimpacted stations. Box plots 


summarizing the water quality characteristics of each zone are shown in Figure 3.22. 


3.3.2 Ultra-Trace Organic Results 


The use of the APLE sampler during the three 1987 sampling periods represents the first 
comprehensive analysis of the concentrations of PBCs, organochlorine pesticides and 
chlorobenzenes undertaken in the Toronto Waterfront. Details of the sampling and 
analytical procedures are summarized in Section 2.0 and references cited therein; spatial 


distributions were obtained for a total of 25 contaminants and 64 PCB congeners. 


Spatial distribution for all contaminants found are plotted by Halfon (1989); data are also 
found in the reports of Eli-Eco (1987) and Zenon (1987). Average concentrations of the 
organochlorine contaminants and total PCBs are given in Table 3.21 for each of the three 
survey periods. The data shown in Table 3.21 indicate that few contaminants are present 
in relatively high concentrations (>1 ng/L). Such contaminants include the 
dichlorobenzenes (DCB), «-BHC, 7-BHC (lindane) and, in October, total PCBs in 23 of 30 
samples collected. This is most significant for total PCBs for which the PWQO is 1 ng/L. 
Levels of 1,2,4-trichlorobenzene (TCB), 1,2,3-TCB, and 1,2,3,4-tetrachloro-benzene (TeCB) 
were between 0.1 and 1 ng/L. Other contaminants averaged below 0.1 ng/L; however, as 
a large number of potentially toxic substances are present, they all are taken into account 
during the impact analysis. 


Halfon (1989) classified the 26 stations according to a non-para- metric ranking procedure 


(Halfon and Reggiani, 1986) based on a vector approach method in which pairs of samples 


are compared according to the relative concentrations of all organochlorines 


3 - 80 


TABLE 3.21: AVERAGE CONCENTRATIONS (ng/L) OF ULTRA-TRACE 
ORGANIC CONTAMINANTS 


Toronto Waterfront, 1987 


re 


Chlorobenzenes June August October 
Wisiel Cosi) 107 Sheet IBITORI Oi) Si) Dolult oo pion eee ee ee 
1,3-dichlorobenzene 17.2 0.46 0.25 
1,4-dichlorobenzene a 8} 9.4 0.0 
1,2-dichlorobenzene Lad 0.47 12 
1,3,5-trichlorobenzene 0.016 0.022 0.028 
1,2,4-trichlorobenzene 0.54 0.19 0.29 
1,2,3-trichlorobenzene 0.13 0.030 0.044 
1,2,3,5-tetrachlorobenzene 0.010 0.005 0.0 
1,2,4,5-tetrachlorobenzene 0.048 0.053 0.057 
1,2,3,4-tetrachlorobenzene 0.14 0.076 0.067 
pentachlorobenzene 0.031 0.036 0.034 
hexachlorobenzene 0.052 0.062 0.074 
Organochlorine pesticides 
Hexachloroethane 0.017 0.030 0.009 
Hexachlorobutadiene 0.006 0.010 0.0 
2,4,5-trichlorotoluene 0.006 0.001 0.003 
2,3,6-trichlorotoluene 0.0 0.002 0.0 
Pentachlorotoluene 0.031 0.019 0.017 
a-BHC 4.2 4.7 233.5 
y-BHC (Lindane) Let 1.0 0.58 
Octachlorostyrene 0.008 0.0 0.005 
y-chlordane 0.010 0.009 0.030 
pp '-DDE 0.030 0.021 0.035 
pp'-DDD 0.015 0.012 0.028 
pp'-DDT 0.0 0.0005 0.004 
Mirex 0.0 0.0004 0.021 
Total PCBs 0.58 0.82 Pile 


Dene en recor tiy ERSTE LETT a ee 


and total PCBs. All possible pairs are intercompared, inconsistencies noted, and similarly 
behaving samples are grouped together. The result of the ranking procedures can be 
displayed as a Hasse diagram in which the most contaminated stations are located at the 
top (e.g. Figure 3.23) and then as a contour map showing areas where most contaminants 
are found (Figures 3.24 - 3.26 for the three cruise periods). These three maps show the 
combined plumes of 25 contaminants, thereby representing the co-occurrence of 
contaminants in any given area. It should be noted that the numerical values on these maps 


are only relative in nature; thus, they do not translate to any quantitative risk for exposure. 


In June 1987 two impact zones were found: one on the western side of the waterfront near 
Humber Bay and one near the Lakeview STP (Figure 3.24; rankings >4). The eastern and 
offshore areas are much cleaner (rankings <2); as the water temperature was found to be 
lower in this area it probably represents dilution of contaminants by cleaner cooler 


hypolimnetic waters. 


In July-August 1987 (Figure 3.25), the contaminant distribution was similar except that the 
plume from the Main STP area was shifted towards the Eastern beaches. Impacted zones 
also occur in the east off Scarborough (Highland Creek STP) and at one station offshore. 


In October 1987 (Figure 3.26), the impacted areas were again similar. Off the Eastern 


beaches, the impact of local sources is reduced by dilution with hypolimnetic waters. 


Positive identification of local sources of these contaminants is difficult as the concentration 
of most contaminants are relatively low and uniform. An exception was 1,4-dichlorobenzene 
(Figure 3.27), which showed an obvious source in Humber Bay during the June and August 
cruises. Local sources were also evident for 1,2-dichlorobenzenes and total PCBs. Several 
other chlorinated benzenes, including hexachlorobenzene, showed a general increase at the 
stations closest to shore. While this suggests multiple sources in the Toronto area, definite 


sources were not pinpointed. 


Figure 3.23: Hasse diagram showing sample ranking of Toronto Waterfront sites. Station 
labels are defined by Halfon (1990). See Halfon (1990) for further details. 
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Figure 3.24: Idealized plumes which summarize the concentrations of 25 contaminants 
according to the results of the ranking analysis (Halfon, 1990), May-June 1987. 
Here stations on the same level in the Hasse diagram (Figure 3.23) are joined 
by a line. In May-June 1987, two polluted areas are seen, one located on the 
western zone of the Toronto Waterfront and one located near the Lakeview 
STP outfall. These stations are ranked at level 4 or higher. The offshore 
areas of the lake and the Scarborough Bluffs are ranked at level 2 or lower. 
In this area the water temperature is low and contaminants might have been 
diluted in cleaner, cooler hypolimnetic waters. 


May 26 - June 11, 1987 
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Figure 8: In May-June 1987, we note two polluted areas; one 
located on the western zone of the Toronto Waterfront and one 
located near the Lakeview STP station. 


Figure 3.25: Idealized plumes which summarize the concentrations of 25 contaminants 
according to the results of the ranking analysis, July-August 1987. The pattern 
of pollution is similar to that observed in May-June. The impact plume off 
the harbour has shifted toward the Eastern Beaches while the plume in 
Humber Bay is practically unchanged. Zones impacted by a combination of 
toxic contaminants exist in the east off Scarborough and at one station 
offshore. Overall, the nearshore of the waterfront is more contaminated than 
the open waters. 


July 27 - August 7, 1987 


Scarborough 


Figure 9: In July-August, the pattern of pollution is similar to 
that one observed in May-June. The. impact plume off the Toronto 
Harbour has shifted toward the eastern beaches while the plume in 
Humber Bay is practically unchanged. 


Figure 3.26: Idealized plumes which summarize the concentrations of 25 contaminants 
according to the results of the ranking analysis, October 1987. The same 
areas identified in the previous two cruises are again observed. The only 
notable difference is an area offshore with rank 9. Off the Eastern Beaches, 
hypolimnetic waters lower the impact by toxic contaminants from shore. 


October 5 - 9, 1987 


Scarborough 
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Figure 10: In October toxic pollutants impact the same ares 
identified in the June and August cruises. 


FIGURE 3.27 


Concentration of 1.4-dichlorobenzene (ng/L) as observed 
by ultra-trace methods, Toronto wdaterfront, June 1987 


MAY 26-JUNE 11.1987 _14-DCB 


As part of Lake Ontario, the Toronto Waterfront is exposed to whole-lake circulation 
patterns, which tend to replenish the water about every six to seven days. Thus, the 
baseline level of contaminants in the Toronto Waterfront is affected by remote sources as 
well as local sources. Only if all these sources were removed, would contaminant 
concentrations decrease rapidly to negligible levels, and rapidly only for water, suspended 
sediments and plankton. Levels in bottom sediments and the higher levels of the biotic 
food chain (e.g. fish) would decrease more slowly due to deposition (sediments) and 


bioaccumulation. These factors are being addressed in the fate modelling study. 


3.3.3 Relative Importance of Various Sources to Toronto 
Waterfront Water Quality 


Besides the Main STP, Lake Ontario in the Metro Toronto waterfront area is impacted by 
two other STPs within Metro Toronto: the Highland Creek STP and Humber STP, and six 
tributaries: the Etobicoke Creek, Mimico Creek, Humber River, Don River, Highland 
Creek and Rouge River. Numerous storm and combined sewers also discharge to the lake 
along the entire waterfront. As well, the atmosphere is another potential source of some 
contaminants such as PCBs. In addition, the Lakeview STP discharges to Lake Ontario in 
Mississauga, immediately west of the Metro Toronto limits. It is also considered here due 


to its potential for impacting upon the Toronto Waterfront. 


Loadings from the four Toronto-area STPs discharging to Lake Ontario were calculated 
from MOE concentration data collected in 1987 and flows supplied by the plants (Table 
3.22-3.23). A comparison of these loadings shows that the Main STP is by far the greatest 
contributor of conventionals (except for NO, + NO,-N) and most heavy metals to the lake. 
In fact, for many conventionals, as well as chromium, copper, mercury and lead, the loading 
from the Main STP is similar to or greater than the combined loadings from the other three 
STPs. Most of this difference is due to the higher effluent flow rate at the Main STP (2.1 
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x Humber STP, 3.3 x Lakeview STP, and 4.7 x Highland Creek STP). For heavy metals, 
Lakeview STP is the greatest contributor of nickel, while Lakeview and Main contribute 


similar amounts of manganese. 


Additional data provided to the MISA advisory committee (MOE, 1989) also indicates that 
in general, the Main STP is the highest contributor of metals of the nine STPs for which 
data were reported, although loadings for several metals from Ottawa (Green Creek) were 
somewhat higher. For most metals, mass balance data showed that the Main STP achieves 
above average removal efficiency; Highland Creek STP was also included in the study and 


showed above average removal efficiency for most metals. 


For bacteria (Table 3.24), on an annual or summer (June and July-August) basis, Lakeview 
STP tends to be the greatest contributor and Highland Creek the least (except for 
P.aeruginosa); however, the day-to-day and month-to-month variability at most locations is 
high. At Lakeview, both scales of temporal variability are less than at other locations, with 
July-August bacterial loadings only slightly less than June or October. In fact, the bacterial 
loadings from the Main STP in October (after cessation of chlorination) are higher than at 
the other STPs. 


Loadings of chlorobenzenes, organochlorine pesticides and PCBs as measured using 16-L 
samples (Sections 2.0, 3.1.1) are summarized in Table 3.25. Except for hexachloroethane, 
the Main STP is by far the greatest contributor of toxic organics, while the Humber STP 
is frequently next highest. While this may appear contradictory to the results of Section 
3.3.2, where the most obvious source of toxic organics was found to be the Humber Bay 
area, it should be re-emphasized that the closest ultra-trace sampling locations to the Main 
STP (upon which section 3.32 was based) were 3 km from the outfall (due to limited water 
depth for the CCIW vessel during the ultra-trace sampling). Thus, the contaminants from 
the Main STP were far more diluted before being sampled, relative to those from the 
Humber STP. At any rate, the high concentration of 1,4-DCB observed in Humber Bay 
(Figure 3.27) is consistent with its loading of 330 g/day (Table 3.25). In addition, restricted 
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TABLE 3.22: LOADINGS OF CONVENTIONALS 
ANNUAL MEANS, 1987 (10° kg/day) 


(a) Ammonia-N 


Main STP 20.0 6.5 8.9 34.6 47 
Humber STP ao 128 har LOSS 48 
Lakeview STP Zak J hep 0.06 yey 49 
Highland Creek STP 3 3 On 7 Oeli 4.7 48 
(b) Total Kjeldahl N 

Main STP 32.1 53) 1 2ie 380. 45 
Humber STP 7,0 Zao Bye) Slat 48 
Lakeview STP SH5 2142 0.6 7.4 49 
Highland Creek STP 4.5 0.8 2.6 672 48 
(c) Nitrite + Nitrate-N 

Main STP 0.86 219 0.03 110. 48 
Humber STP 0.15 0.37 0.01 25 47 
Lakeview STP 2.65 34, 0.16 6.6 49 
Highland Creek STP 4.47 0.84 2] G3 48 
(d) Total Phosphorus 

Main STP 0.79 1.69 0.10 10.89 47 
Humber STP 0.29 0.30 0.07 1.61 48 
Lakeview STP 0.37 0.22 0.12 30 49 
Highland Creek STP 0.29 Onze 0.07 e265 48 
(e) Suspended Solids 

Main STP 271 75.6 33 398. 47 
Humber STP 6.1 6.6 inet: 37 48 
Lakeview STP 4.2 3.8 jets 19 49 
Highland Creek STP 8.8 12249 13 47 47 

(contid) «tom, 


TABLE 3.22 (cont’d): LOADINGS OF CONVENTIONALS 
ANNUAL MEANS, 1987 (10° kg/day) 
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(f) Total Dissolved Solids 
(Estimated as Conductivity x 0.65) 


Main STP Sup. 142. 352. 961. 47 
Humber STP 294. 19. 23%.. 328. 47 
Lakeview STP 184. 20. 132. 218. 49 


Highland Creek STP 145. 34. 99. 211% 47 


(g) Total Filtered Residue 
(Measured Directly by Lab) 


Main STP 475. P22, 1a be ve 734. 46 
Humber STP 247. 26. 201. 334. 48 
Lakeview STP 178. 2s 30. 222% 49 
Highland Creek STP De 26. 78. 165. 47 
(h) Chemical Oxygen Demand (as 0) 

Main STP 87.2 81.0 33.9 442. 41 
Humber STP 24.9 15.2 =e 952 42 
Lakeview STP 24.3 12.0 10.5 64. 46 
Highland Creek STP 20.4 14.5 8.1 70 43 
(1) Dissolved Organic Carbon 

Main STP 13.4 3.4 6.1 pap Nea 47 
Humber STP 329 0.5 263 4.6 46 
Lakeview STP 4.9 2.5 Ze) 13.8 49 
Highland Creek STP 223 0.4 0.6 3.0 47 


TABLE 3.23: LOADINGS OF HEAVY METALS 
ANNUAL MEANS, 1987 (kg/day) 


(a) Chromium 


Main STP Sle 48. 6. 270. 48 
Humber STP 4.7 hed! 0. 36 48 
Lakeview STP Ve) thet 0. 31 50 
Highland Creek STP Sia 4.4 0. Zu 48 
(b) Copper 
Main STP 26. 46. OF 2 GE 48 
Humber STP 35 4.7 OF Zr 48 
Lakeview STP Sal Bio 0. 16. 50 
Highland Creek STP 9 TES 2 al, 50 48 
(c) Iron 
Main STP 12x 0%) 3298 10° On, xel 02m) 21x. 102. aS 
Humber STP 8-4. x6102 6: 7/ax 10? © O@2ex 102 840°x 104gemas 
Lakeview STP 4.7.x, 102 4 Sex 102 @ O84ex 10? 217 xP 10200 
Highland Creek STP 4.8 xf102 5,8fx 1072 €,0¢9 x10? 27 x 10@gm as 
(d) Mercury (g/day) 
Main STP 64. L22% 10. 653. 48 
Humber STP 72.9 PEE 4, on 48 
Lakeview STP 31. 33% 3 130. 48 
Highland Creek STP ihe Ooo 2 42 48 
(cont'd)... 


TABLE 3.23 (cont’d): LOADINGS OF HEAVY METALS 
ANNUAL MEANS, 1987 (kg/day) 


ee wee ee ee ee we ee ee ee ee es ws ws me me wm we we oe ew i ww ww ww ww ww wr wr wn ee wr ern eens esse snwsenrserss_ we eo 


(a) Manganese 


Main STP 74. 28. 36. 156. 48 
Humber STP 47. L3% 26. 79. 48 
Lakeview STP 78. 38. 20. 148. 50 
Highland Creek STP 61. 228 30) 119. 48 
(b) Nickel 

Main STP 34. r33 17. 65. 48 
Humber STP 24 10. 4. 41. 48 
Lakeview STP 60. B2e 19. L295 48 
Highland Creek STP Bh 2 23 125 48 
(c) Lead 

Main STP 12. 19 0. 93 48 
Humber STP 4.0 4.3 0. 14. 48 
Lakeview STP 3.0 2.4 0. 10. 50 
Highland Creek STP Do) 250 0. LE 48 
(d) Zinc 

Main STP 56. 110. 0. 550. 48 
Humber STP 40. 26. 17. 153. 48 
Lakeview STP tae hes bays, 126 83. 50 


Highland Creek STP 14. os 4a... 64. 48 
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TABLE 3.24: LOADINGS OF BACTERIA 
ANNUAL MEANS, 1987 (organisms/day) 


eee eee we ee a a es we we we en we wee we we we we we a we we we on ow a ww wa ww ww ww ww ow sess sess ess se 


Geometric Lower Upper 
Parameter/Location Mean Limit Limit Min. Max. N 


(a) Fecal Coliforms 


Main STP ox low 8.7x1012 9.5x1014 2.0x1034 1.7x1076 46 
Humber STP 3.2x1039 6.4x1014 1.5x10!5  <4.5x10!° = 5.6x1076 45 
Lakeview STP 1.1x10" 3.3x1024 3. 4x1042 7.4x10!2 8.5x101 46 
Highland Creek STP Peale alah 3.0x1022 1.5x102% 78x10) exo 46 
(b) E.coli 

Main STP 220x105. 5e3x10-8 Vogal 127x101) el 5x1 Oe 46 
Humber STP 25x 10t 4.0x101! x1 0c. 2.0x1029 4.3x1046 45 
Lakeview STP 7.8x1039 2.2x1023 2.8x1014 7.0x1022 7.3x1014 46 
Highland Creek STP 1.8x107° 2.431042 3x1 00° 5.6x1012 1.0x10 46 
(c) Fecal Streptococci 

Main STP Tees tL O23 7e.6x10.5 178x10*> 4.5x10!° 9.8x1014 46 
Humber STP Tea x10. 14x10. 1s2 1020 3.8x10!° 3.0x10}° 45 
Lakeview STP 1. .8x10+> 3.8x10}4 8.7x10}3 1.9x1022 4.0x10** 46 
Highland Creek STP 5 xis 1.6x10}4 1. 7x10 1. 7x1024 95. 2xtous 45 
(d) Pseudomonas aeruginosa 

Main STP Sop alee 4.3x102° Be oxlOs- 27 4x10'° weno Ox Oise 46 
Humber STP 1.4x10?4 7.1x10?° 28%L019 .<8.47% 102 Be2x1O.) 45 
Lakeview STP 1.4x1014 Sp alaee 4.3x1012 1.8x101! 1.9x10}5 46 
Highland Creek STP saa 1;8x104 3. 1x10 128x109" 27*0x1022 46 


TABLE 3.25: LOADINGS OF ULTRA-TRACE ORGANIC COMPOUNDS* 
ANNUAL MEANS, 1987 (g/day) 


oe ae a eee a ee ee a ee ee a a = a a a a a a a a a a a a a a a ae aerrerawr awrarrerrrnrnrnraarasanrasrnrsaaeeeeee 
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(a) 1,3-dichlorobenzene 


Main STP 78. 38. Za 15 Se 10 
Humber STP Zor 21% 6. 74. 9 
Lakeview STP 5.4 4.1 0. 123] 10 
Highland Creek STP 2.4 Zine 0 183 10 
(b) 1,4-dichlorobenzene 

Main STP 630. 260. 100. 990. 10 
Humber STP 330. 120. 150. 500. 9 
Lakeview STP 130. 47. 84. 230. 10 
Highland Creek STP 110. 270. 36. 140. 10 
(c) 1,2-dichlorobenzene 

Main STP 280. 130. 63. 480. 10 
Humber STP 95. a5 24. 180. 9 
Lakeview STP 44. Lae 16. 85. 10 
Highland Creek STP 8.3 4.4 249 18. 10 
(d) 1,3,5-trichlorobenzene 

Main STP 19,3) Zr _ 0.4 7.4 10 
Humber STP 0.34 0.12 OeL7 0.52 9 
Lakeview STP 0.13 0.07 0.05 0.23 10 
Highland Creek STP 0.05 0.05 0. 0.17 10 
(e) 1,2,4-trichlorobenzene 

Main STP 116. 118. 136 Sie OF 10 
Humber STP 6.3 1.6 4.4 9.8 9 
Lakeview STP 8.3 16.9 1.4 56. 10 
Highland Creek STP bie 0.03 0.6 via 10 
(£) 1,2,3-trichlorobenzene 

Main STP 54. 79. 6.1 210. 10 
Humber STP 0.96 0.25 0.50 1.44 9 
Lakeview STP 0.50 0.28 0.22 eS 10 
Highland Creek STP 0.20 0.07 0.09 0.28 10 


Note * - data based on large-volume (16L) samples. 


(cont'd)... 


TABLE 3.25 (cont'd): LOADINGS OF ULTRA-TRACE ORGANIC COMPOUNDS 
ANNUAL MEANS, 1987 (g/day) 


oo ee we ee i we i ew we a a a ww a ww a a a a ann eens ses sees 


(g) 1,2,3,5-tetrachlorobenzene 


Main STP OFZ 15.8 Om 42. 10 
Humber STP 0.07 0.05 OF 0.15 9 
Lakeview STP 0. 0. 0. 0. 10 
Highland Creek STP 0.04 0.03 0. 0.08 9 
(h) 1,2,4,5-tetrachlorobenzene 

Main STP 19. 15: 4.3 Ay. 10 
Humber STP 0.35 0.17 0.08 0.59 9 
Lakeview STP 0.08 0.04 0. 0.13 10 
Highland Creek STP 0.26 0.11 0.10 0.44 10 
(1) 1,2,3,4-tetrachlorobenzene 

Main STP 68. 101. ta Zils 10 
Humber STP 0.49 0.22 0. 0.68 9 
Lakeview STP 0.82 0.04 0. 0.13 10 
Highland Creek STP 0.28 OF) 0.06 0.57 10 
(j) Pentachlorobenzene 

Main STP 6.9 De 10. 16 10 
Humber STP Ona2 0.10 0.28 0.61 9 
Lakeview STP 0.11 0.09 0 0.32 10 
Highland Creek STP 0.35 0.20 0 0.69 10 
(k) Hexachlorobenzene 

Main STP (fe il 6.1 0. alge. 10 
Humber STP 0.18 0.04 0.09 One2 9 
Lakeview STP 0.15 0.11 0. 0.39 10 
Highland Creek STP 0.17 ORL? 0. 0.65 10 
(1) Hexachloroethane 

Main STP 0.16 0.33 0. lel 10 
Humber STP 1203 phos hy’ OF 3.0 9 
Lakeview STP 0.22 0.24 OF 0.60 10 
Highland Creek STP 0.16 0.21 0. 0.61 10 


(cont'd) 


TABLE 3.25 (cont’d): LOADINGS OF ULTRA-TRACE ORGANIC COMPOUNDS 
ANNUAL MEANS, 1987 (g/day) 


ba Sa SE SS SS eS ata a ae me maha ar 


Parameter/Location Mean Std. Dev. Min Max N 
(m) a-BHC 
Main STP ss} 0.8 0. 29 10 
Humber STP 1.0 0.3 0.6 1.4 9 
Lakeview STP 0.44 0.33 0. 0.90 10 
Highland Creek STP 0.24 0.17 0. 0.49 10 
(n) y-BHC 
Main STP 21: 17. 0. 50 10 
Humber STP 15. Lon 0. 36 9 
Lakeview STP 11. 10. 0. 28 10 
Highland Creek STP 3.0 2.8 0. 8.5 10 
(o) y-chlordane 
Main STP 0.34 0.18 0 0.67 10 
Humber STP 0.08 0.06 0 0.14 9 
Lakeview STP 0.15 0.17 0 0.49 10 
Highland Creek STP 0.04 0.07 0 0.24 10 
(p) a-chlordane 
Main STP 0.43 0.42 0. 1.39 10 
Humber STP 0.15 0.10 ON. 0.29 9 
Lakeview STP 0.11 0.15 0. 0.51 10 
Highland Creek STP 0.05 0.42 0. ‘Os is) 10 
(q) Total PCBs 
Main STP 16.8 10.3 4.7 40. 10 
Humber STP Shaul | 730) 2 Tet 9 
Lakeview STP 3.5 5.8 0.4 20 10 
Highland Creek STP 1.4 Les 0.2 4.3 10 


circulation in Humber Bay may lower the dilution rate of contaminants relative to the Main 
STP area. Loadings of several other chlorobenzenes from the Humber STP are also 
reflected in the observed concentrations at station 2072. 


3.3.4 Interparameter Correlations 


Interparameter correlations within the receiving water are shown in Table 3.26. The best 
general surrogate for metals appears to be conductivity, although suspended solids, DOC 
and ammonia were also correlated with the same metals (Cr, Cu, Fe, Mn, Ni, Zn). No 
surrogates were found for Cd, Hg or Pb. 


Of the organic parameters measured at conventional detection levels, only lindane (7-BHC) 
and a-BHC were detected in receiving water with sufficient frequency to develop 
meaningful correlations. DOC provides the best general predictor of these two 
contaminants (r = 0.89 and 0.48 for 7-BHC and a-BHC, respectively). However, 
conductivity, COD, ammonia and phosphorus were also correlated. 


Correlation coefficients between surrogate parameters and ultra-trace (APLE) organics were 
also computed within the receiving water. However, as the APLE samples were not 
collected at the same time as the conventionals (because the APLE samples were collected 
by the CCIW survey vessels), it was necessary to use seasonal (June, July-August or 
October) mean values for conventionals at each station. Therefore, these correlations must 
be interpreted with caution. A subset file containing data for all survey dates only at the 
stations from which APLE samples were collected was prepared, and the requisite monthly 
means calculated using the same parameters as used in the STP effluent calculations 
(Section 3.1.4(b)) plus turbidity. The same procedures as described in that section was used 
to transform all data and produce correlations. The resulting correlations are shown in 
Table 3.27 (a-c) for each of the three survey periods. 
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In June and August a large number of significant (P<0.05) correlations are found, while 
in October very few significant correlations were found. In June, the best overall surrogates 
appeared to be total Kjeldahl nitrogen and total phosphorus, each of which had a large 
number of highly significant (P<0.01) correlations. Scatterplots of typical correlations are 
shown in Figures 3.28-3.29. In all these figures, the highest TKN or total P concentration 
was found at station 2072 (Humber STP); this station also represented one of the highest 
concentrations for most parameters and an average value for PCB. The highest PCB 
concentrations during this survey period was found at the far west end of the study area and 
represent "outliers" on Figure 3.26; the source of these PCBs is unknown and appeared to 


originate from outside the study area. 


The June data (Table 3.27(a)) appear to suggest that pH and conductivity would be good 
surrogates for several parameters. However, examination of scatterplots indicated a very 
narrow range for these parameters at most stations, with the correlations entirely 


determined by high conductivity and low pH at the Humber STP station only. 


In August, total phosphorus again was a good surrogate for many of the organic parameters 
(Figure 3.29) but total Kjeldahl N was less successful. Conductivity, and suspended solids 
were useful predictors for some organics (Figure 3.29). High concentration of conventional 
and ultra-trace parameters were found at Stations 2072 (Humber STP), 1375 (inner harbour 
near Don River) and 2891 (3 km south of Main STP). The latter station especially is useful 


with respect to modelling and plume tracking. 


In October, very few successful correlations were found, and no conventional parameter was 
capable of predicting ultra-trace organic contaminant concentrations (Table 3.27c). 
However, a large portion of the contaminants exhibited very low to non-detectable 
concentrations at most stations in October, and consequently no reasonable predictions can 


be made for that survey month. 
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Halfon (1989) also correlated these organic parameters with temperature measured at the 
same time as sample collection. He found a positive correlation with many parameters in 
June, a negative correlation in July-August and no correlation in October. He was able to 
relate both correlations to sources at Humber Bay, which is similar to the results obtained 
here. His impacted zones also showed the effects of other sources including the Main STP. 
Both his work and that performed here suggest that similar transport mechanisms are 


important in determining the distribution of various contaminants. 


Overall, considering all data obtained, total phosphorus appeared to be the best surrogate 
of many ultra-trace organics. However, the correlation patterns are time-dependent and 
must be established separately for each particular case to be studied. While more 
successful than the results for the Main STP effluent, the usefulness of these results is 
limited by the fact that most samples were obtained in the least impacted areas of the 
waterfront (cluster 4; section 3.3.1). The lack of stations close to the Main STP outfall was 
a serious drawback; further work should emphasize this area and other impacted zones in 


order to obtain a better range of values for both conventional and toxic parameters. 


3.3.5 Temporal Trends 


Only a limited comparison of seasonal trends is possible with the given data, by examining 
the cruise means of various parameters within individual clusters and at selected stations 
(1419, 1536, 2029). } 


In general, conductivity, COD, ammonia, TKN and phosphorus tend to decline on the 
waterfront from June through August or September, with an increase again in the fall. 
Similar declines in conductivity and nitrogen parameters have previously been observed 
(Poulton and Griffiths, 1986); however, the latter study showed an irregular seasonal 
variation in total phosphorus in the inner and outer harbours. In contrast, turbidity and 


suspended solids tend to increase slightly from June through August, declining thereafter. 
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At the Main STP outfall (station 1419), deviations from the general trend were often 
observed, either by exaggeration (e.g. turbidity, suspended solids), or by opposition (e.g. 
phosphorus). Variability was often greater at Station 1419, commensurate with the greater 


concentrations here, and with the known variability in discharge rates from the Main STP. 


In comparison to 1986, free chlorine discharge was apparently lower in 1987 (annual means 
at Station 1419 were 1,018 + 1,143 and 270 + 271 ug/L, respectively; medians were 540 and 
324 ug/L, respectively), although the number of sampling days was small in 1986 (n = 4). 
Ammonia and phosphorus concentrations at the outfall were similar between years. Metal 
concentrations were slightly higher in 1987 for iron, chromium, mercury, manganese and 
nickel but lower for zinc (Appendix 2). Only the manganese and nickel increases were 
statistically significant (p<0.05). 


3.3.6 Comparison to Historical Data 


Table 3.28 shows historical water quality data from the Toronto Waterfront General Water 
Quality 1976-1983 Report (Poulton and Griffiths, 1986) and Boyd (1988) for parameters 
of interest. 


The spatial boundaries of the zones used in the 1986 report were calculated based upon the 
same parameters as in this report plus Secchi disk and Chla, which are both, on average, 
spatially invariant. The correspondence between the zones of this report and of Table 3.28 


are roughly as follows: 


3 - 107 


TABLE 3.28: COMPARISON OF CURRENT AND HISTORICAL DATA BY ZONE 


Parameter Zone 1 Zone 2 Zone 3 Zone 4 Zone _5 Zone 6 
1976 

NH3-N 0.08 0.13 0.26 0.56 
Total P 0.031 0.040 0.066 0.10 
Conductivity (umhos/cm) 350 360 430 510 
Turbidity (FTU) 2a 33.0 a yiite! 8.4 
1977-1978 

NH3-N 0.099 0.095 0.094 0.34 2.6 hts) 
Total P 0.02 0.03 0.04 0.06 0.12 0.16 
Conductivity (umhos/cm) 340 360 380 450 450 850 
Turbidity (FTU) 0.94 1.8 Zr: Jem, Jerk 14 
June 1983 

Total P 0.018 0.024 0.027 0.039 0.054 0.19 
Conductivity (umhos/cm) 350 351 370 420 370. 960 
Turbidity (FTU) Iv} yy? Ly 270 1.8 92.2 
1985 

Total P 0.01 0.03 0.03 0.20 0.30 0.3 
Conductivity 324 346 346 429 465 465 
Turbidity (FTU) 251 Ose: Et, bad 7.5 Vise 
1987 

NH3-N 0.05 0.06 0.06 0.56 5.6 

Total P 0.01 0.04 0.04 0.06 OeZ 
Conductivity (umhos/cm) 316 330 330 370 500 


Note: Units are mg/L unless otherwise stated. 


3 - 108 


Zones in Toronto Waterfront 


This Report 

Figure 3,2 Water Quality Report 
Strong Impact Zone 5 
Moderate Impact Zone 4 

Slight Impact | Zones 2 and 3 

No Impact Zone 1 


However, both the stations sampled and the season of sampling have varied from year to 


year. 


Available data suggest that conductivity has been generally decreasing over the past decade, 
but an opposite trend is seen in Zone 5 in the vicinity of the Main STP (Station 1419) 
from 1983 through 1987. Similarly, ammonia has been generally decreasing in 
concentration, but has been increasing in Zone 4 (Don River mouth) and Zone 5 from 
1978 through 1987. Turbidity was apparently decreasing from 1976 through 1983, but 


increased again in 1985. 
Total phosphorus has been generally increasing since 1983, but this trend is much stronger 


in Zones 4 and 5, suggesting that the Don River and STP are major contributors to this 


trend. 
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4.0 DERIVATION OF EFFLUENT LIMITS 


In this section, several different methodologies useful for setting effluent limits are 
discussed. The methods discussed may be useful on their own, or in conjunction with the 


various toxicity tests, the details of which are presented in other component reports. 
4.1 Th ical iction of Dilution D Dispersion 


McCorquodale (1983) described models for calculating the rate of effluent dilution from 
a diffuser outfall in the near field (initial mixing at end of pipe) and far field (dilution due 
to water currents and dispersion). Initial mixing at the end of pipe is dependent upon flow 
rate at each diffuser port (increased flow rate causes poorer dilution), water depth 
(improved mixing at greater depth), and relative plume buoyancy (expressed as a stability 
factor; more buoyant plume results in poorer dilution). Far field mixing is then treated 
using an unsteady state flow model involving solutions of the convective-dispersive equation 
for a volume source generated from the near field mixing. Further details are given by 
McCorquodale (1983) and references cited therein. 


Pending development of MOE’s mixing zone policy, the model was run for estimation of 
initial (near field) dilution (i.e. zero radius of mixing zone) as well as a 500 m radius of 
mixing zone. Dilutions were calculated for the existing outfall location (100 m diffuser 
discharging into a 6 m total water depth), and the proposed offshore location (1000 m 
diffuser discharging into a 15 m total water depth) (Gore and Storrie, 1986b). Three 
current speeds were used: (a) stagnation, initially considered as 0.02 m/s; (b) average 
conditions, 0.1 m/s; and (c) jetting, 0.2 m/s. Three Main STP discharge flow rates were 
also considered: (a) 11 m*/s, which represents the summer 1987 average flow rate; (b) 
15.77 m°/s, which represents the ultimate average design flow, and (c) 35.61 m?’/s, which 
represents the maximum design flow for the input sewer system (Gore and Storrie, 1986b). 
The plume stability factor (f, McCorquodale, 1983) was taken as 0.417 (mixing in top 2.5 
m of water column) with the existing outfall. This produces initial dilutions under average 


4s a 


conditions near to those observed (Section 3.2.2) and is in reasonable agreement with 
profiles taken near the end of the outfall (Section 3.2.5). 


The resulting dilutions are given in Table 4.1. Dilutions were not computed for the 
maximum design flow (35.61 m°*/s) at the existing outfall, as this would represent an 
"overload" condition, which would result in partial discharge through the seawall gates. 
For same cases, the "stagnation" current speed of 0.02 m/s had to be increased as large flow 


rates into a confined space would artificially create faster currents in the near field. 


The dilution factors given in Table 4.1 are expressed as a dilution ratio (X value in Section 
3.2.2, where background value is assumed to be zero) and also as a percentage of source 
concentration. The dilution ratio can be interpreted as 1 part of effluent diluted to the 
stated number of parts of effluent - water mixture, i.e. 1:2.3 dilution for average currents 
and average STP outfall flow equals 1 part effluent plus 1.3 parts of lake water giving 2.3 
parts of mixture. This dilution ratio is slightly lower (i.e. poorer dilution) than the 1987 
observed results at station 1419 (Table 3.16) for conventional parameters. (Range 1:2.9 for 
lab ammonia to 1:7.1 for zinc). The observed results also indicate somewhat more far field 
dilution with distance (i.e. 3 to 15% at 500 m); however, as decay was neglected in these 
dilution calculations, the resulting estimates are conservative. It should be noted that 
parameters with significant decay factors such as fecal coliforms, total phosphorus and 
metals had the highest initial dilution factors in table 3.16; those with less decay generally 
had initial dilutions below 1:4. Nevertheless, it was felt that these should represent a 
realistic range of expected dilutions should the new outfall be built. 


A selected subset of these dilutions are being used in conjunction with laboratory and field 
toxicity tests, to estimate acceptable effluent strengths based on the precept of no rapid 
lethality in the mixing zone and no chronic effects beyond the edge of the mixing zone 
(Policy 5, MOE, 1984). 


TABLE 4.1: RESULTS OF DILUTION MODEL RUNS FOR CURRENT 
AND PROPOSED MAIN STP OUTFALLS 


Main STP Current Mixing Depth Initial Dilution Dilution 
Flow Speed From Surface Dilution at at 
(m3/s) m/s m (no m.z.) 250 m 500 m 


(a) Existing Outfall (100 m diffuser length in 6 m total water depth) 


11 0.05 2.5 1.1 (88%) 2.1 (47%) 3.0 (338%) 
0.1 Lak 2.3 (44%) 3.2 (32%) 4.3 (24%) 
0.2 205 4.5 (22%) 5.6 (18%) 6.2 (16%) 

15.77 0.06 3% 1.1 (88%) 1.3 (76%) 1.7 (588%) 
0.1 5: 1.9 (53%) 2.4 (42%) 3.1 (32%) 
0.2 3. 3.8 (26%) 4.3 (23%) 4.6 (21%) 

(b) Proposed Outfall (1000 m diffuser length in 15 m total water depth) 

Jl up len ey oH 5.5 (18.3%) 15.7 (6.4%) 24.7 (4.18%) 
0.1 SNe 27.3 ( 3.7%) 33.8 (3.0%) 42.9 (2.3%) 
0.2 af 54.5 ( 1.8%) 61.7 (1.6%) 66.4 (1.5%) 

15.77 0.02 3. 3.8 (26.3%) 10.9 (9.1%) 17.2 (5.8%) 
0.1 Sh 19.0 ( 5.3%) 23.5 (4.3%) 29.9 (3.48%) 
0.2 Se 38.0 ( 2.6%) 42.9 (2.3%) 46.3 (2.2%) 

35.6 0.03 She 2.5 (408%) 5.0 (20%) 7.0 (14.3%) 
0.1 o 8.4 (11.9%) 10.4 (9.6%) 13.2 (7.5%) 
0.2 JF 16.8 ( 5.9%) 19.0 (5.3%) 20.5 (4.9%) 


Note: Current speeds of 0.02-0.06 m/s represent “stagnation” speeds of 
0.1 m/s, "average", and 0.2 m/s, "jetting". 


4.2 Use of Theoretical Dilution Model to Derive Effluent Limits 
B n r lity Resul 


The theoretical dilution factors given above may be used to estimate the extent of loading 
reduction required based upon the water quality results in the effluent (section 3.1) and 
Lake Ontario in the vicinity of the Main STP outfall. 


Derivation of water-quality based effluent limits is based upon the premise of achieving 
PWOO at the edge of a defined mixing zone, e.g. zero, 250 or 500 m radius from the outfall 
diffuser. For this purpose, parameters which are most frequently in exceedence are 
discussed. Table 3.17 shows that residual chlorine, ammonia and total phosphorus are 
almost always in exceedence at station 1419 near the end of pipe; of the heavy metals 
surveyed, iron, copper and cadmium are frequently in exceedence. Of the organic 
parameters sampled, lindane (Table 3.18) is frequently in exceedence. In addition, across 
the Toronto Waterfront, PCBs frequently exceed the PWQO of 1 ng/L as determined by 
APLE sampling (section 3.3.2). Within several kilometres of the Main STP (stations 1536, 
2886, 2892 and 2029), the average concentration of total PCBs was found to be 0.3, 0.5 and 
1.0 ng/L during the three 1987 surveys (Halfon, 1990). These figures are lower than the 
values given in Table 3.21 because the latter include results for several points impacted by 
other sources. The average of the above values (0.6 ng/L) has been used as a "background" 
value in this chapter; however, considering the range of values observed, "background" could 
be 1.0 ng/L or higher, in which case required effluent reductions would be 100% for all 
scenarios. This fact should be kept in mind in evaluating the proposed effluent reductions 
for PCBs. However, it can be anticipated that abatement measures taken to reduce the 
loadings of other substances of concern should also reduce PCB loadings and have a 
beneficial effect on concentrations in various components of the ecosystem. 


The Toronto MISA team has recommended that water quality based effluent limits be 
derived on the basis of 95% compliance with objectives or guidelines. For this basis, the 
probability distributions of loadings for parameters of concern (section 3.1.2) are used. For 


substances such as total phosphorus in which the observed background concentration 
(cluster 4; section 3.3.1) is a significant fraction of the objective or guideline, the allowable 
increase in concentration near the outfall is used. In addition, as the dilution characteristics 
of the proposed outfall location are much better than the present location (Table 4.1), 


effluent limits are derived separately for each location. These will be discussed in turn. 


4.2.1 P n fall ion 


Calculations have been performed for mixing zones of zero, 250 and 500 m, assuming 
average currents (0.1 m/s) and average Main STP flow (960 x 10° m°/d = 11 m’/s) 
(dilution factors of 2.3, 3.2 and 4.3, Table 4.1). As we are interested in the maximum 
allowable increase over background, the allowable effluent concentration (c) calculated as 
follows (500 m): 


c = (4.3 x (PWQO - BGD)) + BGD 


The allowable concentration is then multiplied by the corresponding flow to obtain the 
allowable loadings, which is compared with observed 1987 95 percentile loadings (Figures 
3.10 - 3.15) to derive the required loading reductions. The results are given in Table 4.2. 
Similar calculations were done for the ultimate average design flow (15.77 m’/s); as these 
give similar loading reductions, they are not reproduced here. 


Adjacent to the outfall (zero mixing radius), severe reductions in effluent levels are 
predicted, ranging from 68% for lindane to 98% for total phosphorus and 99.7% for 
chlorine residual. If a mixing zone of 500 m radius is permitted, the required reductions 
range from 42% for iron and lindane, to 96% for total phosphorus and 99.4% for chlorine 


residual. Intermediate results are obtained for a 250 m radius mixing zone. 


The model was also used to estimate the mixing zone length required for dilution of the 
total phosphorus concentration to the 20 ug/L guideline for the prevention of nuisance 


algae growth, assuming effluent concentrations of 0.5 and 1.0 mg/L and a background 
concentration of 10 ug/L. (i.e. 1:49 and 1:99 total dilutions, respectively). The results, for 
average currents and average Main STP flow, were 17 and 35 km, respectively. At these 
high distances, the model predictions may only be considered as crude approximations, as 
the model neglects small-scale variability in currents and dispersion, which strongly affect 


the dilutions. 


4.2.2 P utfall ation 


Effluent requirements for the proposed location are also given in Table 4.2. Due to the 
much greater extent of dilution predicted as a combined result of increased water depth and 
diffuser length, the present loadings are predicted to satisfy requirements for all parameters 
but total phosphorus, chlorine residual and PCBs at all distances, and ammonia at zero 
mixing radius only. Required reductions for chlorine residual range from 94 to 97%, and 
for total phosphorus, 65 to 80%, and for PCBs, 58 to 75%, dependent on mixing zone 


radius. 


As with the present outfall above, the model was used to estimate the mixing zone length 
required for dilution of total phosphorus to the PWQG of 20 ug/L from effluent 
concentrations of 0.5 and 1.0 mg/L, assuming a background concentration of 10 ug/L. The 
result indicates that distances of 0.7 and 3.6 km are required for dilution to 20 ng/L under 


these conditions, respectively. 


4.3. Dilution Envelope Modelling 


Drogue-tracking involves the tracking of multiple water-sail drogues (current tracers) 
released simultaneously at a point source. The positions of the drogues then give a 
measure of the direction and the spread of a plume as a function of the time. Current 
direction information from the 1986 drogue cluster tracking study was used in determining 
down-current dilution envelopes (pages 4.8 - 4.9); however, these data were insufficient to 
develop statistical dilution-time relationships due to the frequent presence of non-steady 
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state and eddy conditions. Consequently, information from other STP studies was also used 


to develop these relationships. 


Between 1984 and 1986, a statistical dilution and time relationship for steady, conservative 


sources, was developed for nearshore surface waters (Figure 4.1; Palmer, 1988): 


(Cm 1555 tea (4.2.1) 
C. 


Since the conductivity of the water is related to the amount of solids dissolved in the water, 


it is possible to develop dilution envelopes from Equation 4.2.1, based upon the known 


conductivity: 
Ce ul i ed net (4.2.2) 
C, Ko - Kyae 


where: K, is the value of the surface conductivity at the STP 
outfall, and 


K 


nae is the background conductivity. 


Using the conductivity data collected in 1987, we can estimate the dilution using Equation 
4.2.2 and obtain the value of the time from Equation 4.2.1. Knowing the distance to a 
particular station then defines the velocity. The average velocity was estimated for stations 
1419 to 2891, which lie in a southerly direction from the outfall, and stations 1419 to 2209, 
which lie in an easterly direction from the outfall, and dilution envelopes were developed 
(see Figure 4.2). From the statistical dilution/time relationship of Figure 4.1, we can 
estimate the error involved since, for a particular value of the dilution, the predicted time 
t, ranges betweent,., andt,,,. The error calculated is shown in Figure 4.2 for a 10% 


dilution. 
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FIGURE 4.2 
Conductivity Dilution Isopleths at 0.5 m in 1987 


YU ff} Confidence band on 
Yyy Yy 0.1 isopleth 


An observed concentration isopleth for phosphorus calculated directly from the Main STP 
study data (see Section 3.2.2) and the predicted concentration isopleth for a strictly 
conservative source, are compared in Figure 4.3 for a 0.05 mg/L concentration, along with 
the error for each isopleth. The more rapid observed dilution in Figure 4.3 represents the 


non-conservative nature of phosphorus due to bio-uptake, sedimentation, etc. 


4.4 Effluent Dilution from Plume Tracking Results 


Derivation of water quality based effluent limits is based upon the premise of achieving 
PWOO at the edge of a defined mixing zone, e.g. 500 m radius from the outfall diffuser. 
Tables 3.17 and 3.18 show which PWQOs are exceeded at the outfall (Station 1419). At 
a water depth of 0.5 m, the PWQOs are exceeded for free chlorine, ammonia and 
phosphorus on all sampling days, and for iron, copper and cadmium 40%, 30% and 20% 
of the sampling days. The initial dilution from jetting from the outfall for residual chlorine 
was 1:3.0, for ammonia 1:3.1, for phosphorus 1:5.0, for iron 1:4.1 and copper 1:5.4 (Table 
3.16). The initial dilutions for the outfall diffuser should be the same for all water quality 
parameters. The variations determined for the different parameters are a consequence of 
the different sampling times for the effluent, receiving water and estimating the background 
concentrations in the lake. If the dilutions for the preceding parameters are averaged, the 
mean dilution achieved by the diffuser outfall is 1:4.1. This is somewhat higher than the 
value of 2.3 predicted for average currents (0.1 m/s) and average effluent flow from the 
dilution model discussed above (Table 4.1); it is approximately equal to the predicted 


dilution at 500 m under these conditions. 


Averaged dilution envelopes (Figures 4.1 to 4.3) were calculated from 1987 concentration 
data at stations close to the outfall. The mean velocity in the southerly direction was found 
to be 0.12 m/s (S.D. = 0.05 m/s) and in an easterly direction 0.06 m/s (S.D. = 0.03 m/s). 
The size of these dilution envelopes is directly proportional to the receiving water velocity 
- ie. doubling the velocity would result in a doubling of the size of the non-compliance 
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FIGURE 4.3 
Predicted versus Observed Phosphorus Concentration 
Isopleth at 0.5 m in 1987 
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zone. This fact is also illustrated in the size of non-compliance zones obtained from 
one-day plume tracking experiments in June 1987 (section 3.2.3; Figures 3.36 - 3.38 and 
Appendix 5). 


The rate of far-field dilution for conductivity, ammonia and chlorine residual may be 
estimated by plotting the dilution factors (X or 1/X) against distance from the outfall on 
a semilogarithmic scale, assuming first order decay with distance. Under this assumption, 
the rate constant of decay K (base 10) is simply the negative of the slope of the decay 
curve. The dilution rate, to 100 m (z) may then be expressed as: 


1.0 d 100 
| = = | (-k dx) (4.3.1) 
1.0-z 0 


where c = concentration 
x = distance from outfall. 


Assuming base 10, equation 4.3.1 is simplified to 


log 1.0 - log (1.0 - z) = -100k 


100k 


and z=1-10 (4.3.2) 


The percent dilution is then simply 100z%. 


The Y-intercept of the same plot represents an estimate of initial dilution, as it is the value 


of the dilution factor at zero distance. 


Using this methodology, dilution rates have been calculated for the plume tracking surveys 
undertaken in June 1987 (Figures 3.13-3.15 and Appendix 5). Only down-current 
information in the direction of maximum concentration was used; in addition, information 


obtained at distances at which the effects of eddies and non-steady state conditions occurred 


4 i= -13 


was neglected. Dilution curves for conductivity, ammonia and chlorine residual obtained 
using these assumptions are given in Figures 4.4 - 4.6. These curves emphasize information 
obtained within about 2 km of the outfall, as it is within this zone that we are most 
interested in dilution calculations; furthermore, data collected at greater distances are more 


subject to the effects of eddies and non-steady state conditions. 


Dilution rates and initial dilutions as calculated above are given in Table 4.3 for June 1987 
plume tracking data. Decay rates were not calculated for June 18 (stagnation) due to lack 
of ammonia and chlorine residual data at intermediate points in the current direction. 
Under jetting conditions (June 17 and 19), non-linear decay tended to be present; the decay 
rates calculated above represent the average decay at distances between 0.5 and 2.5 km 
from the outfall, and neglect "pockets" of higher concentration sometimes seen at greater 
distances (Figures 3.14b, 3.15b). 


Within the distances considered above, physical transport appears to be the major 
determinant of the dilution rate, with the highest dilutions being found under stagnation 
(June 18) and the lowest under jetting (June 17, 19). Within this distance and neglecting 
June 18, the dilution rates and initial dilutions for conductivity and ammonia were similar. 
This suggests that processes such as nitrification, bio-uptake and volatilization play only a 
minor role in ammonia decay over this time period. Furthermore, initial dilutions (1:2.9 
to 1:4.8) are similar to or slightly less than the average (1:4.1) predicted from water 
chemistry data. By contrast, chlorine residual has a lower initial dilution but a higher 
dilution rate. While decay is responsible for the higher dilution rate, the reason for the 


lower initial mixing is unknown. 
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Log dilution 


Figure 4.4 


Plot of dilution vs. distance from Main STP outfall 


Conductivity, June 1987 


Distance from Main STP outfall (km) 


June 16 
June 17 
June 18 
June 19 


Log dilution 


Figure 4.5 


Plot of dilution vs. distance from Main STP outfall 
Ammonia, June 1987 
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Figure 4.6 


Plot of dilution vs. distance from Main STP outfall 
Chlorine residual, June 1987 


Legend 
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A limited amount of additional data are available for 1986. Receiving water data are 
available at station 1419 near the outfall, plus stations 1431 to 1434. Drogue data are also 
available to define the downstream direction for dilution calculation; it was found that the 
conductivity at discrete points rather than concentration contours were used in the 
regressions. Station 1419 was also used due to the small number of data points available, 
despite the fact that concentrations at this station are affected by turbulent near-field 
effects. Using only down-current data, the percentage dilution rates per 100 m are given 
in Table 4.4. It can be seen that for ammonia and chlorine residual, the ranges and 
average values of both initial and far-field dilution are similar in both survey years. Data 
for fecal coliforms are also included. The high but variable initial dilutions in August and | 
September reflect the effect of chlorination; when this is absent in November, very low 


initial dilutions but somewhat greater far-field dilutions were found for this parameter. 
4.5 Effluent Limits Based on Plume Tracking 


The plume tracking results presented above may be used in conjunction with the 1987 
effluent data (section 3.1) to derive another estimate for effluent limits. Background data 
are also used in a manner identical to that described in section 4.2. As these calculations 
are based on observed data for the existing outfall, it is not possible to estimate effluent 
limits for the proposed outfall using this method. As it is desirable to deal with average 
to worst case conditions, the "stagnation" results of June 18 are neglected and the following 
calculations are based upon the average of the other three survey days (average to jetting 
current conditions). Conductivity is used as a surrogate for parameters not measured in the 
plume, as it shows significant correlation between the Main STP and area water quality 
stations (Table 3.19). 
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Given the average dilution rates (z) per 100 m for 1987 and average initial dilution at the 
outfall (ID), the average allowable effluent concentration to achieve compliance within a 


specific radius of the outfall may be calculated as follows: 


(a) Concentration at zero radius: 


CG) = (PWQO - BGD) * ID 43.3 
(b) Concentration at radius x in units of 100 m 
colmace Gaon 43.4 


i.e. at 500 m for a dilution of 9%/100 m, and allowable concentration is 1.60 C, = 5.8 
(PWQO-BGD), assuming an average initial dilution of 1:3.6. It should be noted that this 
allowable concentration is the excess over background, therefore the background value is 
added to the above concentration to give the allowable value, i.e. 6.6 (PWQO-BGD) + 
BGD. 


If the goal is that of achieving compliance 95% of the time, then the loading calculated 
from the above concentration should be compared with the 95th percentile of the observed 
loading in order to estimate the required loading reductions. The results are given in Table 
4.5. 


Calculated loading reductions at 500 m range from 22% for iron and lindane, to 98.9% for 
chlorine residual. The stringent loading reductions required reflect the limited amount of 
mixing at the present outfall site, and would be only slightly less if a larger mixing zone (e.g. 
1000 m) were allowed. 


4.6 Fate and Transport Modelling 
A simulation study of the fate and transport of toxic organic contaminants in the Toronto 


Waterfront was conducted by Halfon (1990) using ultra-trace organic contaminant data 
collected in 1987 as part of this study (sections 3.1.1, 3.3.2, 3.3.3) along with external 
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loadings (Niagara River, atmospheric) and currents as simulated by the Rand model (Gore 
and Storrie, 1989). 


The TOXFATE model used in this study simulates the concentrations of organic 
contaminants in water, sediments and several levels of biota. It accounts for processes such 
as volatilization, sorption, bio-uptake and sedimentation, as well as the hydrodynamic 
transport of contaminants input to the lake from Toronto-area sources (STPs and 


tributaries) and distant sources (Niagara River and the atmosphere). 


The objectives of this study included prediction of contaminant concentrations in the area 
surrounding the Toronto Main STP as well as the entire waterfront, estimation of the 
relative contribution of various sources to contaminant concentration, and estimation of 
the extent of loading reductions required to limit the zone of non-compliance with PWQO 
to within 500 m of the outfall. The model scenario that considered the zone near the Main 
STP outfall was used to compare with the estimates of required loading reductions from 


above. 


A load of 100 g/day of a hypothetical contaminant was discharged into the lake at the Main 
STP outfall, and its concentration in the grid square immediately adjacent to the outfall, 
plus those in the next nearest grid squares (i.e. 500 m from the outfall) were followed for 
a 12-day modelling period. As the residence time of a contaminant in such a small zone 
is very small, biochemical decay was neglected and physical transport (advection and 
dispersion) was the only dilution mechanism. Every five minutes, the concentration at the 
outfall grid square and the maximum concentration at the adjacent grid squares (i.e. worst 
case) were saved as two separate time series. These were analyzed by frequency and 
probability methods. A similar run was done for each of the three 1987 modelling periods, 
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both for the present outfall and for the proposed offshore outfall location (Gore and 
Storrie, 1986b). 


Figure 4.7 shows the frequency distribution of concentrations at the present location when 
all thee modelling periods are combined, and Figure 4.8 shows similar information for the 
proposed outfall. It is observed that with a loading of 100 g/day, the most likely 
concentrations at the present outfall are between 1 and 2 ng/L; very little additional 
dilution is observed at 500 m. With the proposed location, however, Figure 4.8 shows most 
likely concentrations of 0.1 to 0.2 ng/L at both distances - i.e. an improvement by a factor 
of 10. Some shifts in maximum concentration and shape of the distribution occur when 
data from individual months are examined; the reader is referred to Halfon (1990) for 


details. 


These results may also be plotted as cumulative frequency distributions - i.e. the probability 
of observing concentration below any value vs. that concentration value (Figure 4.9). With 
a loading of 100 g/day at the present outfall, there is a 95% probability that the 
concentration will be less than 11.5 ng/L at 0 m, and 8 ng/L at 500 m. With the proposed 
outfall, the 95% probability values drop to 3.3 ng/L at 0 m and 1 ng/L at 500 m. 


Both sets of results may be placed in the perspective of an actual contaminant simply by 
multiplying by the relative load, since only physical dilution and no_ background 
contamination are assumed. For example, the observed 1987 PCB loading from the Main 
STP is 16.8 g/day (Table 3.2 (d)). The 95% probability concentrations for the various cases 
are found by multiplying the values above by (16.8/100), i.e. for the present outfall they are 
1.9 and 1.3 ng/L at 0 and 500 m, and for the proposed outfall they are 0.55 and 0.17 ng/L. 
This suggests a further loading reduction by a factor of 1.9 is required to meet PWQO (1 
ng/L) at the outfall in the present case; no loading reduction is required in the proposed 
case. However, it must be remembered that this is assuming no background concentration, 
an assumption which is not accurate for PCBs (see Halfon 1990). While reduction of 
loadings from the Main STP will be most beneficial for some compounds like chloroben- 
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Probability of observing a contaminant concentration 
with a loading of 100 grams per day at the proposed 
offshore Toronto Main STP outfall. 
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zenes which have significant local input, the cleanup of the Toronto Waterfront requires a 
whole-lake approach. Contaminants from other sources, such as from the Niagara River, 
must be reduced to obtain a significant reduction in toxic contaminant concentrations in the 
Toronto Waterfront. 


As conservative behavior was assumed over a distance of 500 m from the outfall, the 
scenario described above is appropriate for deriving effluent limits for parameters currently 
in exceedence of PWQO near the outfall. It is assumed that the predicted concentration 
at zero and 500 m from the outfall is directly proportional to the loading applied. Using 
this assumption and the same (objective - background) values given in Table 4.2, the model 
is used inversely to estimate an allowable load; i.e. if it is desired to achieve a concentration 
of 8.0 ng/L 95% of the time 500 m from the outfall, then the maximum allowable loading 
is 100 kg/d (Figure 4.9). As before, these allowable loadings are compared to the observed 
95 percentile loadings to to derive effluent limits. The results are given in Table 4.6. 


At 500 m from the present outfall, this model predicts that the present loadings should 
satisfy effluent requirements for iron, copper, PCB and lindane. Observed data for this 
distance in 1987 indicate that all samples meet PWQO for iron and copper with the 
exception of one iron sample; thus the predictions are consistent with observations. 
Lindane was not measured at this distance in either survey year. No PCBs were measured 
using ultra- trace methodology at this distance. Required reductions for other parameters 
range from 24% for cadmium to 90% for total phosphorus and 98% for chlorine residual. 


At 500 m from the proposed outfall, all parameters are predicted to be in compliance, 


except for total phosphorus (17% reduction required), PCBs (50% reduction required) and 
chlorine residual (86% reduction required). 
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However, it must be kept in mind that the PCB dilutions were calculated with an average 
background concentration of 0.6 ng/L, which at times may be as high as 1 ng/L (i.e., the 
PWQO). Therefore, the required loading reduction at the outfall location may be as high 
as 100%; furthermore, as Halfon (1990) has stated, contaminant loadings from other sources 
such as Niagara River and atmospheric must be reduced in order to achieve the full benefit 


from local abatement measures. 


4.7 Consensus 


The allowable loads calculated in the previous sections are the load for which the defined 
objective (or guideline) will be satisfied for a specified fraction of the time, given the 
statistical variability of the effluent quality. In performing these calculations, recommended 
allowable loads were based on a compliance target of 95%. In other words, receiving water 
quality at the edge of the mixing zone would only be allowed to be out of compliance 5% 


of the time. A more stringent compliance target requires a more restrictive allowed load. 


In the previous sections, several methodologies were used in deriving effluent limits. Plume 
tracking results combined with inferences based on interparameter correlations allow a 
semi-empirical calculation of allowable loads. Modelling techniques including a simple 
dilution model as well as a more complex fate and transport model provide theoretical 
calculations which also allow for an estimate of allowable loads for the proposed offshore 
outfall location. 


Tables 4.7 and 4.8 present summaries of load allocations for zero and 500 m radius mixing 
zones respectively, as calculated by the methods described in the previous sections. These 
are based on a target of 95% compliance with PWQO at the stated distances from the 
outfall diffuser. 


It is observed that the simple dilution model predicts more stringent limits than the fate 


model; the plume track results are intermediate. As discussed in Section 4.4, the plume 
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tracking results are subject to inherent variability as a consequence of physical conditions 
which change from one day to the next - e.g. from stagnation to average currents to jetting, 
as well as rapidspatial variations in concentration in the initial mixing zone within tens of 


metres from the diffuser. Data in the table represent average to jetting currents. 


Predictions based on the dilution model are also sensitive to the choice of physical 
parameters such as current speed, mixing depth and total water depth (Table 4.1). The 
effluent requirements calculated in Tables 4.2, 4.7 and 4.8. represent average currents (0.1 
m/s); the assumption of jetting would result in higher initial dilution but poorer far-field 
dilution. Observed initial dilutions (Tables 3.16, 4.3 and 4.4) do not support this premise 
of higher initial dilution under jetting conditions; an improved model of initial dilution is 
to be developed in the near future (B. Kohli, pers. comm.). Therefore, the estimates of 
initial dilution at the proposed outfall location can only be considered as a first 
approximation to the actual dilutions which might be achieved at this location. 


The more complex fate and transport model predicts less stringent reductions than that 
predicted by either of the other two methods. This is particularly true for the heavy metals 
and lindane. The major limitation of this model is the size of grid square (S00 m) used; 
therefore, the "zero" mixing zone results actually represent an average concentration 
between 0 and 250 m radius of the outfall. Consequently, the allowable loadings in Tables 
4.6 and 4.7 more accurately represent predictions for a 250 m radius of mixing zone. 
Results for 500 m similarly represent an average of conditions between 250 and 750 m 
radius; however it should be noted that observed data in 1987 at 500 m (stations 2877 and 
2878) showed compliance with PWQO for all heavy metal data except one iron sample; 


lindane and PCBs were not measured at this distance. 


Predictions of dilution and consequent allowable loadings at the proposed outfall from the 
fate model, appear to be unrealistically high compared to the dilution model. However, 
predictions based on this model are strongly dependent upon the RAND model currents 
used as input; it is known that these predictions may not be as reliable at locations close 
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to the model boundary such as this one. These predictions therefore, do not appear to be 
defensible without further validation. It is suggested that further sampling for the 
parameters in question be undertaken, in order to validate the model predictions. This 
should include simultaneous sampling of the effluent, and in-lake at distances of 0 and 500 
m from the outfall in the direction of current movement. Sampling should be done on an 
intensive basis, for example, daily for one week in each month to obtain a statistically 
significant number of data points. See also Section 6.2 for further sampling 
recommendations. 
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5.0 CRITIQUE OF ASSESSMENT PROCEDURES 


In this report, we have illustrated various methods of deriving dilution contours and 
compliance zones from water quality data. Compliance (or non-compliance) zones, and 
effluent limits corresponding to specific compliance zones, are based directly on a dilution 


envelope model. Three types of dilution envelope modelling have been presented: 
ik, A simple graphical approach (Sections 3.2.2, 3.2.3 and 4.4); 


2. A theoretical approach involving calculation of initial mixing at the end of pipe 
from flow rate, water depth and relative plume stability, combined with far-field 


dilution from an unsteady state flow model (section 4.1); and 


3. A combination of an unsteady state hydrodynamic model driven by current meter 


and wind data with a pollutant transport model (Section 4.5). 


In support of the above, an extensive data base has been gathered, involving effluent quality 
at the Main STP, as well as Lake Ontario water quality in the vicinity of the Main STP and 
across the entire Toronto Waterfront. 


A fourth approach, not illustrated here, would make use of concentration and current 
velocity data in all directions, based on drogue cluster releases. With this information, a 
more sophisticated dilution envelope model can be constructed, which changes with current 
conditions. The information gathered in 1986 (Section 4.4; Table 4.4) was useful as an 
adjunct to the plume tracking data gathered in 1987 but was inadequate for full utilization 
of this approach. 


Some of the problems inherent in the modelling techniques have already been discussed 
in Section 4.6. These models are only as good as the data base which supports them, and 
an adequate description of the effluent and receiving water quality at various scales of 
temporal and spatial resolution is vital to the success of the project. 


The techniques used in this study which lead to the effluent limit setting approaches 


described in Section 4 have been summarized in Table 5.1. A few comments follow: 


As the effluent loading shows considerable temporal variability both at short- and long-time 
scales, an integrated approach to address this variability and its impact on the nearshore 
zone is important in establishing initial and far-field dilution rates. This would be best 
achieved by simultaneous time-sequential sampling of both the effluent and the lake 
adjacent to the outfall (station 1419), combined with plume tracking and tracking of drogue 
clusters. Several days of application of this technique in different seasons would help to 
cover seasonal changes combined with variability in the current climate (e.g. stagnation, 
jetting, variable current directions). The largest drawback to this approach is that it is 
extremely manpower intensive as separate crews are required for each component 


undertaken at the same time. 


If surveys of the above type are properly undertaken, they can be augmented with 
small-scale surveys (fixed grid and effluent monitoring) which would provide seasonal 
information on a year-round basis. It would be helpful to have winter data for parameters 
such as ammonia, which nitrifies at a much slower rate and therefore has greater capacity 
to exceed PWQO, despite the lower unionized fraction of total ammonia at lower 


temperatures. 


Given practical data limitations, a combination of plume tracking (including on-board 
measurement of several key parameters) with one or more modelling techniques is likely 


the most useful approach to setting effluent limits with a reasonable level of confidence. 


TABLE 5.1: UTILITY OF VARIOUS TECHNIQUES FOR ESTABLISHING 


EFFLUENT LIMITS 


Level of Effort 
Expended in 
This Study 


Usefulness in 
Technique Limit Derivation 


EFFLUENT CHEMISTRY: 


Sequential sampling Essential 
Daily sampling Essential 
Seasonal sampling Supportive 


AMBIENT WATER CHEMISTRY: 


Sequential sampling at 


outfall Highly supportive 
Plume tracking Essential 
Fixed-grid sampling Supportive 


PHYSICAL MONITORING: 
Drogue tracking Highly supportive 
Current metering Supportive 


Hydrodynamic modelling Supportive 


More than 
Adequate 


Less than 


Less than 
Less than 


More than 


Less than 
Adequate 


Adequate 


adequate 


adequate 


adequate 
adequate 


adequate 


adequate 


15a Oe 


i q 2 ‘ 
a ."@ 4 


t @gunra¥odat sli aay 10.1 Ts 
rv rf MET a 


woe 


ey, 


et = ee a ee» 


: var’ 


1? See 


FAR aeery || oLapy 2 
g Fiat Va 


‘Gon 
gitlcbeet 


& t =e 
* 


ae 
unit iqmas sis dex! 


6.0 CONCLUSIONS AND RECOMMENDATIONS 


6.1 Conclusions 


Calculations based on two models (dilution model and fate model), as well as on field- 
measured plume tracking data, have been used to estimate the extent of loading reductions 
required to meet Provincial Water Quality Objectives at various locations. It should be 
noted that these are only model estimates; the final model for loading reduction 
calculations has not been selected, nor has the allowable distance from the outfall been 


finalized. 


ih. It is important that site-specific water quality based effluent limits be developed in 
accordance with the MISA water quality track. This is evident due to difference in 


effluent reductions required for the present and proposed outfalls. 


2 Based on dilution model results (more conservative than the fate model), it is 
evident that dechlorination of the final effluent will be necessary in order to meet 
Provincial Water Quality Objectives (PWQO) in the receiving water. 


e At the existing outfall location, in order to meet PWQO at the end of the 
pipe 95% of the time, residual chlorine levels in the final effluent must be 
reduced by 99.6% to 5 ug/L. In order to meet PWQO at the edge of a 500m 
mixing zone 95% of the time, final effluent chlorine levels must be reduced 
by 99% to 16 ug/L. 


° At the proposed offshore outfall location, in order to meet PWOQO at the end 
of the pipe 95% of the time, chlorine levels in the effluent must be reduced 
by 97% to 47 ng/L. In order to meet PWQO at the edge of a 500m mixing 
zone 95% of the time, chlorine levels in the effluent must be reduced by 94% 
to 86 pg/L. 


Ammonia removal in the final effluent will be necessary to meet PWQO in the 


receiving water under summer conditions. 


e At the existing outfall, in order to meet PWQO at the end of the pipe 95% 
of the time, ammonia levels in the final effluent must be reduced by 88 to 
93% to 2.3 - 4.3 mg/L. Similarly, for the edge of a 500m mixing zone, the 
levels must be reduced by 79 to 87% to 4.3-6.9 mg/L. 


° At the proposed offshore location, in order to meet PWQO at the end of the 
pipe 95% of the time, ammonia levels in the final effluent must be reduced 
by 28% to 23.5 mg/L. However, no reduction in effluent ammonia levels is 
required to meet PWQO 95% of the time at the edge of a 500m mixing zone. 


Enhanced phosphorus removal (tertiary treatment) of the final effluent will be 
necessary to meet the Provincial Water Quality Guidelines (PWQG) for prevention 
of nuisance algal growth. , 

6 At the existing outfall location, in order to meet the PWQG at the end of the 
pipe 95% of the time, final effluent total phosphorus levels must be reduced 
by 96-98% to 24-46 ug/L. Similarly for the edge of a 500m mixing zone, the 
level must be reduced by 95% to 68 ug/L. Based on dilution model 
predictions, distances of 17 and 35 km are required for dilution to PWQG of 
effluent concentrations of 0.5 and 1.0 mg/L, respectively. 


e At the proposed offshore outfall location, in order to meet PWQG 95% of 
the time at the end of the pipe, final effluent total phosphorus levels must be 


reduced by 80% to 245 ug/L. Similarly for a 500m mixing zone, the levels 
must be reduced by 65% to 440 ug/L. Based on dilution model predictions, 
distances of 0.7 and 3.6 km are required for dilution to PWQG of effluent 
concentrations of 0.5 and 1.0 mg/L, respectively. 


6 f——2 


At the existing outfall location, loading reductions of heavy metals and organic 
contaminants will be necessary in order to meet PWQO in the receiving water. 
However, improved dilution characteristics at the proposed offshore outfall location 


will allow PWQO for these parameters to be achieved after initial dilution. 


° In order to meet PWQO for copper at the end of the pipe 95% of the time 
at the existing outfall location, the final effluent level must be reduced by 65- 
75% to 10-15 pg/L. Similarly, for the edge of a 500m mixing zone, the 
required reduction is 43-58% to 18-24 ug/L. Required reductions for other 
metals and most organics such as lindane are lower and should be achieved 


if the above reductions are met. 


e In the case of PCBs, loading reductions of 93 to 95% are estimated for zero 
and 500 m mixing zones, based on average observed background levels. 
However, the highest background PCB levels are similar to the PWQO of 1 
ng/L; hence required loading reductions may be as high as 100% for either 
outfall location. 


Both plume tracking and modelling studies show that higher current speeds result 
in reduced longitudinal mixing of the effluent plume. Under such conditions, 
increased surface conductivity and ammonia concentrations were observed at 
distances up to 6 km from the Main STP outfall. 


The following general conclusions summarize the most important water quality 
characteristics of the Main STP effluent and Lake Ontario in the vicinity of the 
outfall. 


a. Ammonia and chlorine approach the thresholds of acute toxicity to aquatic 


biota at the Main STP outfall, and are most likely responsible for any acute 


toxicity observed (Flood, 1990, Neville, 1990). Acute toxicity is not likely 
beyond a 500 m radius of the outfall. 


Total ammonia averages approximately 7 mg/L at the outfall (near the 
surface), and the PWQO for unionized ammonia is exceeded in approximately 
75% of samples. Ammonia concentration declines rapidly to about 0.1 mg/L 
(20% exceedence) at a distance of 3 km south of the outfall (background is 
less than 0.04 mg/L). 


Free chlorine averages approximately 330 pg/L at the outfall (near the 
surface) from June through September when chlorination occurs in the Main 
STP. The PWOQO for free chlorine is exceeded in approximately 70% of 
samples over this period. The concentration declines rapidly to about 20 ug/L 
(still 70% exceedence) 3 km south of the outfall (background is zero). 


Total phosphorus averages approximately 0.2 mg/L at the outfall (near the 
surface), and the PWQO is almost always exceeded here. Phosphorus declines 
to about 0.03 mg/L (30% exceedence) 3 km south of the outfall (background 
is about 0.02 mg/L and 15% exceedence). 


A slight pH depression at the outfall (pH 7.3 compared to 8.0 and above at 
surrounding stations) will slightly favour ionization of ammonia thereby 


suppressing ammonia toxicity to some extent. 


Heavy metals exceeding PWQOs include iron (77%), copper (38%) and 
cadmium (31%). Metal concentrations at the outfall near the surface average 
0.47 mg/L, 0.006 mg/L and 0.0005 mg/L, respectively, approaching 
background within 3 km south of the outfall. 


Organochlorine pesticides lindane, heptachlor and dieldrin occasionally 
exceed PWOOs at the outfall; however, average concentrations do not 


approach even chronic toxicity thresholds. 


Main STP effluent generally rises rapidly to the surface, resulting in a surface 
maximum for most contaminants. This pattern persists at stations within 
about 1,000 m of the outfall. | 


Spatial pattern analysis of the Toronto Waterfront indicates a strong water 
quality impact at the Main STP outfall, moderate impact at the Humber Bay 
STP outfall and Coatsworth Cut, and lesser impact zones near the Don River 
and Rouge River mouths and in Toronto Harbour. Long-term temporal 
trends at the Main STP outfall and Don River mouth also suggest ongoing 


degradation of water quality. 


Comparison of 1986 and 1987 water quality suggests slightly higher 
concentrations of most metals near the outfall in 1987; however, only 
manganese and nickel increases were statistically significant. Free chlorine 
discharge was nominally lower in 1987, though the difference was not 
significant, while ammonia and phosphorus concentrations at the outfall were 


similar between years. 


Potential surrogates for trace metals include suspended solids and total 
phosphorus in Main STP effluent, and conductivity and suspended solids in 
receiving water. In the Main STP effluent, the data collected were inadequate 
to define potential surrogates for organic compounds measured at ultra-trace 
levels. In receiving water, total phosphorus was found to be a potential 
surrogate for many of the ultra-trace organic parameters during June and 


August 1987. Total Kjeldahl nitrogen in June, and conductivity, total filtered 


residue and suspended solids in August also showed some usefulness as 


surrogates. No potential surrogates were found in October 1987. 


i Total ammonia in Main STP effluent shows a consistent pattern of increase 
over the course of the day from about 10:00 a.m. to 5:00 p.m. Phosphorus 
is more tightly controlled, though extremes tend to occur later in the day. 
Flow may either increase or decrease over the course of the day. Both 


ammonia and flow exhibit large day-to-day variations. 


m. Long-term prediction of daily average outfall concentrations may be possible 
with some success based on daily average concentrations in the Main STP 
effluent. However, short-term prediction based on instantaneous effluent 
measurements will require a hydrodynamic model which explicitly accounts 
for receiving water conditions. Existing data are inadequate for development 


of these predictive relationships. 
6.2 Recommendations 


Environmental survey design and data handling strategies should be derived from study 
objectives. Pilot study objectives appear to fall into four categories: 


it. Description and prediction of short-term (daily) temporal patterns in receiving water 
quality near point sources. 
2. Description and prediction of long-term (seasonal) temporal patterns in receiving 


water quality near point sources. 


ah Description of spatial patterns over the waterfront (and correspondence with known 
discharges). 
4, Compliance evaluation of effluent discharges and adjacent receiving waters with 


respect to effluent and receiving water objectives. 


The data needs, optimum sampling designs and data handling procedures differ for each 


of these objectives. Any design must consider that, unlike rivers where the downstream 


direction is simply defined, for a lake discharge, downstream’ varies in time and may not 


be definable at some times. There are also transition periods when downstream is changing 


in direction. Furthermore, unlike rivers, the effects of density stratification for lake 


discharges can persist for considerable distances from the outfall. The following design 


principles, data gathering and data handling procedures should be applied: 


la. 


1b. 


lc. 


For short-term (i.e. daily cycle) description and predictive model development, the 
sampling interval should be short (on the order of hours), both in the source effluent 
and the adjacent receiving water. The effluent and receiving water sampling times 
should match as closely as possible. A specific design (i.e. random, stratified 
random, systematic) should be chosen. Gilbert (1987) discusses the advantages and 
disadvantages of each design. 


Because conductivity has been found to be a good surrogate parameter for metals 
and some ultra-trace organics in the receiving water, recording conductivity meters 
could be mounted on floats and located near the outfall at various distances from 
the outfall, as well as in the effluent discharge manhole. This would provide 
complete short-term data on the variabilities in the discharge and the receiving 
water, and ensure a perfectly matched ’real time’ record. 


In order to provide better identification of surrogate parameters for toxic organics 
in the effluent, a more intensive sampling design involving properly matched samples 
for all parameters is required. This should include as large a range of effluent 
conditions as possible, and include improved techniques for handling samples high 
in suspended solids in order that all samples be properly representative of effluent 
conditions. 


1d. 


le. 


Empirical velocity data for the period of sampling would greatly enhance the ability 
to develop predictive models. This must be obtained by releasing a cluster of 


drogues and following the position of each drogue over time. 


Sampling periods need only cover a few days, but each period should represent a 
particular set of well-documented current, weather and wind conditions. With 
sufficient data, dynamic models (such as the Rand Model) can be used for receiving 


water predictions based on time of day and prevailing conditions. 


For long-term description and predictive model development, the primary sampling 
units are days. Again, a specific sampling design should be chosen (see Gilbert, 
1987). Steps should be taken to limit short-term (i.e. time of day) temporal sampling 
bias between days, especially near a source with a known or suspected daily cycle. 
Daily (24-h) composite effluent sampling, supported by more frequent (hourly to 3 
h) grab sampling will help to define the temporal cycle and its probable error limit. 


The same sampling schedule should be used for source and adjacent receiving water. 


For description of broad spatial patterns, every effort should be made to limit 
temporal sampling bias between regions. All stations should be sampled on each 


survey, preferably in a randomized sequence. 


The vertical distribution of samples with depth should follow a specific sampling 
design. The random design is least efficient. A stratified random design is more 
efficient if strata are easily recognized (e.g., samples identified as in-plume, above, 
or below), otherwise a systematic design (fixed depths) is preferred. Arbitrary 
pooling of depths for statistical summary is unwise when depth trends are known to 


exist in receiving water, and no particular sampling design has been followed. 


Any descriptive statistical inference (e.g. mean or variance estimation) requires that 


the sampling units, population boundaries (ie. spatial and temporal domain) and 


sampling design be defined. For example, replicate analyses on a single sample are 
not primary sampling units, but represent sub-units in a two-stage sampling design. 
To use them as primary units, when replication is unequal among the primary units 
of the population, will bias both the mean and variance estimates (Gilbert, 1987). 
Likewise, field replicates (taken at essentially the same spatial point and time) are 


not primary sampling units in a population of broader temporal or spatial domain. 


Preferred methods of handling censored (i.e. less than detection) data in mean and 
variance estimations will depend on the study objectives. For compliance purposes 
where a mean is compared to a water quality guideline, a conservative assumption 
of ’face value’ concentration in each sample is probably preferable. For descriptive 
purposes where a most likely estimate is desired, Cohen’s (1961) maximum likelihood 
method is probably preferable. Each element of a water quality summation should 
be related to one or the other (not both) of these purposes. 
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APPENDIX 1 


WATER QUALITY STATISTICAL SUMMARIES 
FOR TOXICITY TEST PERIODS 
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CONVENTIONALS AND HEAVY METALS 
DURING TOXICITY TESTS 
1987 
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STATION 1419 WATER QUALITY AT 0.5 and 
6 M DEPTHS, 15-19 JUNE 1987 


PAR? N MEAN MEDIAN SD MIN MAX DL? PWQO? N* 
eee 


Time (hours) 71120 9.9 1.9 9.3 13.8 
Temp. (°C) 6 12.9 12.8 0.4 12.4 13.5 
Chlorine (ug/L) 6 513 510 175 320 700 0 6 6 
Ammonia 6 5.46 5.72 2.97 0.56 8.30 0 3 4 
Lab Ammonia 0 = ~ a = = - - - 
pH (units) 5 / ier 4 7.3 0.4 6.6 7.4 = = = 
Phosphorus 7 0.14 0.11 0.11 0.08 0.39 0 7 7 
Cadmium (ug/L) 7 1.45 0.15 2.91 0.15 13.0 5 2 7 
Copper 7 0.005 0.004 0.004 0.001 0.013 O 1 ve 
Iron 7 0.284 0.300 0.093 0.170 0.410 oO 3 7 
Mercury (ug/L) - 7 0.015 0.010 0.009 0.005 0.030 1 0 7 
Zinc 7 0.006 0.007 0.004 0.001 0.011 2 0 7 
Turbidity (ftu) 6 3/3! 3.20 1.22 1.76 5.40 
6.0 m 
Time (hours) 6 10.8 10.1 1.7 oe On Last 
Temp (°C) 58.5 8.4 12 6.8 10.1 
Chlorine (ug/L) 4 0 0 0 0 0 4 4 
Ammonia 5 0.29 0.29 0.20 0.04 0.59 0 0 2 
Lab Ammonia 0 = = = = = = - = 
pH (units) 2 7.8 7.8 0.2 ee. 7.9 = ~ = 
Phosphorus 5 0.04 0.04 0.02 0.02 0.07 0 4 5 
Cadmium (yg/L) 5 2.36 0.15 4.43 0.15 9.00 4 1 5 
Copper 5 0.003 0.001 0.004 0.001 0.010 1 Hock aby 
Iron 5 0.182 0.150 0.110 0.082 0.350 0 3 5 
Mercury (pg/L) 5 0.005 0.005 0 0.005 0.005 1 0 5 
Zinc 5 0.004 0.002 0.006 0.001 0.014 1 0 5 
Turbidity (ftu) 5 2.60 2.20 0.75 2.00 3.60 
eee eee ee ee ee ee eT 


1 Parameter units are mg/L unless otherwise indicated. 

2 Number of samples below detection limit. 

* Number of samples above Provincial Water Quality Objectives (PWQO). 
* Number of samples compared to PWQO. 


2 0 KN 


STATION 1419 WATER QUALITY AT 0.5 and 
6 M DEPTHS, 16 JUNE-08 JULY 1987 


PAR! N MEAN MEDIAN SD MIN MAX DL2 PWQO> N‘* 
0.5m 

Time (hours) 6 ele 2 10.4 220 9.3 £13.58 

Temp. (°C) S 13.0 12.9 O74 & 12.488 13.5 

Chlorine (ug/L) 5 496 420 190 320 700 0 5 5 
Ammonia 5 4.91 4.50 2.96 0.56 | 8.30 Z 0 3 
Lab Ammonia 0 = = = oo = = = = 
pH (units) 4 = ~ 3 = = = > = 
Phosphorus 6 0.10 0.11 0.02 0.08 0.13 6 0 6 
Cadmium 6 2 4 6 
Copper 6 0.005 0.004 0.004 0.001 0.013 1 0 6 
Iron 6 0.277 0.265 0.100 0.170 0.410 2 0 6 
Mercury (ug/L) 6 0.013 0.010 0.006 0.005 0.020 0 1 6 
Zinc 6 0.006 0.006 0.004 0.001 0.011 0O 2 6 
Turbidity (ftu) 5 3e23 3.20 13 eee) 4.0 

6.0 m 

Time (hours) 5 11.0 10.1 1.8 9.5 14.1 

Temp. (°C) 4 8.9 8.7 0.8 8.3 10.1 

Chlorine (ug/L) 4 0 0 0 0 0 4 0 é 
Ammonia 4 0.36 0.31 0.16 0.23 0.59 0 0 2 
Lab Ammonia 0 = = = = > = - ~ 
pH (units) 2 = = = = = = S = 
Phosphorus 5 0.04 0.04 0.02 0.02 0.07 4 0 5 
Cadmium 4 0 4 4 
Copper 4 0.001 0.001 0 0.001 0.001 0 1 “ 
Iron 4 0.191 0.165 0.125 0.082 0.350 1 0 4 
Mercury (pg/L) 5 0.005 0.005 0 0.005 0.005 0 5 5 
Zinc 4 0.002 0.002 0.001 0.001 0.002 0 1 4A 
Turbidity (ftu) 5 2.60 2.20 0.75 2.00 3.60 


Parameter units are mg/L unless otherwise indicated. 
Number of samples below detection limit. 
Number of samples above Provincial Water Quality Objectives (PWQO). 
Number of samples compared to PWQO. 


APPENDIX 2 


Water Quality Statistical Summaries 
for Each Survey Year 
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ANNUAL MEAN 
HEAVY METALS 
1986 


i 


MAIN STP (STATION 22) FINAL EFFLUENT 
QUALITY, ANNUAL MEAN, 1986 


PAR? N MEAN MEDIAN SO MIN MAX DL? 
Time (hours) 16 12.7 13.0 13 10.0 14.5 

Cadmium (pg/L) 15 0.793 0.561 0.628 0.150 2.34 2 
Chromium 15 0.025 0.023 0.015 0.012 0.059 0 
Copper 15 0.027 0.015 0.024 0.006 0.081 0 
Iron 15 0.844 0.615 0.585 0.323 2.292 0 
Mercury (pg/L) 15 0.075 0.062 0.053 0.020 0.186 0 
Manganese 15 0.054 0.053 0.005 0.049 0.068 0 
Nickel 11S) 0.027 0.025 0.008 0.015 0.039 0 
Zinc 15 0.070 0.049 0.050 0.018 0.195 0 


1 Parameter units are mg/L unless otherwise indicated. 
2 Number of samples below detection limit. 
3 Number of samples above Industrial Waste Discharge Guidelines (IWDG). 


oOo oO +> 


STATION 1419 WATER QUALITY (ug/L) AT SURFACE AND 
BOTTOM, ANNUAL MEAN, 1986 


(a) Surface (0 to 1 m water depth) 


PAR. MEAN MEDIAN OEY MIN MAX N DL PwQO0 
Cd 0.18 0.15 0.05 Onis 0.27 10 0 2 
Cr 5.0 5.3 ces 1.8 SS 10 0 0 
Cu 4.4 KN) 2.4 15 8.9 10 0 3 
Fe 281. 244. 213. 66. TES 10 0 4 
Hg 0.02 0.02 0.01 <0.01 0.04 9 1 1 
Mn Wee 16. 6. 6.8 28.5 10 0 = 
Ni 7. 7. Ze 4.3 Liz 10 0 0 
Zn Le 11. 7. 3.8 Zen 10 0 0 

(b) Bottom or depth of turbidity maximum (4 to 6 m water depth) 
STD 

VAR. MEAN MEDIAN DEV MIN MAX N DL PwQO 
Cd 0.30 0.15 0.44 Onis 1.46 9 0 1 
Cr 3.0 ore Zies 0.8 8.2 9 0 0 
Cu 3.0 rye 2.0 1.0 * fa’ 9 0 1 
Fe 272. 89. 345. 46. 1037. 9 0 2 
Hg 0.01 0.01 0.00 <0.01 0.02 8 3 0 
Mn 9. 4. oF 2s 25. 9 0 = 
Ni 4. 32 4. 1, 15% 9 0 0 
Zn Wn 5 6. 2 Zin 9 0 0 

Note: DL = Number of samples below detection limit. 
PwQO = Number of samples above Provincial Water Quality 


Objective (PWQO). 


STATION 1536 WATER QUALITY (g/L) AT SURFACE AND 
BOTTOM, ANNUAL MEAN, 1986 


(a) Surface (0 to 1 m water depth) 


PAR. MEAN — MEDIAN DEV. «MIN MAX ON DL PwQO 
Cd <0.3 <0.3 8 8 0 
Cr <1.0 1.0 8 5 0 
Cu 1.1 1.0 0.5 <1.0 2.0 8 2 0 
Fe 42. 44, 22. 14. 80. 8 0 0 
Hg <0.01 0.06 8 5 0 
Mn 2.6 2.0 1.6 1.0 5.0 8 0 = 
Ni ~<2.0 2.0 8 7 0 
Zn 3.2 2.0 3.6 <1.0 11.0 8 2 0 

(b) Bottom or depth of turbidity maximum (4 to 20 m water depth) 
STD 

PAR. MEAN MEDIAN DEV MIN MAX N DL PWQO 
Cd <0.3 <0.3 7 7 0 
Cr <1.0 2.0 7 4 
Cu 0.9 1.0 0.2  ° °&<1.0 1.0 7 1 0 
Fe Ben 31. 10. ~ 13. 47. 7 0 0 
Hg <0.01 <0.01 7 7 0 
Mn Ned 2.0 0.8 1.0 3.0 7 0 = 
Ni <2.0 4.0 7 6 
Zn 2.1 1.0 1.6 <1.0 5.0 7 2 0 

Note: DL = Number of samples below detection limit. 


PwQO 


Number of samples above Provincial Water Quality 
Objective (PWQO). 


STATION 2029 WATER QUALITY (pg/L) AT SURFACE AND 
BOTTOM, ANNUAL MEAN, 1986 


(a) Surface (0 to 1 m water depth) 


PAR. MEAN MEDIAN DEV. MIN MAX oN DL PWQO 
Cd <0.3 <0.3 7 7 0 
Cr <1.0 1.0 7 4 0 
Cu 0.9 1.0 02 “<Aip 10.6 7 2 0 
Fe 21.0 14.0 19.3 35) 50.08R 7 0 0 
Hg <0.01 <0.01 6 6 0 
Mn 1.3 1.0 0.8 0.8 308 7 1 - 
Ni <2.0 DG 7 
Zn 1.4 1.0 0.8 0.5 28 7 2 0 


SSS Ss ss SSS 


(b) Bottom or depth of turbidity maximum (6 to 15 m water depth) 
a I ae he a Ae)” 55s 


PAR. MEAN — MEDIAN BEV MIN MAX N DL  PWQO 
a ane NR ee 1c Ne a tet ate A oo 
Cd <0.3 <0.3 6 6 0 
Cr <1.0 1.0 é 5 0 
Cu 1.0 1.0 0.0 1.0 1.0 6 0 0 
Fe 25. 20 nL OF 6 60. 6 0 0 
Hg <0.01 <0.01 5 5 0 
Mn 12 1.0 0.4 1.0 2.0 6 0 - 

Ni <2.0 2.0 6 5 0 
Zn 2.4 mT 1.2 1.0 4.0 6 0 0 


Note: DL = Number of samples below detection limit. 
PWQO = Number of samples above Provincial Water Quality 
Objective (PWQO). 
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Qual? 


MAIN STP (STATION 22) FINAL EFFLUENT 
QUALITY, ANNUAL MEAN, 1987 


PAR? N MEAN MEDIAN SD MIN MAX DL? IwDG? 
Time (hours) 25 10.6 10.7 1.8 8.3 f13h7 

Lab Ammonia 25 21.0 20.5 4.9 8.9 ans es 

pH (units) 24 7.4 7.4 0.4 6.6 8.1 

Phosphorus 25 0.83 0.34 AO) 0512 eins 0 2 
Nitrite/Nitrate 23 0.85 0.25 2209" 0705 #210. 2 

TKN 23 30.3 2550 LUG me liiccemleles 
Conductivity 25 924 938 70 705 1,027 

(umhos/cm) 

DOC 23 14.0 14.3 2.6 9-8" 718.5 

Suspended Solids 25 27.5 8.7 72.9 371 362:5 0 4 
Fecal Coliforms 24 3.48 2.94 1704 1.37) 6717 


(log #/100 mL) 


1 Parameter units are mg/L unless otherwise indicated. 


2 Number of samples below detection limit. 
3 Number of samples above Industrial Waste Discharge Guidelines (IWDG). 
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APPENDIX 3 


CONCENTRATION MAPS AND 
DILUTION FACTORS FOR 1986 
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TABLE A3.1: DILUTION FACTORS? FOR DIFFERENT CHEMICALS AT SELECTED 
STATIONS IN 1986 


a ee | eee, 
PARAMETER — 1419 _ 1431 1432 1433 1434 


i 


Field Ammonia: 
o 0.5 m depth 0.33 0.12 0.03 0.09 0.07 
o .6.0 m depth 0.08 0.005 0 0 0.005 


Lab Ammonia: 


o 0.5 m depth O31 0.11 —0.03 0.08 0.07 
o 6.0 m depth 0.06 0.01 0 0.01 0.005 
Phosphorus: , 

o 0.5 m depth 0.22 0.04 0.01 0.05 0.06 
o 6.0 m depth 0 0 0 0 0 


* Dilution Factor = (Concentration-Background)/(Source-Background) 
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APPENDIX 4 


COMPARISON OF METHODS FOR 
CENSORED DATA HANDLING 
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APPENDIX 4: COMPARISON OF METHODS FOR CENSORED DATA HANDLING 


A number of methods have been used for estimation of mean and variance 
from left- censored data sets. A data set is left-censored when 
contaminant concentrations in some of the samples are lower than the 
detection limit of the analytical instrument used to measure them. 

Gilbert (1987) describes most of the methods which have been used. Helsel 
(1986) compares the performance of some of these methods on artificial 
data sets from different known distributions. We have compared four 
methods using Main STP effluent data. 


Three of the methods compared here involve replacement of the 'non-detect' 
observation with an estimate of the true value: 


1. DL -- replacement by face value of detection limit. 
2. 4#DL - replacement by half the face value. 
3. LR - linear regression of measured values on rankits. 


The fourth method provides a direct estimate of true mean and variance 
based on knowledge of the total number of observations and the proportion 
below the detection limit: 


4. ML —- maximum likelihood estimation. 


Methods 1 and 2 are arbitrary, while 3 and 4 assume normality and strive 
to optimize the solution to the problem. A consensus seems to be emerging 
that the maximum likelihood method (Cohen, 1961) is superior, if 
assumptions are satisfied; however, it is quite sensitive to non-normality. 
Helsel (1986) notes that the regression method is more robust. Every 
effort should be made to normalize the data by transformation prior to 
application of either method. A log transformation is generally suitable 
for log-normal data. 


Table A4.1 shows the estimate of both geometric and arithmetic means for 
six chemical parameters in Main STP effluent (1987 data) using the four 
different methods. Original distributions, and distributions after LR 
replacement, are illustrated in Figures A4.1 through A4.6. The table 
illustrates several general points. 


Firstly, the geometric mean is a biased estimator of the u parameter of 
the log-normal distribution. There is a low bias which can be quantified 
as follows (Gilbert, 1987): 


u | 296 oe ov) | 


where: E(GM) expected value of geometric mean, 
u = true mean, and 


E(GM) 


7 y? = variance of In transformed data. 
The term in square brackets is the bias factor. 


Secondly, the methods are in reasonably close agreement for the degree of 
censoring encountered in the Main STP effluent. However, the initial 
assumption of normality or log normality is critical. For example, when 
log-normality is assumed for ammonia, which is actually closer to normal 
in distribution, the maximum likelihood method produces a u estimate out 
of line with other methods, and with a very large variance. 


This last point raises some question as to the quantity of data needed to 
produce a reliable estimate. Helsel (1986) recommends that little can be 
done with data that is more than 80% censored with n = 10. With the 
quantity of data available at most individual receiving water stations, 
after depth partitioning, it will be very difficult to decide with any 
certainty on the form of the underlying distribution. We suggest that at 
least 20 measured observations should be available to make this 
determination. 


TABLE A4.1: COMPARISON OF METHODS FOR CENSORED DATA HANDLING OF 
MAIN STP EFFLUENT CONCENTRATIONS 


ee 


Geometric Mean2 Log-normal$ 
PARAMETER h? OL 40L LR ML u V 


Zinc (mg/L) 0.08 0.026 0.026 0.027 0.026 0.05 0.01 
Copper (mg/L) 0.12 0.017 0.016 0.016 0.017 0.02 0.003 
NO,/NO, (mg/L) 0.13 0.24 0.24 0°23) 0:24 0.77) 5.54 
Ammonia* (mg/L) 0.02 D/P OSL 78) 205081820 26.7 858 
Dichloromethane (ug/L) 0.11 4.0 4.0 4.6 4.5 5.5 14.8 
Chloroform (pg/L) 0.13 2.1 2.04 2.2 are 2.8 a0 


Arithmetic Mean2 Normal? 
PARAMETER h DL 40L LR ML u V 


Se SSS 


Zinc (mg/L) 0.08 0.06 0.06 - 0.05 0.05 0.01 
Copper (mg/L) 0.12 0.03 0.03 - 0.03 0.03 0.003 
NO,/NO, (mg/L) 0513 0.82 0.82 Hs 0.75 0.75 4.79 
Ammonia* (mg/L) 0.02 20.4 20.4 20.6 20.4 20.4 32.2 
Dichloromethane (pg/L) 0.11 5.4 aes) - 5.3 Spe.) bg) 
Chloroform (pg/L) 0.13 2.8 2.7 = 2.6 2.6 3.9 


+h = proportion of distribution below detection limit. | 


2d Geometric mean computed from mean of logs, arithmetic mean from 
untransformed data, with 'less thans' treated as: 


DL = face value of detection limit. 

40L = one-half face value of detection limit. 

LR = Jinear regression estimate of true value. 

ML = input to maximum likelihood estimate of mean. 


3 u and V are maximum likelihood estimates of true mean and variance for 
log-normal or normal distribtuion. 


* Ammonia is approximately normal, not log-normal in distribution. 


DISTRIBUTION OF ZINC CONCENTRATION (mg/L) IN MAIN STP 


FIGURE A4.1: 
EFFLUENT BEFORE AND AFTER LINEAR REGRESSION 
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FIGURE A4.2: DISTRIBUTION OF COPPER CONCENTRATION (mg/L) IN MAIN STP 
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EFFLUENT BEFORE AND AFTER LINEAR REGRESSION 
REPLACEMENT OF DATA BELOW THE DETECTION LIMIT 
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FIGURE A4.3: 


Frequency 


Frequency 


DISTRIBUTION OF NO5/NO3 CONCENTRATION (mg/L) IN MAIN 
STP EFFLUENT BEFORE AND AFTER LINEAR REGRESSION 
REPLACEMENT OF DATA BELOW THE DETECTION LIMIT 
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DISTRIBUTION OF AMMONIA CONCENTRATION (mg/L) IN MAIN 
STP EFFLUENT BEFORE AND AFTER LINEAR REGRESSION 
REPLACEMENT OF DATA BELOW THE DETECTION LIMIT 


FIGURE A4.4: 
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DISTRIBUTION OF DICHLOROMETHANE CONCENTRATION (ug/L) 
IN MAIN STP EFFLUENT BEFORE AND AFTER LINEAR 
REGRESSION REPLACEMENT OF DATA BELOW THE DETECTION 


FIGURE A4.5: 
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FIGURE A4.6: DISTRIBUTION OF CHLOROFORM CONCENTRATION (ug/L) IN 
MAIN STP EFFLUENT BEFORE AND AFTER LINEAR REGRESSION 
REPLACEMENT OF DATA BELOW THE DETECTION LIMIT 
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APPENDIX 5 


PLUME TRACKING RESULTS 


June 16 - 17, 1987 
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